Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,631,894 C→T 11.8% S49N (AGC→AAC)  gnsB ← Qin prophage; multicopy suppressor of secG(Cs) and fabA6(Ts)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,631,8940CT11.8% 1266.7 / 113.6 416S49N (AGC→AAC) gnsBQin prophage; multicopy suppressor of secG(Cs) and fabA6(Ts)
Reads supporting (aligned to +/- strand):  ref base C (189/177);  new base T (31/18);  total (221/195)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 418 total aligned reads displayed.

GCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAGTTTTCAATATTCATCAT  >  CP009273/1631745‑1632042
                                                                                                                                                     |                                                                                                                                                    
 ccGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCt                                                                                                                                                     <  2:4062306/150‑1 (MQ=255)
 ccGGACTGCATTGGTTTAATAAAAACCATCAACAAGTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCt                                                                                                                                                     <  2:4431047/150‑1 (MQ=255)
    gACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTc                                                                                                                                                  >  1:1000974/1‑150 (MQ=255)
    gACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATACGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTc                                                                                                                                                  >  1:2872050/1‑150 (MQ=255)
     aCTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCa                                                                                                                                                 >  2:1441864/1‑150 (MQ=255)
     aCTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCa                                                                                                                                                 >  1:3132248/1‑150 (MQ=255)
      cTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCaa                                                                                                                                                <  1:5169307/150‑1 (MQ=255)
       tGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAg                                                                                                                                               <  1:3251098/150‑1 (MQ=255)
       tGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAg                                                                                                                                               >  1:1226897/1‑150 (MQ=255)
          aTTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGAcc                                                                                                                                            >  1:2229554/1‑150 (MQ=255)
           ttGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGAccc                                                                                                                                           >  1:998908/1‑150 (MQ=255)
             ggTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGt                                                                                                                                         <  1:3662020/150‑1 (MQ=255)
                ttAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCAt                                                                                                                                      <  1:1277617/150‑1 (MQ=255)
                 taataaAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCAtt                                                                                                                                     >  1:3702491/1‑150 (MQ=255)
                 taataaAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCAtt                                                                                                                                     >  2:916245/1‑150 (MQ=255)
                   ataaAAACAATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTATATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCAtaat                                                                                                                                   >  2:1506181/1‑147 (MQ=255)
                      aaaaCCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAAAGAGGTTTATGGATAAAATTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACACGTCATTTTTTc                                                                                                                                >  2:2787581/1‑150 (MQ=255)
                       aaaCCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTcc                                                                                                                               <  2:3518259/150‑1 (MQ=255)
                        aaCCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTccc                                                                                                                              >  2:1309225/1‑150 (MQ=255)
                           cATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACCCGTCATTTTTTCCCGTg                                                                                                                           >  2:1680342/1‑150 (MQ=255)
                              caacaaTTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCgg                                                                                                                        >  1:3193519/1‑150 (MQ=255)
                              caacaaTTGTGATTTAGATATTCGGAACCATTCAAAAATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATAAATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACCCGTCATTTTTTCACGTGCgg                                                                                                                        >  1:2619865/1‑150 (MQ=255)
                               aacaaTTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGt                                                                                                                       >  1:2113196/1‑150 (MQ=255)
                                 cacttGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTaa                                                                                                                     <  1:681903/147‑1 (MQ=255)
                                        gATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAAt                                                                                                              <  2:3041575/150‑1 (MQ=255)
                                         aTTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATa                                                                                                             <  2:3284281/150‑1 (MQ=255)
                                           ttAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAATACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACt                                                                                                           <  2:3136740/150‑1 (MQ=255)
                                               atatTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTa                                                                                                       >  2:2867140/1‑150 (MQ=255)
                                                   tCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAActt                                                                                                   <  1:1685768/150‑1 (MQ=255)
                                                    cGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTAGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAActtc                                                                                                  >  1:1392206/1‑150 (MQ=255)
                                                     ggAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAActtct                                                                                                 <  2:4897268/150‑1 (MQ=255)
                                                        aCCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTc                                                                                              >  2:1604240/1‑150 (MQ=255)
                                                         ccATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTcc                                                                                             >  2:3729994/1‑150 (MQ=255)
                                                         ccATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCAGTCATTTTTTCCCGTGCGGTAAACTGAATAATGGTAACTTCTGTcc                                                                                             >  2:4508380/1‑150 (MQ=255)
                                                          cATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAAATTCTTTccc                                                                                            >  1:924185/1‑150 (MQ=255)
                                                          cATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTccc                                                                                            <  1:809599/150‑1 (MQ=255)
                                                          cATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTccc                                                                                            >  2:1265051/1‑150 (MQ=255)
                                                              cAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTc                                                                                        >  1:3805955/1‑150 (MQ=255)
                                                                aaTATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTCGCTTCTGTCCCGGTCtt                                                                                      >  1:4223730/1‑150 (MQ=255)
                                                                  tataACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCtttt                                                                                    <  1:3095102/150‑1 (MQ=255)
                                                                  tataACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCtttt                                                                                    <  1:2369841/150‑1 (MQ=255)
                                                                    taACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTc                                                                                  <  2:1945445/150‑1 (MQ=255)
                                                                     aaCAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCt                                                                                 >  1:4816133/1‑150 (MQ=255)
                                                                     aaCAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCt                                                                                 <  2:2383053/150‑1 (MQ=255)
                                                                      aCAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGTTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCtt                                                                                <  2:5259733/150‑1 (MQ=255)
                                                                      aCAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCtt                                                                                <  2:677677/150‑1 (MQ=255)
                                                                       cAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGATCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTTTGTTCTTTCCCGGTCTTTTTCTTa                                                                               >  1:867641/1‑150 (MQ=255)
                                                                       cAAAACCCCGTAAAAACGAGGTTTATGGATAAATATTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTa                                                                               <  2:2468985/150‑1 (MQ=255)
                                                                          aaCCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGt                                                                            <  2:3930710/150‑1 (MQ=255)
                                                                            ccccGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTc                                                                          <  1:1584204/150‑1 (MQ=255)
                                                                            ccccCTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTc                                                                          <  2:2904856/150‑1 (MQ=255)
                                                                                  aaaaaCGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCaataat                                                                    >  1:4586701/1‑150 (MQ=255)
                                                                                  aaaaaCGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCaataat                                                                    <  1:2802499/150‑1 (MQ=255)
                                                                                  aaaaaCGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCaataat                                                                    >  2:3297904/1‑150 (MQ=255)
                                                                                    aaaCGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTGACTTCTTTCCCGGTCTTTTTATTAAGTTAAATAAttt                                                                  >  1:1855369/1‑150 (MQ=255)
                                                                                    aaaCGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAAttt                                                                  >  2:4779919/1‑150 (MQ=255)
                                                                                      aCGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAAttttt                                                                >  1:1945696/1‑150 (MQ=255)
                                                                                        gAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAAttttttt                                                              <  1:1777938/150‑1 (MQ=255)
                                                                                         aGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAAtttttttt                                                             >  2:5012471/1‑150 (MQ=255)
                                                                                          ggTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTg                                                            <  1:2468172/150‑1 (MQ=255)
                                                                                           gTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGt                                                           >  1:153261/1‑150 (MQ=255)
                                                                                            tttATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGtt                                                          <  2:2294481/150‑1 (MQ=255)
                                                                                             ttATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATGGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTta                                                         <  2:3420313/150‑1 (MQ=255)
                                                                                              tATGGATAAATTTTATTATTGAATACATCAGCTTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTtat                                                        <  2:4291970/150‑1 (MQ=255)
                                                                                              tATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTtat                                                        >  1:3355864/1‑150 (MQ=255)
                                                                                              tATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTtat                                                        <  1:4199458/150‑1 (MQ=255)
                                                                                                tGGATAAATTTTATTATTGAAGACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTATTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTtatat                                                      <  2:4107717/150‑1 (MQ=255)
                                                                                                   aTAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAAGTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATAtt                                                   >  1:477049/1‑150 (MQ=255)
                                                                                                        ttttATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGaa                                              <  1:4416843/150‑1 (MQ=255)
                                                                                                          ttattaTTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAat                                            <  2:4566781/150‑1 (MQ=255)
                                                                                                          ttattaTTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAat                                            <  2:2005127/150‑1 (MQ=255)
                                                                                                           tattatTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAata                                           <  2:2681388/150‑1 (MQ=255)
                                                                                                            attattGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAGCTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAatat                                          >  1:4252391/1‑150 (MQ=255)
                                                                                                            attattGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTACCTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAatat                                          >  1:1174447/1‑150 (MQ=255)
                                                                                                            attattGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAatat                                          <  2:1159415/150‑1 (MQ=255)
                                                                                                            attattGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTAACGTGCGGTAAACTGAATACTGGTAACTTCTTTACCGGTCTTTTTATTAAGTTCAATAATTTTTTATGTTATATATTCAGAAataa                                          >  2:3591917/1‑149 (MQ=255)
                                                                                                             ttattGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtc                                         <  2:4359029/150‑1 (MQ=255)
                                                                                                               attGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctg                                       >  1:651642/1‑150 (MQ=255)
                                                                                                               attGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctg                                       >  2:4190279/1‑150 (MQ=255)
                                                                                                               attGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctg                                       >  2:86659/1‑150 (MQ=255)
                                                                                                                 tGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgct                                     <  1:4911720/150‑1 (MQ=255)
                                                                                                                 tGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgct                                     <  1:5010479/150‑1 (MQ=255)
                                                                                                                 tGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgct                                     <  1:5283624/150‑1 (MQ=255)
                                                                                                                   aaTACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgcttc                                   >  2:32212/1‑150 (MQ=255)
                                                                                                                    aTACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgcttct                                  <  2:3648187/150‑1 (MQ=255)
                                                                                                                    aTACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgcttct                                  <  2:2113196/150‑1 (MQ=255)
                                                                                                                     tACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgcttctg                                 >  2:3875188/1‑150 (MQ=255)
                                                                                                                     tACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgcttctg                                 >  1:5198753/1‑150 (MQ=255)
                                                                                                                      aCATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATAtctgcttctgc                                >  2:809160/1‑150 (MQ=255)
                                                                                                                          cAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCtttt                            >  2:4542379/1‑150 (MQ=255)
                                                                                                                          cAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCtttt                            <  1:3888234/150‑1 (MQ=255)
                                                                                                                               taaattaaTCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGtttt                       <  2:310925/150‑1 (MQ=255)
                                                                                                                                 aattaaTCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTa                     >  1:3601112/1‑150 (MQ=255)
                                                                                                                                 aattaaTCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTa                     <  2:845574/150‑1 (MQ=255)
                                                                                                                                  attaaTCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTaa                    >  1:4693367/1‑150 (MQ=255)
                                                                                                                                   ttaaTCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAg                   <  1:394663/150‑1 (MQ=255)
                                                                                                                                     aaTCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAGtt                 <  2:5148378/150‑1 (MQ=255)
                                                                                                                                         ttGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAGTTTTCa             <  1:1074184/150‑1 (MQ=255)
                                                                                                                                          tGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAGTTTTCaa            <  2:591076/150‑1 (MQ=255)
                                                                                                                                                  tAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTATTAATTTAAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAGTTTTCAATATtcata    >  2:468028/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCGGACTGCATTGGTTTAATAAAAACCATCAACAATTGTGATTTAGATATTCGGAACCATTCAAATATAACAAAACCCCGTAAAAACGAGGTTTATGGATAAATTTTATTATTGAATACATCAGATTAAATTAATCTTGACATCATAGCTTTCAAGACCCGTCATTTTTTCCCGTGCGGTAAACTGAATACTGGTAACTTCTTTCCCGGTCTTTTTCTTAAGTTCAATAATTTTTTTTGTTATATATTCAGAAATATCTGCTTCTGCTTTTGTTTTTAAGTTTTCAATATTCATCAT  >  CP009273/1631745‑1632042

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