Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,075,114 C→A 94.3% intergenic (‑120/‑158) tktA ← / → loiP transketolase 1, thiamine triphosphate‑binding/Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; heat shock protein that binds Era

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,075,1140CA94.3% 910.6 / 39.2 280intergenic (‑120/‑158)tktA/loiPtransketolase 1, thiamine triphosphate‑binding/Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; heat shock protein that binds Era
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base A (115/149);  total (122/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 280 total aligned reads displayed.

CATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCG  >  CP009273/3074967‑3075260
                                                                                                                                                   |                                                                                                                                                  
cATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  >  1:2995178/1‑150 (MQ=255)
cATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  <  2:4320517/150‑1 (MQ=255)
cATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  >  1:2740219/1‑150 (MQ=255)
cATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  >  1:2018733/1‑150 (MQ=255)
 aTTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATc                                                                                                                                                 <  2:3439698/150‑1 (MQ=255)
  ttGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATcc                                                                                                                                                >  2:529919/1‑150 (MQ=255)
    ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCtt                                                                                                                                              >  2:1904679/1‑150 (MQ=255)
       aaGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTg                                                                                                                                           >  2:248998/1‑150 (MQ=255)
       aaGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTg                                                                                                                                           >  1:3866559/1‑150 (MQ=255)
          ctctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCg                                                                                                                                        <  2:3637764/150‑1 (MQ=255)
           tctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGc                                                                                                                                       >  2:3857299/1‑150 (MQ=255)
            ctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCt                                                                                                                                      >  1:2912709/1‑150 (MQ=255)
             tttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCtt                                                                                                                                     >  2:151388/1‑150 (MQ=255)
                aCGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAg                                                                                                                                  <  1:3778669/150‑1 (MQ=255)
                aCGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAg                                                                                                                                  <  2:2381709/150‑1 (MQ=255)
                 cGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGa                                                                                                                                 >  2:827372/1‑150 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaa                                                                                                                                >  1:2365114/1‑150 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaa                                                                                                                                <  2:1949829/150‑1 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTAGCCGCTTCAGaa                                                                                                                                >  1:1307429/1‑150 (MQ=255)
                   tGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaaa                                                                                                                               >  1:1397241/1‑150 (MQ=255)
                   tGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaaa                                                                                                                               >  2:2112272/1‑150 (MQ=255)
                     aGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtc                                                                                                                             >  2:3094356/1‑150 (MQ=255)
                     aGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtc                                                                                                                             >  1:212333/1‑150 (MQ=255)
                     aGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtc                                                                                                                             <  2:3998602/150‑1 (MQ=255)
                       gACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACTGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtctc                                                                                                                           <  2:4096685/150‑1 (MQ=255)
                       gACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtctc                                                                                                                           <  2:220841/150‑1 (MQ=255)
                         cATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTg                                                                                                                         <  1:1917912/150‑1 (MQ=255)
                         cATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTg                                                                                                                         <  2:272946/150‑1 (MQ=255)
                            tttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGAt                                                                                                                      <  2:765047/150‑1 (MQ=255)
                            tttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGAt                                                                                                                      <  2:2006584/150‑1 (MQ=255)
                             ttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATc                                                                                                                     <  2:266816/150‑1 (MQ=255)
                             ttGACTCCAGATCGGATGAGGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATc                                                                                                                     <  2:2076308/150‑1 (MQ=255)
                                   ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATATCTGGATCATGCTc                                                                                                               <  1:1027407/150‑1 (MQ=255)
                                      gATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCa                                                                                                            >  1:1991349/1‑150 (MQ=255)
                                       aTCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCAt                                                                                                           >  1:3701547/1‑150 (MQ=255)
                                             tgatgaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTgcg                                                                                                     >  2:1374149/1‑150 (MQ=255)
                                               atgaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCa                                                                                                   >  1:2555949/1‑150 (MQ=255)
                                                  aaGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATc                                                                                                >  2:3957028/1‑150 (MQ=255)
                                                     ggCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTAc                                                                                             <  1:3545635/150‑1 (MQ=255)
                                                     ggCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTAc                                                                                             >  2:3189200/1‑150 (MQ=255)
                                                       cACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTc                                                                                           <  1:1945981/150‑1 (MQ=255)
                                                        aCGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCg                                                                                          <  1:3490581/150‑1 (MQ=255)
                                                         cGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGc                                                                                         >  1:1308883/1‑150 (MQ=255)
                                                             cTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGt                                                                                     >  1:3251743/1‑150 (MQ=255)
                                                              ttAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTc                                                                                    <  1:2454877/150‑1 (MQ=255)
                                                                aaCGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCg                                                                                  <  1:1840264/150‑1 (MQ=255)
                                                                   gACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTg                                                                               >  2:1121681/1‑150 (MQ=255)
                                                                    aCTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGTCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTgc                                                                              <  2:4266112/150‑1 (MQ=255)
                                                                     cTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTgcg                                                                             <  1:3304557/150‑1 (MQ=255)
                                                                     cTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTgcg                                                                             >  2:3896098/1‑150 (MQ=255)
                                                                         acgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCttt                                                                         <  2:3330929/150‑1 (MQ=255)
                                                                         acgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCttt                                                                         <  1:915877/150‑1 (MQ=255)
                                                                          cgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCtttt                                                                        <  2:3237607/150‑1 (MQ=255)
                                                                          cgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCtttt                                                                        <  1:1090685/150‑1 (MQ=255)
                                                                            acAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTcc                                                                      <  2:864598/150‑1 (MQ=255)
                                                                             cAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCt                                                                     <  2:732499/150‑1 (MQ=255)
                                                                             cAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCt                                                                     >  1:3959892/1‑150 (MQ=255)
                                                                                cgcgTTTTGGGCTACGCCTGAAAATTTGCCAACAATTGACCGCCAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCTCATGTTTCGCAATATACTCGCCCGTCCGCTGCGCTTTTCATtat                                                                  <  2:901572/150‑1 (MQ=255)
                                                                                cgcgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTtat                                                                  >  2:1221438/1‑150 (MQ=255)
                                                                                  cgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATac                                                                <  1:1754115/150‑1 (MQ=255)
                                                                                  cgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATac                                                                <  1:3638982/150‑1 (MQ=255)
                                                                                  cgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATac                                                                <  1:542529/150‑1 (MQ=255)
                                                                                   gTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATact                                                               <  1:2588156/150‑1 (MQ=255)
                                                                                      ttGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgag                                                            <  2:3419845/150‑1 (MQ=255)
                                                                                         ggCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagact                                                         <  2:1849885/150‑1 (MQ=255)
                                                                                           cTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagactga                                                       <  2:683667/150‑1 (MQ=255)
                                                                                            tACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagactgag                                                      >  2:2954206/1‑150 (MQ=255)
                                                                                                ccGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCGTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTc                                                  <  2:2828030/150‑1 (MQ=255)
                                                                                                 cGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCGGCGCTTTTCCTTATACTGAGACTGAGCGTCg                                                 <  2:2035187/150‑1 (MQ=255)
                                                                                                  ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGa                                                <  1:4149538/150‑1 (MQ=255)
                                                                                                   gAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGAt                                               >  1:4327506/1‑150 (MQ=255)
                                                                                                   gAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGAt                                               <  1:673498/150‑1 (MQ=255)
                                                                                                    aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGAtt                                              <  1:3215018/150‑1 (MQ=255)
                                                                                                    aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGAtt                                              >  1:313902/1‑150 (MQ=255)
                                                                                                     aaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTc                                             <  2:3684033/150‑1 (MQ=255)
                                                                                                       aTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCAc                                           >  1:1387248/1‑150 (MQ=255)
                                                                                                       aTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCAc                                           <  2:2363049/150‑1 (MQ=255)
                                                                                                         ttGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCt                                         >  1:1117015/1‑150 (MQ=255)
                                                                                                            ccAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCa                                      >  1:1444128/1‑150 (MQ=255)
                                                                                                             caacaaTGTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCaa                                     <  2:3582065/150‑1 (MQ=255)
                                                                                                                caaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACg                                  <  1:390291/150‑1 (MQ=255)
                                                                                                                 aaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACgg                                 <  1:859680/150‑1 (MQ=255)
                                                                                                                 aaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACgg                                 <  2:2055142/150‑1 (MQ=255)
                                                                                                                  aTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGgc                                >  2:2934870/1‑150 (MQ=255)
                                                                                                                        cGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAttt                          >  1:1718222/1‑150 (MQ=255)
                                                                                                                          cAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAttttt                        <  1:3856163/150‑1 (MQ=255)
                                                                                                                                 cgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaataa                 >  1:1666287/1‑150 (MQ=255)
                                                                                                                                  gcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaataat                >  1:3352674/1‑150 (MQ=255)
                                                                                                                                   cgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATc               <  2:3883490/150‑1 (MQ=255)
                                                                                                                                    gTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATcc              <  1:1750148/150‑1 (MQ=255)
                                                                                                                                     tCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCt             >  1:2224441/1‑150 (MQ=255)
                                                                                                                                     tCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCt             <  1:2293262/150‑1 (MQ=255)
                                                                                                                                       aTGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGa           >  2:3824924/1‑150 (MQ=255)
                                                                                                                                       aTGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGa           <  2:14485/150‑1 (MQ=255)
                                                                                                                                       aTGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGa           >  2:1636466/1‑150 (MQ=255)
                                                                                                                                            cATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtg      >  2:4169026/1‑150 (MQ=255)
                                                                                                                                            cATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtg      >  2:3055194/1‑150 (MQ=255)
                                                                                                                                             aTGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtgt     <  1:4079770/150‑1 (MQ=255)
                                                                                                                                             aTGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtgt     <  1:4242607/150‑1 (MQ=255)
                                                                                                                                                gAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCg  >  1:3836199/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCG  >  CP009273/3074967‑3075260

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: