Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,171,215 C→A 7.0% V111V (GTC→GTA lolC → lipoprotein‑releasing system transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,171,2150CA 7.0% 628.3 / 16.8 170V111V (GTC→GTAlolClipoprotein‑releasing system transmembrane protein
Reads supporting (aligned to +/- strand):  ref base C (76/82);  new base A (5/7);  total (81/89)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 172 total aligned reads displayed.

GATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAATCGCGGTGAT  >  CP009273/1171068‑1171362
                                                                                                                                                   |                                                                                                                                                   
 aTGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCg                                                                                                                                                    >  2:2653244/1‑148 (MQ=255)
    ccACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCgggg                                                                                                                                                 <  1:1357562/148‑1 (MQ=255)
     cacaGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGt                                                                                                                                                >  1:3108814/1‑148 (MQ=255)
      acaGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTg                                                                                                                                               <  2:2721714/148‑1 (MQ=255)
      acaGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTg                                                                                                                                               <  1:2460019/148‑1 (MQ=255)
      acaGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTg                                                                                                                                               <  2:80100/148‑1 (MQ=255)
         ggCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATg                                                                                                                                            <  1:3028192/148‑1 (MQ=255)
           cAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCt                                                                                                                                          <  2:2735242/148‑1 (MQ=255)
           cAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCt                                                                                                                                          <  2:1546666/148‑1 (MQ=255)
           cAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCt                                                                                                                                          <  2:1546689/148‑1 (MQ=255)
                ctctctTCTGAGCATGGCTCTCTTACCCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTa                                                                                                                                     <  1:3149044/148‑1 (MQ=255)
                ctctctTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTa                                                                                                                                     <  2:2171688/148‑1 (MQ=255)
                      tctGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGAcc                                                                                                                               >  1:39848/1‑148 (MQ=255)
                      tctGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGAcc                                                                                                                               >  1:1254816/1‑148 (MQ=255)
                      tctGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGAcc                                                                                                                               >  2:855468/1‑148 (MQ=255)
                      tctGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGAcc                                                                                                                               >  2:855577/1‑148 (MQ=255)
                      tctGAGCAGGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGAcc                                                                                                                               <  1:2311058/148‑1 (MQ=255)
                              tGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCaa                                                                                                                       >  2:1732653/1‑148 (MQ=255)
                                gCTCTCTTACCCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGACCCGGCGCaaaa                                                                                                                     <  1:2792535/148‑1 (MQ=255)
                                gCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGCTGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCaaaa                                                                                                                     >  2:2050144/1‑148 (MQ=255)
                                gCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCaaaa                                                                                                                     <  2:2703951/148‑1 (MQ=255)
                                  tctctTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAg                                                                                                                   >  2:134549/1‑148 (MQ=255)
                                   ctctTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGa                                                                                                                  <  2:402714/148‑1 (MQ=255)
                                    tctTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGAt                                                                                                                 >  2:2212440/1‑148 (MQ=255)
                                     ctTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATc                                                                                                                >  2:1673615/1‑148 (MQ=255)
                                       tAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCa                                                                                                              <  1:2399384/148‑1 (MQ=255)
                                                gcaACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCg                                                                                                     <  1:3142098/148‑1 (MQ=255)
                                                   aCTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTAt                                                                                                  >  2:742265/1‑148 (MQ=255)
                                                    cTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATc                                                                                                 <  2:286788/148‑1 (MQ=255)
                                                    cTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATc                                                                                                 >  1:822597/1‑148 (MQ=255)
                                                     tCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCt                                                                                                >  2:185730/1‑148 (MQ=255)
                                                     tCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCt                                                                                                <  2:1320513/148‑1 (MQ=255)
                                                      cccAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTg                                                                                               >  1:2930960/1‑148 (MQ=255)
                                                       ccAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTgg                                                                                              <  1:1257208/148‑1 (MQ=255)
                                                        cAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGt                                                                                             >  1:2812903/1‑148 (MQ=255)
                                                             aCGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAAtg                                                                                        <  2:2661639/148‑1 (MQ=255)
                                                             aCGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAAtg                                                                                        <  2:445239/148‑1 (MQ=255)
                                                                  aGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGaaa                                                                                   <  1:761503/148‑1 (MQ=255)
                                                                  aGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGaaa                                                                                   <  2:205852/148‑1 (MQ=255)
                                                                   gTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGaaac                                                                                  >  2:2522124/1‑148 (MQ=255)
                                                                   gTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGaaac                                                                                  >  1:2696210/1‑148 (MQ=255)
                                                                    tCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTCGTCAATGTGaaaca                                                                                 >  2:720698/1‑148 (MQ=255)
                                                                        aCTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACt                                                                             >  2:2268374/1‑148 (MQ=255)
                                                                        aCTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACt                                                                             >  1:1543172/1‑148 (MQ=255)
                                                                        aCTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACt                                                                             >  1:1542795/1‑148 (MQ=255)
                                                                        aCTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACt                                                                             >  1:1077561/1‑148 (MQ=255)
                                                                          tGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGa                                                                           <  2:3257734/148‑1 (MQ=255)
                                                                          tGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGa                                                                           <  2:1115202/148‑1 (MQ=255)
                                                                          tGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGa                                                                           <  2:3166795/148‑1 (MQ=255)
                                                                               ggCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCg                                                                      <  1:1579803/148‑1 (MQ=255)
                                                                               ggCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCg                                                                      <  1:2631102/148‑1 (MQ=255)
                                                                                 cGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAg                                                                    <  2:1596983/148‑1 (MQ=255)
                                                                                 cGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAg                                                                    <  1:467778/148‑1 (MQ=255)
                                                                                   ttAATCGCGTCGCACCTGTTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGcc                                                                  >  2:895192/1‑148 (MQ=255)
                                                                                     aaTCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCgg                                                                >  2:2352396/1‑148 (MQ=255)
                                                                                     aaTCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCgg                                                                >  2:454586/1‑148 (MQ=255)
                                                                                     aaTCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCgg                                                                <  1:524166/148‑1 (MQ=255)
                                                                                      aTCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCggg                                                               <  1:3293664/148‑1 (MQ=255)
                                                                                      aTCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCggg                                                               <  1:3293638/148‑1 (MQ=255)
                                                                                      aTCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCggg                                                               >  1:907051/1‑148 (MQ=255)
                                                                                          cgTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGaaa                                                           <  1:2889284/148‑1 (MQ=255)
                                                                                          cgTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGaaa                                                           <  1:864441/148‑1 (MQ=255)
                                                                                          cgTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGaaa                                                           >  1:3188869/1‑148 (MQ=255)
                                                                                           gTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAat                                                          <  2:170317/148‑1 (MQ=255)
                                                                                             cGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAatat                                                        <  1:633403/148‑1 (MQ=255)
                                                                                             cGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAatat                                                        <  1:132747/148‑1 (MQ=255)
                                                                                                 ccTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATg                                                    >  1:1921523/1‑148 (MQ=255)
                                                                                                  cTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGt                                                   >  2:2514602/1‑148 (MQ=255)
                                                                                                   tATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTc                                                  >  2:1571653/1‑148 (MQ=255)
                                                                                                    aTTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCa                                                 >  1:1816285/1‑148 (MQ=255)
                                                                                                     ttACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCAt                                                <  2:292291/148‑1 (MQ=255)
                                                                                                      tactacCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATc                                               <  1:4324/148‑1 (MQ=255)
                                                                                                      tactacCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATc                                               <  2:1559340/148‑1 (MQ=255)
                                                                                                       actacCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATcc                                              <  1:2816106/148‑1 (MQ=255)
                                                                                                       actacCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATcc                                              <  2:3070104/148‑1 (MQ=255)
                                                                                                          acCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCg                                           >  2:939537/1‑148 (MQ=255)
                                                                                                             ggTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCg                                        >  1:2221758/1‑148 (MQ=255)
                                                                                                                 aTGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGaaca                                    >  2:1321466/1‑148 (MQ=255)
                                                                                                                 aTGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGaaca                                    <  2:111592/148‑1 (MQ=255)
                                                                                                                 aTGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGaaca                                    >  1:44435/1‑148 (MQ=255)
                                                                                                                      gTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTg                               >  1:332650/1‑148 (MQ=255)
                                                                                                                       tACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGc                              <  2:822597/148‑1 (MQ=255)
                                                                                                                          tGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTc                           >  2:2754091/1‑148 (MQ=255)
                                                                                                                            cAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTcac                         <  1:1694509/148‑1 (MQ=255)
                                                                                                                              aaGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAg                       <  2:1463317/148‑1 (MQ=255)
                                                                                                                               aGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGc                      >  2:2967856/1‑148 (MQ=255)
                                                                                                                               aGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGc                      >  2:314891/1‑148 (MQ=255)
                                                                                                                               aGCGCGCGCAGCGTGGCGGTAGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGc                      >  1:1086760/1‑148 (MQ=255)
                                                                                                                                        aGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTa             <  2:420197/148‑1 (MQ=255)
                                                                                                                                         gCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGTCGAACAACTTGACTCAAAGCTAGGCGTTaa            <  2:2972346/148‑1 (MQ=255)
                                                                                                                                         gCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTaa            <  1:1685267/148‑1 (MQ=255)
                                                                                                                                          cGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGACATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAAt           >  2:1868691/1‑148 (MQ=255)
                                                                                                                                          cGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAAt           <  1:269261/148‑1 (MQ=255)
                                                                                                                                          cGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAAt           >  1:524240/1‑148 (MQ=255)
                                                                                                                                          cGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAAt           >  2:1032675/1‑148 (MQ=255)
                                                                                                                                          cGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAAt           >  1:1265891/1‑148 (MQ=255)
                                                                                                                                            tGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAATcg         >  1:207837/1‑148 (MQ=255)
                                                                                                                                                ggTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTGGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAATCGCGGt     <  1:2653244/148‑1 (MQ=255)
                                                                                                                                                  tCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAATCGCGGTGa   <  1:2212440/148‑1 (MQ=255)
                                                                                                                                                  tCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAATCGCGGTGa   <  1:2052275/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTGGCGGTCGGGGTGATGCTCGGTATCGACCCGGCGCAAAAAGATCCACTTACACCGTATCTGGTCAATGTGAAACAAACTGACCTCGAGCCGGGGAAATATAATGTCATCCTCGGCGAACAACTTGCCTCACAGCTAGGCGTTAATCGCGGTGAT  >  CP009273/1171068‑1171362

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: