Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,075,114 C→A 85.7% 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,1140CA85.7% 744.1 / 127.4 287intergenic (‑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 (22/19);  new base A (133/112);  total (156/131)
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.89e-01

Only 100 of 288 total aligned reads displayed.

CATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAA  >  CP009273/3074967‑3075263
                                                                                                                                                   |                                                                                                                                                     
 aTTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATc                                                                                                                                                    <  1:3163356/150‑1 (MQ=255)
   tGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCt                                                                                                                                                  >  2:3767164/1‑150 (MQ=255)
    ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCtt                                                                                                                                                 >  2:4203205/1‑150 (MQ=255)
    ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtt                                                                                                                                                 <  1:4012768/150‑1 (MQ=255)
    ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtt                                                                                                                                                 <  1:2822009/150‑1 (MQ=255)
      cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtttt                                                                                                                                               >  2:654591/1‑150 (MQ=255)
      cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtttt                                                                                                                                               <  2:2397641/150‑1 (MQ=255)
      cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtttt                                                                                                                                               >  1:1489218/1‑150 (MQ=255)
      cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtttt                                                                                                                                               <  2:138451/150‑1 (MQ=255)
      cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtgtt                                                                                                                                               >  2:2351725/1‑150 (MQ=255)
       aaGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTg                                                                                                                                              >  2:946641/1‑150 (MQ=255)
        aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGc                                                                                                                                             >  2:1680009/1‑150 (MQ=255)
          ctctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCg                                                                                                                                           <  1:1365318/150‑1 (MQ=255)
          ctctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCg                                                                                                                                           <  1:1753339/150‑1 (MQ=255)
                aCGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAg                                                                                                                                     >  1:164681/1‑150 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGaa                                                                                                                                   <  2:1435343/150‑1 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaa                                                                                                                                   <  2:1324473/150‑1 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaa                                                                                                                                   >  1:3462472/1‑150 (MQ=255)
                       gACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAAtctc                                                                                                                              >  1:3338972/1‑150 (MQ=255)
                       gACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtctc                                                                                                                              >  1:1121665/1‑150 (MQ=255)
                       gACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtctc                                                                                                                              >  2:4068268/1‑150 (MQ=255)
                          aTTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTgg                                                                                                                           <  1:1468395/150‑1 (MQ=255)
                           ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGa                                                                                                                          <  1:3650709/150‑1 (MQ=255)
                            tttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGAt                                                                                                                         <  2:1869003/150‑1 (MQ=255)
                               gACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTGTTGCCGCTTCAGAAATCTCTGGATCAt                                                                                                                      >  1:3558154/1‑150 (MQ=255)
                                aCTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATg                                                                                                                     >  2:2798277/1‑150 (MQ=255)
                                 cTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGc                                                                                                                    >  2:2917467/1‑150 (MQ=255)
                                 cTCCAGATCGGATGATGAAGGGCACGCCCTAAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTAACCGCAAGCCGCGCGTCATGTACATGGAAAATCCATTAGCCGCTTCAGACAACTCTCGATCATGc                                                                                                                    >  1:4406132/1‑150 (MQ=255)
                                  tCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCt                                                                                                                   <  2:3286022/150‑1 (MQ=255)
                                   ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTc                                                                                                                  >  2:3183184/1‑150 (MQ=255)
                                   ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTc                                                                                                                  <  2:4490611/150‑1 (MQ=255)
                                   ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTc                                                                                                                  <  2:2111325/150‑1 (MQ=255)
                                     aGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGc                                                                                                                <  1:296573/150‑1 (MQ=255)
                                        tCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATg                                                                                                             <  2:3019295/150‑1 (MQ=255)
                                         cGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGt                                                                                                            >  2:2071525/1‑150 (MQ=255)
                                              gatgaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTgcgc                                                                                                       >  1:746979/1‑150 (MQ=255)
                                                 gaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAAt                                                                                                    >  2:823023/1‑150 (MQ=255)
                                                  aaGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATc                                                                                                   >  1:3095794/1‑150 (MQ=255)
                                                   aGGGAACGCACTTAACTACTTGACGACATCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCTCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCt                                                                                                  <  2:2719645/150‑1 (MQ=255)
                                                    gggCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTa                                                                                                 <  2:3998049/150‑1 (MQ=255)
                                                    gggCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTa                                                                                                 <  2:1267648/150‑1 (MQ=255)
                                                               tAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTcc                                                                                      <  2:1489372/150‑1 (MQ=255)
                                                                 aCGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGc                                                                                    >  1:2558698/1‑150 (MQ=255)
                                                                   gACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTg                                                                                  >  1:2245415/1‑150 (MQ=255)
                                                                        gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCtt                                                                             >  2:300421/1‑150 (MQ=255)
                                                                        gacgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCtt                                                                             <  2:358932/150‑1 (MQ=255)
                                                                         acgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCttt                                                                            >  1:4376786/1‑150 (MQ=255)
                                                                         acgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCttt                                                                            <  1:1262490/150‑1 (MQ=255)
                                                                           gacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTc                                                                          <  2:314702/150‑1 (MQ=255)
                                                                           gacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTc                                                                          >  2:2412141/1‑150 (MQ=255)
                                                                              aGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCtt                                                                       >  2:4455249/1‑150 (MQ=255)
                                                                               gcgcGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCtca                                                                      <  1:4054464/150‑3 (MQ=255)
                                                                                cgcgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTtat                                                                     >  2:1614606/1‑150 (MQ=255)
                                                                                  cgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATact                                                                  >  2:4497841/1‑150 (MQ=255)
                                                                                   gTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATact                                                                  <  2:3570432/150‑1 (MQ=255)
                                                                                   gTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATact                                                                  >  1:1206839/1‑150 (MQ=255)
                                                                                    ttttGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactg                                                                 <  1:2798514/150‑1 (MQ=255)
                                                                                       tGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgaga                                                              <  2:535942/150‑1 (MQ=255)
                                                                                       tGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgaga                                                              >  1:688207/1‑150 (MQ=255)
                                                                                        gggCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagac                                                             <  1:2125926/150‑1 (MQ=255)
                                                                                         ggCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagact                                                            <  1:827803/150‑1 (MQ=255)
                                                                                         ggCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGATGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagact                                                            >  2:547410/1‑150 (MQ=255)
                                                                                          gCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagactg                                                           <  2:1276883/150‑1 (MQ=255)
                                                                                           cTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagactga                                                          >  1:405685/1‑150 (MQ=255)
                                                                                            tACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagactgag                                                         >  2:2479915/1‑150 (MQ=255)
                                                                                             aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGc                                                        <  2:1555199/150‑1 (MQ=255)
                                                                                             aCGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGc                                                        <  1:1239628/150‑1 (MQ=255)
                                                                                               gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGt                                                      >  2:4053830/1‑150 (MQ=255)
                                                                                                ccGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTc                                                     >  2:2157340/1‑150 (MQ=255)
                                                                                                 cGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCg                                                    >  1:288719/1‑150 (MQ=255)
                                                                                                  ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGa                                                   >  1:3528919/1‑150 (MQ=255)
                                                                                                      aaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCa                                               >  2:2561026/1‑150 (MQ=255)
                                                                                                      aaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCa                                               <  2:1733912/150‑1 (MQ=255)
                                                                                                         ttGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCt                                            <  1:963726/150‑1 (MQ=255)
                                                                                                          tGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTg                                           <  2:489405/150‑1 (MQ=255)
                                                                                                            ccAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCa                                         >  2:1076392/1‑150 (MQ=255)
                                                                                                             caacaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCaa                                        <  1:4233174/150‑1 (MQ=255)
                                                                                                             caacaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCaa                                        <  1:1600725/150‑1 (MQ=255)
                                                                                                              aacaaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCaaa                                       >  2:2936865/1‑150 (MQ=255)
                                                                                                              aacaaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCaaa                                       <  2:3515648/150‑1 (MQ=255)
                                                                                                               acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAAc                                      <  2:457003/150‑1 (MQ=255)
                                                                                                                 aaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACgg                                    >  2:2859828/1‑150 (MQ=255)
                                                                                                                  aTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGgc                                   >  2:3827108/1‑150 (MQ=255)
                                                                                                                   tttACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGgcg                                  >  1:3404500/1‑150 (MQ=255)
                                                                                                                    ttACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGgcgc                                 >  1:1694161/1‑150 (MQ=255)
                                                                                                                     tACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCa                                <  1:578517/150‑1 (MQ=255)
                                                                                                                      aCCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAt                               >  1:205867/1‑150 (MQ=255)
                                                                                                                        cGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAttt                             >  1:3863251/1‑150 (MQ=255)
                                                                                                                         gCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAtttt                            >  2:2209945/1‑150 (MQ=255)
                                                                                                                           aaGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTa                          >  2:4299124/1‑150 (MQ=255)
                                                                                                                            aGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAg                         <  1:4466637/150‑1 (MQ=255)
                                                                                                                              ccGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaa                       >  1:471256/1‑150 (MQ=255)
                                                                                                                                 cgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaataa                    <  1:420353/150‑1 (MQ=255)
                                                                                                                                    gTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATcc                 >  1:3707548/1‑150 (MQ=255)
                                                                                                                                    gTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATcc                 <  2:2034451/150‑1 (MQ=255)
                                                                                                                                        tGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGAc             <  1:3549775/150‑1 (MQ=255)
                                                                                                                                         gTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGAcc            >  1:4079126/1‑150 (MQ=255)
                                                                                                                                          tACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCt           <  2:4159904/150‑1 (MQ=255)
                                                                                                                                              tGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtgtg       >  2:44144/1‑150 (MQ=255)
                                                                                                                                                gAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCg     >  1:2667989/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAA  >  CP009273/3074967‑3075263

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: