Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,174,132 G→T 6.8% E187* (GAA→TAA)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,174,1320GT 6.8% 1469.6 / 50.5 379E187* (GAA→TAA) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (163/190);  new base T (15/11);  total (178/201)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 391 total aligned reads displayed.

TCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGC  >  CP009273/2173983‑2174281
                                                                                                                                                     |                                                                                                                                                     
 ctcGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGgaga                                                                                                                                                      >  1:4594393/1‑150 (MQ=255)
 ctcGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGgaga                                                                                                                                                      >  1:3375973/1‑150 (MQ=255)
 ctcGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGgaga                                                                                                                                                      <  2:2554592/150‑1 (MQ=255)
  tcGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGaa                                                                                                                                                     >  2:867863/1‑150 (MQ=255)
  tcGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGaa                                                                                                                                                     <  1:2299172/150‑1 (MQ=255)
  tcGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGaa                                                                                                                                                     >  1:248393/1‑150 (MQ=255)
  tcGCCATGAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGaa                                                                                                                                                     <  1:4305031/150‑1 (MQ=255)
   cGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGaaa                                                                                                                                                    >  2:2677368/1‑150 (MQ=255)
    gCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAg                                                                                                                                                   <  1:4308632/150‑1 (MQ=255)
    gCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAg                                                                                                                                                   >  1:2744204/1‑150 (MQ=255)
    gCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTGATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAg                                                                                                                                                   <  1:3101872/150‑1 (MQ=255)
     ccATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGc                                                                                                                                                  <  2:308271/150‑1 (MQ=255)
      cATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCt                                                                                                                                                 <  1:3911568/150‑1 (MQ=255)
      cATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCt                                                                                                                                                 >  1:4217893/1‑150 (MQ=255)
      cATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCt                                                                                                                                                 <  2:5097175/150‑1 (MQ=255)
       aTTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTg                                                                                                                                                >  2:3260120/1‑150 (MQ=255)
       aTTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGGTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTg                                                                                                                                                <  1:2818045/150‑1 (MQ=255)
         tAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGt                                                                                                                                              <  1:2693182/150‑1 (MQ=255)
          aGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTc                                                                                                                                             >  2:3412085/1‑150 (MQ=255)
           ggTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCa                                                                                                                                            >  2:3809281/1‑150 (MQ=255)
           ggTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCa                                                                                                                                            >  2:421985/1‑150 (MQ=255)
            gttgttTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAg                                                                                                                                           <  2:4294636/150‑1 (MQ=255)
                 ttGCTTCCGGCGCGTGCTGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGGTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGAtt                                                                                                                                      <  2:4972779/150‑1 (MQ=255)
                 ttGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGAtt                                                                                                                                      >  2:5017842/1‑150 (MQ=255)
                 ttGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGAtt                                                                                                                                      >  2:2474477/1‑150 (MQ=255)
                     ttCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGAttgttg                                                                                                                                  >  2:373032/1‑150 (MQ=255)
                       ccGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGAt                                                                                                                                >  2:1204419/1‑150 (MQ=255)
                         ggCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGATAAAGCTGGTCAGCGATTGTTGATTc                                                                                                                              >  2:4031882/1‑150 (MQ=255)
                         ggCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTc                                                                                                                              <  1:1262257/150‑1 (MQ=255)
                            gcgTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCt                                                                                                                           >  2:3982331/1‑150 (MQ=255)
                              gTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGc                                                                                                                         >  1:2950233/1‑150 (MQ=255)
                                gctGTTGATTCCTTTATTGATGCTGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCt                                                                                                                       <  2:2860246/148‑1 (MQ=255)
                                gATGTTGATTCCTTTATTGATGATGTAGCGCAGGCGTCCAGCCCGACCCAAAAATCCGATCTTGCGGTAGCAGGCGCTGAAGACATCGCAATGGCTAATTCACGTGCAATTGAAGCAGAAAGCAGGTCAGCCATAGATGATTcccaccct                                                                                                                       >  1:117135/1‑145 (MQ=255)
                                 aTGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCtt                                                                                                                      <  2:3319657/150‑1 (MQ=255)
                                    ttGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCaa                                                                                                                   <  1:530976/150‑1 (MQ=255)
                                      gATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAAtt                                                                                                                 <  1:2169568/150‑1 (MQ=255)
                                       aTTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTg                                                                                                                <  1:3130245/150‑1 (MQ=255)
                                             ttATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGAc                                                                                                          <  2:807851/150‑1 (MQ=255)
                                              tATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACa                                                                                                         <  2:4480852/150‑1 (MQ=255)
                                                  gatgatGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGca                                                                                                     >  2:2357864/1‑150 (MQ=255)
                                                   atgatgTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCTCGACATGcac                                                                                                    >  2:2761084/1‑150 (MQ=255)
                                                    tgatgTAGCGCTGGCGTCCAGCCCTACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAATCTGGTCAGCGATTGTTGATTCCATGCCTTCAATTGCGCGACATGcaca                                                                                                   <  1:851577/150‑1 (MQ=255)
                                                    tgatgTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGATAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGcaca                                                                                                   >  1:2883367/1‑150 (MQ=255)
                                                      atgTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGATAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGc                                                                                                 >  1:3614982/1‑150 (MQ=255)
                                                      atgTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGc                                                                                                 >  1:1544009/1‑150 (MQ=255)
                                                       tgTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGATAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCa                                                                                                <  2:4360578/150‑1 (MQ=255)
                                                        gTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCaa                                                                                               <  2:1764583/150‑1 (MQ=255)
                                                          aGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCaaaa                                                                                             <  1:4027372/150‑1 (MQ=255)
                                                           gcgcTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCaaaaa                                                                                            <  1:1835454/150‑1 (MQ=255)
                                                            cgcTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAAc                                                                                           <  2:1333623/150‑1 (MQ=255)
                                                             gcTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACg                                                                                          >  2:914050/1‑150 (MQ=255)
                                                               tGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCa                                                                                        >  1:3587709/1‑150 (MQ=255)
                                                                   gTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCAc                                                                                    >  2:2268332/1‑150 (MQ=255)
                                                                    tCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACt                                                                                   >  2:4147505/1‑150 (MQ=255)
                                                                      cAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTAc                                                                                 >  2:158954/1‑150 (MQ=255)
                                                                      cAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTAc                                                                                 >  2:1875015/1‑150 (MQ=255)
                                                                          ccGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTc                                                                             >  2:2038265/1‑150 (MQ=255)
                                                                            gACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAg                                                                           >  1:5031697/1‑150 (MQ=255)
                                                                            gACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAg                                                                           >  1:3710380/1‑150 (MQ=255)
                                                                             aCCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAgc                                                                          <  2:1040073/150‑1 (MQ=255)
                                                                             aCCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAgc                                                                          >  1:1880163/1‑150 (MQ=255)
                                                                             aCCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAgc                                                                          <  1:948157/150‑1 (MQ=255)
                                                                               ccATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAgcgc                                                                        >  2:1709671/1‑150 (MQ=255)
                                                                                cATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGcc                                                                       >  1:656688/1‑150 (MQ=255)
                                                                                     aTCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACTGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACt                                                                  <  2:1791720/150‑1 (MQ=255)
                                                                                      tCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTg                                                                 >  2:1540643/1‑150 (MQ=255)
                                                                                        cGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTggc                                                               >  1:1631911/1‑150 (MQ=255)
                                                                                         gATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTggcg                                                              >  1:1059090/1‑150 (MQ=255)
                                                                                         gATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTggcg                                                              <  1:2275266/150‑1 (MQ=255)
                                                                                              ggCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGgcggcgcg                                                         <  2:4615874/150‑1 (MQ=255)
                                                                                                cGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGgcggcgcggc                                                       <  1:3652017/150‑1 (MQ=255)
                                                                                                 ggTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGATAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCt                                                      <  1:2455940/150‑1 (MQ=255)
                                                                                                 ggTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCt                                                      <  1:4415903/150‑1 (MQ=255)
                                                                                                 ggTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCt                                                      <  1:4903307/150‑1 (MQ=255)
                                                                                                 ggTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCt                                                      <  1:4021639/150‑1 (MQ=255)
                                                                                                    tGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGga                                                   <  2:3541021/150‑1 (MQ=255)
                                                                                                    tGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGga                                                   <  1:1699825/150‑1 (MQ=255)
                                                                                                    tGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGga                                                   >  2:4173974/1‑150 (MQ=255)
                                                                                                        ggCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGAtt                                               <  1:808742/150‑1 (MQ=255)
                                                                                                         gcgcGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTg                                              <  1:2654262/150‑1 (MQ=255)
                                                                                                          cgcgGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTgc                                             >  2:4721393/1‑150 (MQ=255)
                                                                                                             ggTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGATAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCt                                          >  2:1599894/1‑150 (MQ=255)
                                                                                                             ggTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCt                                          >  1:4117934/1‑150 (MQ=255)
                                                                                                                aGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGaa                                       >  2:389902/1‑150 (MQ=255)
                                                                                                                     tCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATtg                                  >  1:3218213/1‑150 (MQ=255)
                                                                                                                       gCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATtgtg                                <  2:1610431/150‑1 (MQ=255)
                                                                                                                       gCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATtgtg                                >  2:3608367/1‑150 (MQ=255)
                                                                                                                           gggCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCa                            <  1:2038265/150‑1 (MQ=255)
                                                                                                                             gCGATTTCTCGTGCCATTGACGGAGAAAGCTGGACAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAAt                          <  2:1845406/150‑1 (MQ=255)
                                                                                                                              cGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATg                         <  1:3095744/150‑1 (MQ=255)
                                                                                                                              cGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATg                         <  1:3000338/150‑1 (MQ=255)
                                                                                                                               gATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGc                        >  2:617877/1‑150 (MQ=255)
                                                                                                                                   tctcGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGAc                    >  2:484266/1‑150 (MQ=255)
                                                                                                                                     tcGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACgg                  >  2:2010355/1‑150 (MQ=255)
                                                                                                                                      cGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGc                 <  2:3468237/150‑1 (MQ=255)
                                                                                                                                              tGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATc         <  1:3533352/150‑1 (MQ=255)
                                                                                                                                                 cGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGc      >  1:3822166/1‑150 (MQ=255)
                                                                                                                                                  ggAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGcc     >  2:5140605/1‑150 (MQ=255)
                                                                                                                                                  ggAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGcc     >  2:2598532/1‑150 (MQ=255)
                                                                                                                                                     gaAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGc  >  1:1416239/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGC  >  CP009273/2173983‑2174281

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: