Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 321,223 C→A 59.9% intergenic (‑148/+65) ykgH ← / ← betA putative inner membrane protein/choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273321,2230CA59.9% 88.0 / 274.6 188intergenic (‑148/+65)ykgH/betAputative inner membrane protein/choline dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (47/28);  new base A (72/40);  total (119/69)
Fisher's exact test for biased strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 188 total aligned reads displayed.

AATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATCA  >  CP009273/321076‑321370
                                                                                                                                                   |                                                                                                                                                   
aataTATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAgg                                                                                                                                                   <  2:3270050/150‑1 (MQ=255)
           agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                        <  2:3511130/150‑1 (MQ=255)
           agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                        <  2:3511117/150‑1 (MQ=255)
           agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                        <  2:2022626/150‑1 (MQ=255)
             aaTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTgg                                                                                                                                      >  1:1206088/1‑150 (MQ=255)
              aTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGa                                                                                                                                     <  1:3013607/150‑1 (MQ=255)
              aTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGa                                                                                                                                     >  1:1176087/1‑150 (MQ=255)
              aTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGa                                                                                                                                     <  1:581708/150‑1 (MQ=255)
                aaaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAAc                                                                                                                                   >  2:3401105/1‑150 (MQ=255)
                  aaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACgg                                                                                                                                 >  2:3425885/1‑150 (MQ=255)
                   aTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGt                                                                                                                                >  1:1722887/1‑150 (MQ=255)
                    tCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTa                                                                                                                               >  1:794081/1‑150 (MQ=255)
                      tGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAAc                                                                                                                             <  1:2227875/150‑1 (MQ=255)
                           cAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCa                                                                                                                        >  2:2722669/1‑150 (MQ=255)
                           cAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATAACAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCa                                                                                                                        >  2:293718/1‑150 (MQ=255)
                            aGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCaa                                                                                                                       >  1:4646155/1‑150 (MQ=255)
                            aGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCaa                                                                                                                       >  1:1781296/1‑150 (MQ=255)
                              aTAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATc                                                                                                                     >  1:1735735/1‑150 (MQ=255)
                                aaGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTg                                                                                                                   >  1:4015893/1‑150 (MQ=255)
                                 aGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGa                                                                                                                  >  1:33041/1‑150 (MQ=255)
                                  ggTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGAt                                                                                                                 <  1:3747040/150‑1 (MQ=255)
                                     tGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCgg                                                                                                              >  1:1376166/1‑150 (MQ=255)
                                        aTTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTc                                                                                                           <  1:2396886/150‑1 (MQ=255)
                                          tAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCt                                                                                                         <  1:464297/150‑1 (MQ=255)
                                           aaTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTg                                                                                                        >  1:631898/1‑150 (MQ=255)
                                            aTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGc                                                                                                       >  1:4562438/1‑150 (MQ=255)
                                              aGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGt                                                                                                     >  1:2114185/1‑150 (MQ=255)
                                                tCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTa                                                                                                   >  2:4450320/1‑150 (MQ=255)
                                                  tGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGt                                                                                                 >  2:140187/1‑150 (MQ=255)
                                                    tATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTc                                                                                               >  1:2440093/1‑150 (MQ=255)
                                                     aTCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTca                                                                                              >  2:2575302/1‑150 (MQ=255)
                                                     aTCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTca                                                                                              <  1:684899/150‑1 (MQ=255)
                                                        aGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACAt                                                                                           <  2:4581923/150‑1 (MQ=255)
                                                           aaTGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCAc                                                                                        >  2:4641000/1‑150 (MQ=255)
                                                                tCGGTTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTc                                                                                   <  2:1614/150‑1 (MQ=255)
                                                                    gTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttt                                                                               >  2:4745003/1‑150 (MQ=255)
                                                                      tAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttttt                                                                             >  2:411598/1‑150 (MQ=255)
                                                                      tAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttttt                                                                             >  2:1395479/1‑150 (MQ=255)
                                                                        aaTATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATT‑GGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                                          <  2:3467965/150‑1 (MQ=255)
                                                                         ctatCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                                          <  2:566272/149‑1 (MQ=255)
                                                                         atatCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                                          <  2:3324656/150‑1 (MQ=255)
                                                                           atCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGct                                                                        <  1:4586602/150‑1 (MQ=255)
                                                                            tCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGctc                                                                       >  1:575168/1‑150 (MQ=255)
                                                                             cAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGctct                                                                      >  2:1438463/1‑150 (MQ=255)
                                                                             cAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGctct                                                                      >  1:862258/1‑150 (MQ=255)
                                                                              aGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCtctc                                                                     >  1:2920508/1‑150 (MQ=255)
                                                                                 aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                  <  1:3425885/150‑1 (MQ=255)
                                                                                 aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                  >  2:4724568/1‑150 (MQ=255)
                                                                                 aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                  >  2:3112973/1‑150 (MQ=255)
                                                                                  aaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCg                                                                 >  2:2322258/1‑150 (MQ=255)
                                                                                    aaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGc                                                               >  2:2177463/1‑150 (MQ=255)
                                                                                     aGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCa                                                              >  2:2420329/1‑150 (MQ=255)
                                                                                           catcatGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAt                                                        <  2:1608441/150‑1 (MQ=255)
                                                                                              catGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTg                                                     <  1:4233689/150‑1 (MQ=255)
                                                                                              catGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTg                                                     <  1:4233512/150‑1 (MQ=255)
                                                                                                tGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGcc                                                   >  1:3923452/1‑150 (MQ=255)
                                                                                                tGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGcc                                                   >  2:26468/1‑150 (MQ=255)
                                                                                                 gAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCca                                                  <  2:383228/150‑1 (MQ=255)
                                                                                                 gAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCca                                                  <  2:2318767/150‑1 (MQ=255)
                                                                                                         gATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATGTGCCACAAAatat                                          >  2:4321833/1‑150 (MQ=255)
                                                                                                         gATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAatat                                          >  2:1928298/1‑150 (MQ=255)
                                                                                                          atatGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATc                                         >  2:4821412/1‑150 (MQ=255)
                                                                                                            atGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATccc                                       >  2:2520424/1‑150 (MQ=255)
                                                                                                            atGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATccc                                       >  1:4867531/1‑150 (MQ=255)
                                                                                                            atGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATccc                                       >  1:4865979/1‑150 (MQ=255)
                                                                                                             tGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccg                                      <  1:358248/150‑1 (MQ=255)
                                                                                                             tGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccg                                      >  2:2713299/1‑150 (MQ=255)
                                                                                                              gAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgc                                     >  2:3955601/1‑150 (MQ=255)
                                                                                                               aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                    >  2:886358/1‑150 (MQ=255)
                                                                                                               aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                    >  1:3138330/1‑150 (MQ=255)
                                                                                                                aGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgccg                                   >  1:519928/1‑150 (MQ=255)
                                                                                                                 gCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGt                                  >  1:124511/1‑150 (MQ=255)
                                                                                                                   aTGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGc                                >  2:2796543/1‑150 (MQ=255)
                                                                                                                       gagTTACCTCAGTGTTTATATATGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATAGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctcct                            <  1:2308948/150‑1 (MQ=255)
                                                                                                                       gagTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctcct                            <  1:333144/150‑1 (MQ=255)
                                                                                                                           tACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGc                        >  2:2100780/1‑150 (MQ=255)
                                                                                                                            aCCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCa                       <  1:1318786/150‑1 (MQ=255)
                                                                                                                              cTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgc                     >  2:1188954/1‑150 (MQ=255)
                                                                                                                              cTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgc                     >  2:621184/1‑150 (MQ=255)
                                                                                                                              cTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgc                     >  1:1659464/1‑150 (MQ=255)
                                                                                                                                cAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgcgc                   >  2:1707961/1‑150 (MQ=255)
                                                                                                                                 aGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                  <  2:1745703/150‑1 (MQ=255)
                                                                                                                                    gtTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTcc               >  1:3215000/1‑150 (MQ=255)
                                                                                                                                      ttATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTg             <  1:2107381/150‑1 (MQ=255)
                                                                                                                                      ttATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTg             >  2:4160342/1‑150 (MQ=255)
                                                                                                                                      ttATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTg             >  2:3039947/1‑150 (MQ=255)
                                                                                                                                       tatataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGt            <  1:963513/150‑1 (MQ=255)
                                                                                                                                       tatataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGt            <  1:4078469/150‑1 (MQ=255)
                                                                                                                                         tataAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTcc          >  1:114102/1‑150 (MQ=255)
                                                                                                                                         tataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTcc          >  1:1005197/1‑150 (MQ=255)
                                                                                                                                          ataAGGATTCGGTCCCCCTCTCTGGAACGTTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGTCTCATTTTTTCGCTCTCACCGGCATCCCTTTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCa         <  1:3646883/150‑1 (MQ=255)
                                                                                                                                           taAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCAc        >  2:3799253/1‑150 (MQ=255)
                                                                                                                                            aaGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACg       >  2:207517/1‑150 (MQ=255)
                                                                                                                                            aaGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACg       <  2:1914038/150‑1 (MQ=255)
                                                                                                                                             aGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGa      >  1:255592/1‑150 (MQ=255)
                                                                                                                                               gATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAAt    >  1:1832747/1‑150 (MQ=255)
                                                                                                                                                aTTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATc   >  2:1018826/1‑150 (MQ=255)
                                                                                                                                                aTTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATc   <  1:2575302/150‑1 (MQ=255)
                                                                                                                                                aTTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATc   >  1:2942786/1‑150 (MQ=255)
                                                                                                                                                 ttCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATCa  <  1:1999289/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATCA  >  CP009273/321076‑321370

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