Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,076,055 T→G 12.4% Q22P (CAG→CCG)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,076,0550TG12.4% 709.4 / 8.5 326Q22P (CAG→CCG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (77/176);  new base G (43/3);  total (147/179)
Fisher's exact test for biased strand distribution p-value = 2.06e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 329 total aligned reads displayed.

TCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGC  >  CP009273/4075908‑4076204
                                                                                                                                                   |                                                                                                                                                      
  aTTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAgg                                                                                                                                                    >  1:4201476/1‑150 (MQ=255)
    tGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcg                                                                                                                                                  >  2:3252762/1‑150 (MQ=255)
     ggggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                                                                                                                 >  1:634413/1‑150 (MQ=255)
       ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                                                                                                               >  1:1090043/1‑150 (MQ=255)
        gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCGGGGGCGCggg                                                                                                                                              >  1:2871905/1‑150 (MQ=255)
          tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTGGCCCGCCGGAGGCGCcgggg                                                                                                                                            >  1:4421451/1‑150 (MQ=25)
           gTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGt                                                                                                                                           <  1:1877645/150‑1 (MQ=255)
             cGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCGGATGCGCGGGGGTCa                                                                                                                                         >  1:4488343/1‑150 (MQ=255)
              gTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAg                                                                                                                                        >  2:2692166/1‑150 (MQ=255)
               tACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGa                                                                                                                                       >  1:3065865/1‑150 (MQ=255)
                acacCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCCGAc                                                                                                                                      >  2:3298514/1‑150 (MQ=255)
                acacCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGAc                                                                                                                                      >  2:3853039/1‑150 (MQ=255)
                acacCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCCGAGGCCCGGGGGGCAgcc                                                                                                                                      >  2:3616995/1‑148 (MQ=255)
                  acCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                                                                                                    >  2:4526838/1‑150 (MQ=255)
                    cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTGCCCCGCCCGGGGCGCGGGGGGCAGACCggt                                                                                                                                  >  1:3307263/1‑148 (MQ=25)
                     ccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                                                                                                                 <  2:724737/150‑1 (MQ=255)
                      cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCACAcccgtag                                                                                                                                >  1:3114101/1‑145 (MQ=255)
                       aaTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTGCCGCGGCGGAGGCGCGGGGGGGAGAcccgtttt                                                                                                                               >  1:82242/1‑144 (MQ=25)
                       aaTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                                                                                                               <  2:503952/150‑1 (MQ=255)
                       aaTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCAGACCGTTAGt                                                                                                                               >  2:4321560/1‑150 (MQ=255)
                        aTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                                                                                                              <  1:1714991/150‑1 (MQ=255)
                        aTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTGAGTc                                                                                                                              >  2:2373533/1‑150 (MQ=255)
                          gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                                                                                                            >  2:3491149/1‑150 (MQ=255)
                           ttATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgcg                                                                                                                           <  1:1032183/150‑1 (MQ=255)
                             aTTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGTGGCGCGGGGGTCAGACCGTTAGTCGCGGa                                                                                                                         >  2:4707031/1‑150 (MQ=255)
                              ttGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACCGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGaa                                                                                                                        <  1:3346686/150‑1 (MQ=255)
                                  cGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCCTTAGTCGCGGgacgac                                                                                                                    >  1:2391162/1‑150 (MQ=255)
                                   gACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgacg                                                                                                                   <  2:1970752/150‑1 (MQ=255)
                                           cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                                                                                                           <  2:301587/150‑1 (MQ=255)
                                             cgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGGCAGCCCGTTAGTCGCGGAACGACGACTGatatag                                                                                                         >  2:4125514/1‑148 (MQ=255)
                                              gAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                        <  1:2627382/150‑1 (MQ=255)
                                               aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTGAGACCGGTAGTCGCGGAACGACGAATGATATATTg                                                                                                       >  2:530103/1‑150 (MQ=255)
                                                 tGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCTATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgag                                                                                                     <  1:1188384/150‑1 (MQ=255)
                                                 tGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCAGACCGGTAGTCGCGGAACGACGAATGATATATTgtg                                                                                                     >  1:1542313/1‑148 (MQ=255)
                                                  gTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                                    <  2:2078911/150‑1 (MQ=255)
                                                  gTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                                    >  2:3948331/1‑150 (MQ=255)
                                                   tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAt                                                                                                   <  1:1951553/150‑1 (MQ=255)
                                                   tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAt                                                                                                   >  2:4346210/1‑150 (MQ=255)
                                                   tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGGTAGTCGCGGAACGACGAAAGATATCTTGAGAt                                                                                                   >  2:3902012/1‑150 (MQ=255)
                                                    cGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                                  <  2:2022468/150‑1 (MQ=255)
                                                    cGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                                  <  2:1078112/150‑1 (MQ=255)
                                                    cGTGCCACTCTGAACACGCCCCCTGACAGGCTTTACAGCCGATACGGGTGGTAACGTCGAGGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                                  <  1:3222291/150‑1 (MQ=38)
                                                        ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGTGATTGATa                                                                                             >  2:2516844/1‑150 (MQ=255)
                                                        ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                             <  2:2883528/150‑1 (MQ=255)
                                                        ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                             <  1:4166750/150‑1 (MQ=255)
                                                         cACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTTAGATTGATaa                                                                                            >  1:4728128/1‑150 (MQ=255)
                                                          aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg                                                                                           <  1:3976141/150‑1 (MQ=255)
                                                           ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc                                                                                          <  2:1514297/150‑1 (MQ=255)
                                                            tctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGcc                                                                                         <  1:98776/150‑1 (MQ=255)
                                                             ctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCa                                                                                        <  2:658043/150‑1 (MQ=255)
                                                             ctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACCTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCa                                                                                        <  1:1404120/150‑1 (MQ=255)
                                                              tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCat                                                                                       <  2:2256318/150‑1 (MQ=255)
                                                              tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCGGTGGCGCGGGGGGCAGACCGGTAGTCGCGGGACCCCGAAAGATATATTGTGATTGATAAGCCaa                                                                                       >  2:817667/1‑149 (MQ=25)
                                                               gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGTTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata                                                                                      <  2:1163503/150‑1 (MQ=255)
                                                               gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata                                                                                      <  1:190125/150‑1 (MQ=255)
                                                                aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGTGGCGCGGGGGTAAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAACCaaat                                                                                     >  2:4157271/1‑150 (MQ=255)
                                                                aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCTGAGGCGCGGGGGGCAGACCCTTAGTCGCGGGACCACGCATGAAATCTTGAGAGTGATAAACCaaaa                                                                                     >  1:4287074/1‑147 (MQ=25)
                                                                 acacGCCACCTTACAGGATTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg                                                                                    <  1:3646472/150‑1 (MQ=255)
                                                                   acGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtc                                                                                  <  2:2320213/150‑1 (MQ=255)
                                                                    cGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGATGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcg                                                                                 <  2:791511/150‑1 (MQ=255)
                                                                      ccACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgtc                                                                               <  2:1771168/150‑1 (MQ=255)
                                                                        aCCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtctc                                                                             <  1:400455/150‑1 (MQ=255)
                                                                         ccTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                                                                            <  2:3754374/150‑1 (MQ=255)
                                                                           tGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGTGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGTTTGTTAATCCATATTTCGTCTCCGt                                                                          >  2:2568728/1‑150 (MQ=39)
                                                                             aCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTa                                                                        <  1:3619891/150‑1 (MQ=255)
                                                                                 gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCGGAGGCGCGGGGGTCAGAACCGTAGTAGCGGAACCACCAAAGATATATTTTGATTGATAAAACCAAATTCGTCTCCCGTACAcc                                                                    >  1:2154950/1‑150 (MQ=19)
                                                                                   tttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtt                                                                  <  1:4543903/150‑1 (MQ=255)
                                                                                    ttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCttt                                                                 <  1:703465/150‑1 (MQ=255)
                                                                                    ttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAAGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCttt                                                                 >  2:1732232/1‑150 (MQ=255)
                                                                                    ttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCCGTAGTCGCGGAACGACGAATGGTATGTTGGGGTTGATAAGCCAAATGTCGTCTCCGTTACACCgtt                                                                 >  2:3275861/1‑150 (MQ=255)
                                                                                     tACAGCCGAGACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt                                                                <  1:4573631/150‑1 (MQ=255)
                                                                                      aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTc                                                               <  1:2263017/150‑1 (MQ=255)
                                                                                            gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt                                                         <  2:2995504/150‑1 (MQ=255)
                                                                                            gATACAGGTGGTAACGTCGATGAGTTTCGCCACGTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCCGAAAGACGCAAGATATCGTGGGATTGTTAAGCCCAAATTCGTCTCCGTTACACCATTTCcacaat                                                         >  2:1720401/1‑148 (MQ=255)
                                                                                             aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc                                                        <  1:1173553/150‑1 (MQ=255)
                                                                                             aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc                                                        <  1:3555169/150‑1 (MQ=255)
                                                                                             aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc                                                        <  2:4028360/150‑1 (MQ=255)
                                                                                                cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAATCCATATGTCGTCTCCGTTACACCTTTTCCACATTcaca                                                     >  2:1140126/1‑150 (MQ=255)
                                                                                                 aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTTAGATTGATAAGCCATATGTAGTCTACGTTACAACATTTCCACATTCACaa                                                    >  2:317943/1‑150 (MQ=255)
                                                                                                  ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAATCCAAATGTCGTCTCCGTTACACCATTTCCACATTAACaag                                                   >  2:1483672/1‑150 (MQ=255)
                                                                                                  ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAAGATATCTTGAGATTGATAAGCCAAATGTCGTCTCCGTTACACCATTTACACAATAACaag                                                   >  2:125877/1‑150 (MQ=255)
                                                                                                   gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaagg                                                  <  2:4201476/150‑1 (MQ=255)
                                                                                                     ggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaa                                                <  1:2658230/150‑1 (MQ=255)
                                                                                                        aaCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAaggaagga                                             <  2:1545033/150‑1 (MQ=255)
                                                                                                            tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTa                                         <  2:2503733/150‑1 (MQ=255)
                                                                                                             cGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                        <  1:2733364/150‑1 (MQ=255)
                                                                                                                   gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCg                                  <  1:912156/150‑1 (MQ=255)
                                                                                                                            cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGGTAGTCGCGGAACGACGAATGATATCGTGAGATTGATAAGCCAAATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                         >  2:3126652/1‑150 (MQ=255)
                                                                                                                                ttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGtt                     <  2:1603759/150‑1 (MQ=255)
                                                                                                                                      gAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAGTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATc               <  2:3652175/150‑1 (MQ=255)
                                                                                                                                        aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAc             >  2:2714430/1‑150 (MQ=255)
                                                                                                                                           cccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCAGTTACACCTTTTACACATTCACAATGAAGGACTTAAACTCCGGCGTCTGCGTGTTCTCAACACCaa          >  1:54329/1‑150 (MQ=255)
                                                                                                                                            ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAAc         <  1:3446493/150‑1 (MQ=255)
                                                                                                                                            ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAAc         <  1:3853039/150‑1 (MQ=255)
                                                                                                                                             cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACg        >  2:2120431/1‑150 (MQ=255)
                                                                                                                                               cgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGaa      >  2:4407482/1‑150 (MQ=255)
                                                                                                                                               cgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGaa      >  1:3029824/1‑150 (MQ=255)
                                                                                                                                                gcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAAt     <  2:4488166/150‑1 (MQ=255)
                                                                                                                                                 ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg    >  2:4633048/1‑150 (MQ=255)
                                                                                                                                                 ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg    >  1:1567138/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGC  >  CP009273/4075908‑4076204

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