Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,632,204 C→A 31.5% intergenic (‑165/+7) gnsB ← / ← ynfN Qin prophage; multicopy suppressor of secG(Cs) and fabA6(Ts)/Qin prophage; cold shock‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,632,2040CA31.5% 133.4 / 80.6 79intergenic (‑165/+7)gnsB/ynfNQin prophage; multicopy suppressor of secG(Cs) and fabA6(Ts)/Qin prophage; cold shock‑induced protein
Reads supporting (aligned to +/- strand):  ref base C (34/20);  new base A (21/3);  total (55/24)
Fisher's exact test for biased strand distribution p-value = 3.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

TTATGACTTTCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCGCCATT  >  CP009273/1632059‑1632351
                                                                                                                                                 |                                                                                                                                                   
ttATGACTTTCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCatca                                                                                                                                                 >  2:1166489/1‑150 (MQ=255)
ttATGACTTTCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCatca                                                                                                                                                 >  2:1166827/1‑150 (MQ=255)
    gACTTTCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCatcaaatc                                                                                                                                             <  2:4322945/150‑1 (MQ=255)
    gACTTTCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCatcaaatc                                                                                                                                             <  2:4321492/150‑1 (MQ=255)
       tttCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAAc                                                                                                                                          >  1:1240199/1‑150 (MQ=255)
         tCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTc                                                                                                                                        <  1:5157461/150‑1 (MQ=255)
               tACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGc                                                                                                                                  >  1:3506806/1‑150 (MQ=255)
                                atatGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATATCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAg                                                                                                                 >  1:2420104/1‑150 (MQ=255)
                                atatGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAg                                                                                                                 <  2:852599/150‑1 (MQ=255)
                                  atGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGaa                                                                                                               >  1:1403552/1‑150 (MQ=255)
                                          aaCGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTtcttct                                                                                                       <  1:1223939/150‑1 (MQ=255)
                                             gtgtAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGACTTATTGCTGCTTCCAGAATTTCTTCTGaa                                                                                                    >  2:2753819/1‑150 (MQ=255)
                                               gtAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGt                                                                                                  >  2:1153624/1‑150 (MQ=255)
                                                   aagataagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACa                                                                                              <  1:4947989/150‑1 (MQ=255)
                                                   aagataagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACa                                                                                              <  1:2320814/150‑1 (MQ=255)
                                                     gataagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATg                                                                                            >  1:2069867/1‑150 (MQ=255)
                                                     gataagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATg                                                                                            >  1:2070470/1‑150 (MQ=255)
                                                     gataagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATg                                                                                            >  2:3971710/1‑150 (MQ=255)
                                                        aagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCg                                                                                         >  2:399449/1‑150 (MQ=255)
                                                        aagatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCg                                                                                         >  1:2818325/1‑150 (MQ=255)
                                                         agatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGa                                                                                        >  1:2830619/1‑150 (MQ=255)
                                                          gatGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGAt                                                                                       >  2:2255230/1‑150 (MQ=255)
                                                               aaCATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCcgcg                                                                                  >  2:2272111/1‑150 (MQ=255)
                                                               aaCATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCcgcg                                                                                  >  2:554641/1‑150 (MQ=255)
                                                                 cATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGc                                                                                >  1:4651895/1‑150 (MQ=255)
                                                                 cATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAAGTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACAGGTCGATCCGCGGc                                                                                >  2:3206956/1‑150 (MQ=255)
                                                                          aatCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAAt                                                                       >  1:2508774/1‑150 (MQ=255)
                                                                          aatCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAAt                                                                       >  2:1133055/1‑150 (MQ=255)
                                                                          aatCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAAt                                                                       >  2:1132553/1‑150 (MQ=255)
                                                                             cAACATACGATAGTCTAAATTTTACATAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGAc                                                                    >  1:3883387/1‑150 (MQ=255)
                                                                             cAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGAc                                                                    >  2:2491250/1‑150 (MQ=255)
                                                                             cAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGAc                                                                    >  2:2490323/1‑150 (MQ=255)
                                                                               aCATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACtt                                                                  <  1:1397712/150‑1 (MQ=255)
                                                                                  tACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTAt                                                               >  1:2741681/1‑150 (MQ=255)
                                                                                  tACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTAt                                                               >  1:2741608/1‑150 (MQ=255)
                                                                                  tACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTAt                                                               >  1:2741595/1‑150 (MQ=255)
                                                                                       tAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGAtct                                                          >  2:3858201/1‑150 (MQ=255)
                                                                                                ttttACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCaga                                                 <  1:2723624/150‑1 (MQ=255)
                                                                                                      acaATCAGAAAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACCTATATGTCTTTATGATCTCCGGTCAGAGATGGa                                           <  1:214319/150‑1 (MQ=255)
                                                                                                      acaAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGa                                           <  1:2838128/150‑1 (MQ=255)
                                                                                                       caAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGaa                                          <  1:2428683/150‑1 (MQ=255)
                                                                                                       caAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGaa                                          <  1:2429223/150‑1 (MQ=255)
                                                                                                        aaaCAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGaaa                                         >  1:4824041/1‑150 (MQ=255)
                                                                                                        aaaCAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGaaa                                         >  1:5174531/1‑150 (MQ=255)
                                                                                                        aaaCAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGaaa                                         >  1:1283609/1‑150 (MQ=255)
                                                                                                          aCAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAAcc                                       >  2:3839831/1‑150 (MQ=255)
                                                                                                            aGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCt                                     >  1:1777253/1‑150 (MQ=255)
                                                                                                            aGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCt                                     >  1:1777266/1‑150 (MQ=255)
                                                                                                              aCAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGc                                   >  1:175901/1‑150 (MQ=255)
                                                                                                               cAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCg                                  <  1:1166827/150‑1 (MQ=255)
                                                                                                               cAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCg                                  <  1:1166489/150‑1 (MQ=255)
                                                                                                                aaaGAGAATTTTCCTGAATTATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCgg                                 <  1:4577975/150‑1 (MQ=255)
                                                                                                                    agaATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCAt                             >  1:5197072/1‑150 (MQ=255)
                                                                                                                     gaATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCAtt                            <  2:1240199/150‑1 (MQ=255)
                                                                                                                      aaTTTTCCTGAATTATCAATGCAATAGCAACAAATCAACTCAAGAGCCGTATTGCTGCTTCCAGAAATTCTTCTGAAGTAACATTTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGAGGCCATTa                           >  1:3804939/1‑150 (MQ=255)
                                                                                                                           tacTGAATTATCAATTGATTATTATCATTTCATTTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATTTCGATCCGCGGATACTTATATGACTTTATGATCTCCTTTAAGAGATGGAATCCCTGCGGCCATTACAGTa                      <  2:3342762/148‑1 (MQ=255)
                                                                                                                           tCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTa                      >  2:1232939/1‑150 (MQ=255)
                                                                                                                               gAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGgt                  >  1:1451925/1‑150 (MQ=255)
                                                                                                                               gAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGgt                  >  1:1448683/1‑150 (MQ=255)
                                                                                                                               gAATTATCAATGCAATACCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGgt                  >  2:4893948/1‑150 (MQ=255)
                                                                                                                                  ttATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGgtgtg               <  1:2247863/150‑1 (MQ=255)
                                                                                                                                   tATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGgtgtgt              <  2:3506806/150‑1 (MQ=255)
                                                                                                                                   tATCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGgtgtgt              <  2:5107017/150‑1 (MQ=255)
                                                                                                                                    aTCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGtt             >  1:4035654/1‑150 (MQ=255)
                                                                                                                                    aTCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGtt             >  1:4035664/1‑150 (MQ=255)
                                                                                                                                    aTCAATGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGtt             >  1:4036000/1‑150 (MQ=255)
                                                                                                                                     tCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGttt            >  1:1940389/1‑150 (MQ=255)
                                                                                                                                      cAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGtttt           >  1:4065663/1‑150 (MQ=255)
                                                                                                                                       aaTGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGttttt          <  2:436761/150‑1 (MQ=255)
                                                                                                                                       aaTGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGttttt          <  1:755045/150‑1 (MQ=255)
                                                                                                                                       aaTGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGttttt          >  2:1569752/1‑150 (MQ=255)
                                                                                                                                       aaTGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGttttt          >  2:4208175/1‑150 (MQ=255)
                                                                                                                                       aaTGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGttttt          >  2:1569906/1‑150 (MQ=255)
                                                                                                                                         tGCAATAGAATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTc        >  1:1203074/1‑150 (MQ=255)
                                                                                                                                          gCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCg       >  1:254084/1‑150 (MQ=255)
                                                                                                                                            aaTAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCGcc     >  2:3984335/1‑150 (MQ=255)
                                                                                                                                               aGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCGCCAtt  <  1:3206956/150‑1 (MQ=255)
                                                                                                                                               aGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCGCCAtt  >  2:2516961/1‑150 (MQ=255)
                                                                                                                                               aGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCGCCAtt  >  2:2516797/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTATGACTTTCCAGTTACACAGTAAGTCGATTATATGGTGCAAACGTGTAAAAGATAAGATGAAACATCGCAATAATCAACATACGATAGTCTAAATTTTACACAAACAGACAAAGAGAATTTTCCTGAATTATCAATGCAATAGCATCAAATCAACTCAAGAGCCTTATTGCTGCTTCCAGAATTTCTTCTGAAGTAACATGTCGATCCGCGGCTACATAAATGACTTTATGATCTCCGGTCAGAGATGGAAACCCTGCGGCCATTACAGTAAGGTGTGTTTTTTCGCCATT  >  CP009273/1632059‑1632351

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: