Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,228,879 A→C 16.8% L432R (CTG→CGG)  nhaB ← sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,228,8790AC16.8% 784.0 / 167.3 308L432R (CTG→CGG) nhaBsodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (117/139);  new base C (29/23);  total (146/162)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

TAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGATGC  >  CP009273/1228732‑1229027
                                                                                                                                                   |                                                                                                                                                    
tAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCcagcag                                                                                                                                                     >  1:4374382/1‑149 (MQ=255)
  aGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCt                                                                                                                                                   <  2:2843636/149‑1 (MQ=255)
   ggCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGAGTGGTACCGGTATTAATGGCGACCGCCAGCCGCTc                                                                                                                                                  >  2:1221613/1‑149 (MQ=255)
      aGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTa                                                                                                                                               <  1:1676221/149‑1 (MQ=255)
       gggCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTAt                                                                                                                                              <  2:4822549/149‑1 (MQ=255)
          ccATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGc                                                                                                                                           <  1:5111327/149‑1 (MQ=255)
          ccATCCACACCATGCGGCCATAAGAGAGGCGAATCAATAGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGc                                                                                                                                           <  1:2157657/149‑1 (MQ=255)
             tCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTc                                                                                                                                        <  1:174947/149‑1 (MQ=255)
             tCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTc                                                                                                                                        >  2:3397342/1‑149 (MQ=255)
              ccACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCa                                                                                                                                       >  2:3732021/1‑149 (MQ=255)
                  acCATGCGGCCATAAGAGAGGCGAATCACTGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGt                                                                                                                                   <  1:3542116/149‑1 (MQ=255)
                    cATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGa                                                                                                                                 >  2:3314026/1‑149 (MQ=255)
                      tGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATa                                                                                                                               <  1:1044885/149‑1 (MQ=255)
                      tGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATa                                                                                                                               <  2:677698/149‑1 (MQ=255)
                       gCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAg                                                                                                                              <  1:322635/149‑1 (MQ=255)
                              tAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACt                                                                                                                       <  1:2565506/149‑1 (MQ=255)
                              tAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACt                                                                                                                       >  1:3804529/1‑149 (MQ=255)
                               aaGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACtt                                                                                                                      <  2:685616/149‑1 (MQ=255)
                                agagagGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACttt                                                                                                                     >  2:1789550/1‑149 (MQ=255)
                                 gagagGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTc                                                                                                                    <  1:3437705/149‑1 (MQ=255)
                                   gagGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCa                                                                                                                  >  2:43174/1‑149 (MQ=255)
                                   gagGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCa                                                                                                                  >  2:43369/1‑149 (MQ=255)
                                    agGCGAATCAATGGCGCGAGTGCAGAGGTCAGAAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTCCCGGAATTAATGGCGACCGCCAGCAGCTCGTATAGCTTCAACGTGATAGCGCCACTTTCCaa                                                                                                                 >  1:1242490/1‑148 (MQ=255)
                                       cGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGc                                                                                                              <  1:958977/149‑1 (MQ=255)
                                         aaTCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCg                                                                                                            >  2:672904/1‑149 (MQ=255)
                                          aTCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGc                                                                                                           <  1:1062613/149‑1 (MQ=255)
                                          aTCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGAAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGc                                                                                                           <  1:632801/149‑1 (MQ=255)
                                           tCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCt                                                                                                          >  2:2584747/1‑149 (MQ=255)
                                           tCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCt                                                                                                          >  2:2363447/1‑149 (MQ=255)
                                             aaTGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCttt                                                                                                        >  2:4683470/1‑149 (MQ=255)
                                             aaTGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCttt                                                                                                        <  1:651342/149‑1 (MQ=255)
                                              aTGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTc                                                                                                       >  2:3465043/1‑149 (MQ=255)
                                              aTGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTc                                                                                                       >  1:3144621/1‑149 (MQ=255)
                                                   gcgAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTc                                                                                                  <  2:4602685/149‑1 (MQ=255)
                                                           agagGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGatat                                                                                          >  2:4439971/1‑149 (MQ=255)
                                                           agagGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGatat                                                                                          >  2:2312878/1‑149 (MQ=255)
                                                            gagGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGAtata                                                                                         >  1:484021/1‑149 (MQ=255)
                                                              ggTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATaaa                                                                                       <  2:3402310/149‑1 (MQ=255)
                                                              ggTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATaaa                                                                                       <  2:3552260/149‑1 (MQ=255)
                                                               gTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGGTGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAAt                                                                                      <  1:4709864/149‑1 (MQ=255)
                                                               gTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAAt                                                                                      <  1:2817358/149‑1 (MQ=255)
                                                                tCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATc                                                                                     <  1:586907/149‑1 (MQ=255)
                                                                tCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATc                                                                                     <  2:145716/149‑1 (MQ=255)
                                                                  agcagGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGt                                                                                   >  1:3275002/1‑149 (MQ=255)
                                                                   gcagGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTc                                                                                  >  2:3085621/1‑149 (MQ=255)
                                                                   gcagGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATAGta                                                                                  >  1:3987930/1‑148 (MQ=255)
                                                                    caggaacaggaaCGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTcc                                                                                 <  1:2730531/149‑1 (MQ=255)
                                                                    caggaacaggaaCGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTcc                                                                                 <  1:1842583/149‑1 (MQ=255)
                                                                     aggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTccc                                                                                <  1:3752797/149‑1 (MQ=255)
                                                                     aggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTccc                                                                                <  1:1158630/149‑1 (MQ=255)
                                                                     aggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTccc                                                                                >  2:4248130/1‑149 (MQ=255)
                                                                     aggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTccc                                                                                <  2:2211081/149‑1 (MQ=255)
                                                                      ggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTcccc                                                                               >  1:4382175/1‑149 (MQ=255)
                                                                      ggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTcccc                                                                               >  1:2990298/1‑149 (MQ=255)
                                                                      ggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTcccc                                                                               <  2:120881/149‑1 (MQ=255)
                                                                      ggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTcccc                                                                               >  2:4071542/1‑149 (MQ=255)
                                                                      ggaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACAGCCAGCAGCTCGAATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTATCGCTTCGTTGATATAAATCGTcccc                                                                               >  1:2750869/1‑149 (MQ=255)
                                                                       gaacaggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCa                                                                              >  2:3213043/1‑149 (MQ=255)
                                                                           aggaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGaa                                                                          >  2:4916024/1‑149 (MQ=255)
                                                                             gaacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGa                                                                        <  2:1942756/149‑1 (MQ=255)
                                                                              aacGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGAc                                                                       <  1:4608661/149‑1 (MQ=255)
                                                                                cGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGt                                                                     <  2:2422230/149‑1 (MQ=255)
                                                                                 gCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGtt                                                                    <  1:2247797/149‑1 (MQ=255)
                                                                                  cAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTa                                                                   <  1:895225/149‑1 (MQ=255)
                                                                                     ccTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATcc                                                                <  1:450181/149‑1 (MQ=255)
                                                                                         aCCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGaaa                                                            <  1:1302960/149‑1 (MQ=255)
                                                                                         aCCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGaaa                                                            >  1:246838/1‑149 (MQ=255)
                                                                                         aCCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGaaa                                                            >  2:777068/1‑149 (MQ=255)
                                                                                             ttCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGa                                                        >  1:4862304/1‑149 (MQ=255)
                                                                                              tCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGAt                                                       >  2:172224/1‑149 (MQ=255)
                                                                                                     aGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGAcagca                                                >  1:3373189/1‑149 (MQ=255)
                                                                                                        gACGGAGGTCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGAc                                             <  2:4517559/149‑1 (MQ=255)
                                                                                                         aCGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCTTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGAcc                                            >  2:820402/1‑149 (MQ=255)
                                                                                                         aCGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGAcc                                            <  1:196726/149‑1 (MQ=255)
                                                                                                           ggAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGt                                          >  2:1042163/1‑149 (MQ=255)
                                                                                                             aGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTTATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTg                                        >  2:1648614/1‑149 (MQ=255)
                                                                                                                  aGATTGGTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAat                                   >  1:4229614/1‑149 (MQ=255)
                                                                                                                    aTTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAatat                                 <  2:2998859/149‑1 (MQ=255)
                                                                                                                     ttGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAAtata                                >  2:958143/1‑149 (MQ=255)
                                                                                                                       ggTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAACATATAGa                              >  2:2765724/1‑149 (MQ=255)
                                                                                                                       ggTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGa                              >  2:3721501/1‑149 (MQ=255)
                                                                                                                       ggTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGa                              >  2:2765381/1‑149 (MQ=255)
                                                                                                                       ggTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGa                              >  1:36363/1‑149 (MQ=255)
                                                                                                                       ggTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATAAGAAATGGATGACAGCAGACCGTTGAAAATATAGa                              >  1:147221/1‑149 (MQ=255)
                                                                                                                        gTACCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGaa                             >  2:4022068/1‑149 (MQ=255)
                                                                                                                        gTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGaa                             <  2:4375105/149‑1 (MQ=255)
                                                                                                                          aCCGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACa                           <  1:1541776/149‑1 (MQ=255)
                                                                                                                            cGGTATTAATGGCGACCGCCAGCCGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGc                         <  2:473056/149‑1 (MQ=255)
                                                                                                                                       gCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGc              >  2:858376/1‑149 (MQ=255)
                                                                                                                                         gACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCAt            >  1:348358/1‑149 (MQ=255)
                                                                                                                                           ccGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGt          <  2:643068/149‑1 (MQ=255)
                                                                                                                                              ccAGCCGCTCGTATTGCTTCAACGTGATCGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGGCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACCGCGACAGCTGAGCATGTTcc       <  1:3077375/149‑1 (MQ=255)
                                                                                                                                              ccAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTcc       <  1:2313191/149‑1 (MQ=255)
                                                                                                                                              ccAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTcc       <  1:4188447/149‑1 (MQ=255)
                                                                                                                                               cagcagCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCg      >  2:3818571/1‑149 (MQ=255)
                                                                                                                                                agccgcTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGa     >  2:2030778/1‑149 (MQ=255)
                                                                                                                                                 gcagcTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGAt    >  1:760588/1‑149 (MQ=255)
                                                                                                                                                 gcagcTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGAt    <  2:2230588/149‑1 (MQ=255)
                                                                                                                                                  cagcTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGATg   <  2:2886648/149‑1 (MQ=255)
                                                                                                                                                   agcTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGATGc  <  2:2535282/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGCGACCGCCAGCAGCTCGTATTGCTTCAACGTGATAGCGCCACTTTCCATTGCCGCTTTCGCTTCGTTGATATAAATCGTCCCCACGAAGACGTTATCCGAAATGGATGACAGCAGACCGTTGAAAATATAGAACAGCGACAGCTGAGCATGTTCCGATGC  >  CP009273/1228732‑1229027

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: