Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,833,730 C→A 100% R3S (CGT→AGT)  ynjE → molybdopterin synthase sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,833,7300CA96.7% 994.0 / 20.6 307R3S (CGT→AGT) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base A (173/124);  total (177/130)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 311 total aligned reads displayed.

CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833583‑1833876
                                                                                                                                                   |                                                                                                                                                  
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                  >  1:2700453/1‑150 (MQ=255)
 tcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAACCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtg                                                                                                                                                 <  2:340988/150‑1 (MQ=255)
 tcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtg                                                                                                                                                 <  2:1214249/150‑1 (MQ=255)
 tcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtg                                                                                                                                                 <  2:3714812/150‑1 (MQ=255)
   cAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGtt                                                                                                                                               <  2:2143961/150‑1 (MQ=255)
    aGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGttt                                                                                                                                              >  2:1822204/1‑150 (MQ=255)
        tGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCa                                                                                                                                          >  2:4481780/1‑150 (MQ=255)
         gTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACCAACTGCGATATCGGAAAGTTTTTCTCAATGCGTCAGTTCAGATTGGCGCACTGTATAGTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaa                                                                                                                                         <  1:286742/150‑1 (MQ=255)
         gTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaa                                                                                                                                         >  2:2562732/1‑150 (MQ=255)
          ttCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaaa                                                                                                                                        <  1:1822204/150‑1 (MQ=255)
            ccTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATg                                                                                                                                      >  1:4375978/1‑150 (MQ=255)
            ccTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATg                                                                                                                                      >  2:1993967/1‑150 (MQ=255)
              tgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGAc                                                                                                                                    <  1:1070568/150‑1 (MQ=255)
              tgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGAc                                                                                                                                    >  1:3275093/1‑150 (MQ=255)
              tgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGAc                                                                                                                                    >  2:3954972/1‑150 (MQ=255)
              tgatgATAGTTCTGTTCTGGAAATGGCGCAGTGGTCAGAAAATTACAACAAACTGCTATAACGCAAAGTTTTTCACAATGCGTCACTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGAc                                                                                                                                    >  2:271056/4‑150 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTTAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                  >  2:4008366/1‑150 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                  >  2:2677278/1‑150 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                  >  2:46184/1‑150 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                  <  2:705727/150‑1 (MQ=255)
                    tAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACCACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCt                                                                                                                              <  2:3210338/150‑1 (MQ=255)
                    tAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCt                                                                                                                              >  1:3977215/1‑150 (MQ=255)
                     aGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTg                                                                                                                             <  1:559203/150‑1 (MQ=255)
                      gttctgttctGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTgg                                                                                                                            <  1:3221813/150‑1 (MQ=255)
                       ttctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGc                                                                                                                           >  1:3571956/1‑150 (MQ=255)
                       ttctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGACGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGc                                                                                                                           <  2:3764135/150‑1 (MQ=255)
                        tctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCa                                                                                                                          <  1:483344/150‑1 (MQ=255)
                         ctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCaa                                                                                                                         >  2:743178/1‑150 (MQ=255)
                            ttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATgg                                                                                                                      >  1:515981/1‑150 (MQ=255)
                             tctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGc                                                                                                                     >  2:2291202/1‑150 (MQ=255)
                               tgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCtt                                                                                                                   >  2:533224/1‑150 (MQ=255)
                                 gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTa                                                                                                                 <  1:3389132/150‑1 (MQ=255)
                                         gcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggc                                                                                                         >  2:2443564/1‑150 (MQ=255)
                                                tCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGct                                                                                                  >  1:3850363/1‑150 (MQ=255)
                                                 cAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGctt                                                                                                 <  1:748100/150‑1 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              >  2:4338954/1‑150 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              <  2:3734558/150‑1 (MQ=255)
                                                       aTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGt                                                                                           <  1:88861/150‑1 (MQ=255)
                                                       aTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGt                                                                                           >  2:2638775/1‑150 (MQ=255)
                                                        ttACATCAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTg                                                                                          >  1:4481461/1‑150 (MQ=255)
                                                         tACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                         >  2:4511663/1‑150 (MQ=255)
                                                            aacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGcc                                                                                      >  1:3227603/1‑150 (MQ=255)
                                                             acaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCg                                                                                     >  1:2304577/1‑150 (MQ=255)
                                                               aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTATGGCTGGCTTGCGCTTCTTCGTGGGCCGCt                                                                                   >  2:3126658/1‑150 (MQ=255)
                                                               aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCt                                                                                   >  1:1484116/1‑150 (MQ=255)
                                                               aaaCAGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCt                                                                                   >  1:3062554/1‑150 (MQ=255)
                                                                aaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTg                                                                                  >  2:4246626/1‑150 (MQ=255)
                                                                aaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTg                                                                                  >  2:3515524/1‑150 (MQ=255)
                                                                aaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTg                                                                                  >  2:162972/1‑150 (MQ=255)
                                                                 aCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGa                                                                                 >  2:1971976/1‑150 (MQ=255)
                                                                  cTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGaa                                                                                >  1:358989/1‑150 (MQ=255)
                                                                      gATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAAAGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTgg                                                                            >  2:1031455/1‑150 (MQ=255)
                                                                         aaCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGa                                                                         >  1:2185462/1‑150 (MQ=255)
                                                                            gCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGc                                                                      >  2:1226538/1‑150 (MQ=255)
                                                                            gCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGc                                                                      <  1:1914241/150‑1 (MQ=255)
                                                                              aaaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCct                                                                    <  2:2830004/150‑1 (MQ=255)
                                                                               aaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCctc                                                                   >  1:472931/1‑150 (MQ=255)
                                                                                aGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCctct                                                                  >  1:684754/1‑150 (MQ=255)
                                                                                aGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCctct                                                                  >  2:2930283/1‑150 (MQ=255)
                                                                                 gTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCtt                                                                 >  1:389322/1‑150 (MQ=255)
                                                                                  tttttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTa                                                                >  2:736879/1‑150 (MQ=255)
                                                                                       ctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                           <  2:4103677/150‑1 (MQ=255)
                                                                                         cAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGa                                                         >  1:1005925/1‑150 (MQ=255)
                                                                                         cAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGa                                                         >  2:1670215/1‑150 (MQ=255)
                                                                                            tGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCa                                                      <  2:1294597/150‑1 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                   >  2:3196895/1‑150 (MQ=255)
                                                                                                 cAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTc                                                 >  1:2771507/1‑150 (MQ=255)
                                                                                                  aGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTca                                                >  1:840027/1‑150 (MQ=255)
                                                                                                   gTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaa                                               >  2:1131739/1‑150 (MQ=255)
                                                                                                    ttCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaac                                              <  1:1586700/150‑1 (MQ=255)
                                                                                                      cAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaa                                            >  1:4423728/1‑150 (MQ=255)
                                                                                                       aGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaac                                           >  1:4308062/1‑150 (MQ=255)
                                                                                                         aaTGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaacaa                                         >  1:1467112/1‑150 (MQ=255)
                                                                                                          aTGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaa                                        <  1:4089390/150‑1 (MQ=255)
                                                                                                           tGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaaa                                       >  2:595976/1‑150 (MQ=255)
                                                                                                           tGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaaa                                       >  2:4187567/1‑150 (MQ=255)
                                                                                                               gcACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGc                                   <  2:2796775/150‑1 (MQ=255)
                                                                                                               gcACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGc                                   >  1:3340471/1‑150 (MQ=255)
                                                                                                               gcACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGc                                   <  1:504012/150‑1 (MQ=255)
                                                                                                                  cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaaa                                <  1:263103/150‑1 (MQ=255)
                                                                                                                  cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaaa                                <  2:1352011/150‑1 (MQ=255)
                                                                                                                   tCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAg                               <  1:4432041/150‑1 (MQ=255)
                                                                                                                    cAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGc                              <  1:2152347/150‑1 (MQ=255)
                                                                                                                     aaaaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCg                             >  1:2411926/1‑150 (MQ=255)
                                                                                                                     aaaaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCg                             <  2:2522740/150‑1 (MQ=255)
                                                                                                                         cTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatag                         <  1:4213096/150‑1 (MQ=255)
                                                                                                                          tACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgataga                        <  2:2058756/150‑1 (MQ=255)
                                                                                                                           aCAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagat                       >  1:2537986/1‑150 (MQ=255)
                                                                                                                            cAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagata                      >  1:4141283/1‑150 (MQ=255)
                                                                                                                                gTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCg                  >  1:4052534/1‑150 (MQ=255)
                                                                                                                                gTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCg                  >  2:1729001/1‑150 (MQ=255)
                                                                                                                                gTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCg                  <  2:2228524/150‑1 (MQ=255)
                                                                                                                                 tCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGc                 <  1:1768060/150‑1 (MQ=255)
                                                                                                                                        tttCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAgcg          >  1:1887058/1‑150 (MQ=255)
                                                                                                                                        tttCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAgcg          >  2:3338661/1‑150 (MQ=255)
                                                                                                                                           cGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGt       >  2:4576814/1‑150 (MQ=255)
                                                                                                                                             aTGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGttt     >  1:1412213/1‑150 (MQ=255)
                                                                                                                                             aTGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGttt     <  2:2329135/150‑1 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    >  1:208197/1‑150 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    >  2:151173/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833583‑1833876

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: ‑