Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,483,675 T→C S208S (TCT→TCC)  unknown → Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,483,6750TC100.0% 235.9 / NA 85S208S (TCT→TCC) unknownType II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/47);  total (38/47)

CGCCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCAGTAAGAATAAAAAAACATTGCAAGAGCTTGTAAATTTTGGTTCTATGACTCGCGAAATGGCGAATTTCTTAATTATTGCTGCCAGATCCCGGGTAAATATAATTGTGTCAGGTGGTACAGGATCTGGTAAAACAACCTTACTTAATGCGCTATCTAACTATATTTCCCATACCGAACGTGTCATCACACTAGAAGATACTGCCGAGCTGCGTTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAA  >  1/1483433‑1483916
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
cgtcGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaac                                                                                                                                                                                                                                           <  1:618678/248‑1 (MQ=255)
cgtcGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaac                                                                                                                                                                                                                                           <  1:180538/248‑1 (MQ=255)
  tcGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaa                                                                                                                                                                                                                                         >  1:323857/2‑251 (MQ=255)
    gTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaa                                                                                                                                                                                                                                         >  1:273647/1‑249 (MQ=255)
              tgGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactctccttaatg                                                                                                                                                                                                                             <  1:814888/251‑12 (MQ=255)
              tgGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactctccttaatg                                                                                                                                                                                                                             <  1:819003/251‑12 (MQ=255)
              tgGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactctccttaatg                                                                                                                                                                                                                             >  1:519038/1‑240 (MQ=255)
                gTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATgc                                                                                                                                                                                                                            <  1:635676/250‑1 (MQ=255)
                   gTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTa                                                                                                                                                                                                                        >  1:601913/1‑251 (MQ=255)
                   gTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTa                                                                                                                                                                                                                        >  1:40588/1‑251 (MQ=255)
                           gTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                >  1:408684/1‑248 (MQ=255)
                            tGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTTGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                <  1:448071/250‑4 (MQ=255)
                            tGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                <  1:697782/250‑4 (MQ=255)
                            tGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                <  1:422294/250‑4 (MQ=255)
                              gACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATtata                                                                                                                                                                                                              <  1:2897/250‑1 (MQ=255)
                                     gACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCAt                                                                                                                                                                                                      >  1:249121/1‑251 (MQ=255)
                                       cTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATAc                                                                                                                                                                                                    <  1:235047/251‑1 (MQ=255)
                                         acccgATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                >  1:640523/5‑234 (MQ=255)
                                          cccgATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGa                                                                                                                                                                                                 >  1:761072/4‑251 (MQ=255)
                                            cgATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAAc                                                                                                                                                                                               <  1:268674/250‑1 (MQ=255)
                                             gATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCt                                                                                                                                                                                                                     >  1:50244/1‑228 (MQ=255)
                                              aTGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactctcctt                                                                                                                                                                                                                                 <  1:89988/215‑8 (MQ=255)
                                                      cGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATcac                                                                                                                                                                                     <  1:418083/251‑1 (MQ=255)
                                                        tttAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATcaca                                                                                                                                                                                    <  1:516466/250‑1 (MQ=255)
                                                              tgttgtAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTagaag                                                                                                                                                                             <  1:489297/251‑1 (MQ=255)
                                                                ttgtaATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                <  1:571432/214‑4 (MQ=255)
                                                                ttgtaATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGAt                                                                                                                                                                           <  1:544612/251‑1 (MQ=255)
                                                                  gtaaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATgc                                                                                                                                                                                                                            <  1:619802/197‑1 (MQ=255)
                                                                  gtaaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCt                                                                                                                                                                                                                     >  1:731078/4‑207 (MQ=255)
                                                                   taaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactctcctt                                                                                                                                                                                                                                 <  1:479978/192‑8 (MQ=255)
                                                                   taaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactc                                                                                                                                                                                                                                      <  1:616421/187‑3 (MQ=255)
                                                                   taaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATAc                                                                                                                                                                         <  1:401642/248‑1 (MQ=255)
                                                                   taaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATAc                                                                                                                                                                         <  1:521795/248‑1 (MQ=255)
                                                                                 gcgttAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTagaag                                                                                                                                                                             >  1:595869/4‑232 (MQ=255)
                                                                                 gcgttAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGttt                                                                                                                                                          >  1:437599/4‑251 (MQ=255)
                                                                                   gttAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAGCTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                >  1:429363/2‑192 (MQ=255)
                                                                                      aGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGaa                                                                                                                                                      >  1:883240/1‑250 (MQ=255)
                                                                                        aTGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAAc                                                                                                                                                     <  1:321550/249‑1 (MQ=255)
                                                                                        aTGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACaa                                                                                                                                                   >  1:144450/1‑250 (MQ=255)
                                                                                         tggtacaTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAAc                                                                                                                                                     >  1:652590/5‑248 (MQ=255)
                                                                                         tggtacaTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTAAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGa                                                                                                                                                       >  1:333926/5‑246 (MQ=255)
                                                                                           gtacaTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACaac                                                                                                                                                  >  1:277005/3‑247 (MQ=255)
                                                                                               aTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATAttt                                                                                                                                                                                                           >  1:377932/1‑188 (MQ=255)
                                                                                               aTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAAc                                                                                                                                                                                               >  1:49684/1‑200 (MQ=255)
                                                                                                tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                           <  1:843206/251‑1 (MQ=255)
                                                                                                tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                           <  1:290863/251‑1 (MQ=255)
                                                                                                tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                           <  1:450213/251‑1 (MQ=255)
                                                                                                tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                           <  1:207061/251‑1 (MQ=255)
                                                                                                       ccaTTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGc                                                                                                                                    <  1:508393/249‑1 (MQ=255)
                                                                                                                 gttCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATgc                                                                                                                                                                                                                            <  1:340412/152‑1 (MQ=255)
                                                                                                                 gttCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAAca                                                                                                                           >  1:384333/2‑249 (MQ=255)
                                                                                                                          gagtAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGcc                                                                                                                  >  1:628033/4‑250 (MQ=255)
                                                                                                                            gtAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCgg                                                                                                                >  1:46410/2‑250 (MQ=255)
                                                                                                                              aaaaaaaGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAAc                                                                                                                                                                                               >  1:161846/1‑169 (MQ=255)
                                                                                                                               aaaaaagcttGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGaaa                                                                                                                             <  1:626718/226‑1 (MQ=255)
                                                                                                                                                cTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGaggt                                                                                           >  1:579715/1‑248 (MQ=255)
                                                                                                                                                      aaTTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAaacaactctcctt                                                                                                                                                                                                                                 >  1:857036/1‑104 (MQ=255)
                                                                                                                                                       aTTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATg                                                                                    <  1:815454/251‑1 (MQ=255)
                                                                                                                                                        tttcggTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGc                                                                                   <  1:354902/247‑1 (MQ=255)
                                                                                                                                                             gTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                           >  1:560235/1‑190 (MQ=255)
                                                                                                                                                               tCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATAc                                                                                                                                                                                                    >  1:765229/1‑131 (MQ=255)
                                                                                                                                                               tCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgttg                                                                                                                                         <  1:656595/190‑1 (MQ=255)
                                                                                                                                                                          cgataTGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAAt                                                                  <  1:554324/246‑1 (MQ=255)
                                                                                                                                                                                   aaaTTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCactacgta                                                        <  1:96476/250‑8 (MQ=255)
                                                                                                                                                                                     aTTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                      <  1:241731/251‑1 (MQ=255)
                                                                                                                                                                                       tttttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                    <  1:608606/248‑1 (MQ=255)
                                                                                                                                                                                       tttttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                    <  1:389156/248‑1 (MQ=255)
                                                                                                                                                                                          ttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTcc                                                 <  1:480766/251‑1 (MQ=255)
                                                                                                                                                                                                      gcCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                      >  1:311426/1‑234 (MQ=255)
                                                                                                                                                                                                      gcCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                    <  1:818038/236‑1 (MQ=255)
                                                                                                                                                                                                      gcCAGATCTCGAGTAAATATAATTGTGTCAGGTGGCACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattat                                     >  1:63340/1‑251 (MQ=255)
                                                                                                                                                                                                               cgagTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTAc                                                                                                                                      <  1:409918/142‑1 (MQ=255)
                                                                                                                                                                                                                 agTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGc                                                                                                                                    >  1:722098/2‑145 (MQ=255)
                                                                                                                                                                                                                    aaaTATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAAc                                                                                                                            <  1:100438/150‑1 (MQ=255)
                                                                                                                                                                                                                    aaaTATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGcc                       <  1:493036/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaaTATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGc                        >  1:889650/1‑250 (MQ=255)
                                                                                                                                                                                                                      atatAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCgt                     >  1:644994/1‑251 (MQ=255)
                                                                                                                                                                                                                        ataATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCgtgg                   <  1:144208/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gtggtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCAtt       >  1:829725/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGc                                                                                   <  1:558827/171‑1 (MQ=255)
                                                                                                                                                                                                                                        gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                      >  1:794791/1‑200 (MQ=255)
                                                                                                                                                                                                                                        gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                    <  1:146158/202‑1 (MQ=255)
                                                                                                                                                                                                                                        gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattatcgtaggt                              >  1:361783/1‑217 (MQ=255)
                                                                                                                                                                                                                                        gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCag   >  1:334492/1‑250 (MQ=255)
                                                                                                                                                                                                                                         tACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCaga  <  1:494249/251‑3 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
CGCCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCAGTAAGAATAAAAAAACATTGCAAGAGCTTGTAAATTTTGGTTCTATGACTCGCGAAATGGCGAATTTCTTAATTATTGCTGCCAGATCCCGGGTAAATATAATTGTGTCAGGTGGTACAGGATCTGGTAAAACAACCTTACTTAATGCGCTATCTAACTATATTTCCCATACCGAACGTGTCATCACACTAGAAGATACTGCCGAGCTGCGTTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAA  >  1/1483433‑1483916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑