Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,177 C→G 13.4% intergenic (+42/‑165) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,068,1770CG13.4% 668.2 / 77.7 222intergenic (+42/‑165)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (77/115);  new base G (8/21);  total (85/137)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 232 total aligned reads displayed.

AGCCCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACT  >  NC_002947/2068028‑2068326
                                                                                                                                                     |                                                                                                                                                      
aGCCCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTcc                                                                                                                                                        <  1:6771821/150‑1 (MQ=255)
aGCCCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTcc                                                                                                                                                        <  1:3332523/150‑1 (MQ=255)
  cccGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGca                                                                                                                                                      <  1:19254/150‑1 (MQ=255)
  cccGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCca                                                                                                                                                      <  1:6548029/150‑1 (MQ=255)
   ccGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCcac                                                                                                                                                     <  1:5707090/150‑1 (MQ=255)
    cGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                                    >  2:2800435/1‑150 (MQ=255)
     ggACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAg                                                                                                                                                   <  1:1029589/150‑1 (MQ=255)
       aCCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGc                                                                                                                                                 <  1:4459867/150‑1 (MQ=255)
         cAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAggacc                                                                                                                                               >  2:5805128/1‑147 (MQ=255)
          aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTGGCCGGGGTGGCCGCGGCCACCgggctg                                                                                                                                              >  1:2528685/1‑146 (MQ=255)
          aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                              <  2:4259911/150‑1 (MQ=255)
          aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                              >  2:5323741/1‑150 (MQ=255)
               gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCTCCCCCACGCATTTcaca                                                                                                                                         >  2:1518411/1‑150 (MQ=255)
                aaaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAg                                                                                                                                        <  2:1814940/150‑1 (MQ=255)
                 aaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAgg                                                                                                                                       <  1:3052834/150‑1 (MQ=255)
                   tGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGt                                                                                                                                     >  2:1437731/1‑150 (MQ=255)
                    gCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGtt                                                                                                                                    <  2:7987904/150‑1 (MQ=255)
                     cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                                   <  2:2512785/150‑1 (MQ=255)
                      tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGttag                                                                                                                                  >  2:6004274/1‑148 (MQ=255)
                        aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                               >  2:2645690/1‑150 (MQ=255)
                        aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                               >  2:5458152/1‑150 (MQ=255)
                        aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAGTGCACAGGGTTCGGa                                                                                                                               >  2:5786568/1‑150 (MQ=255)
                        aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCGCCCACAGGCATTGCACAGGGGTCGGa                                                                                                                               >  1:3612747/1‑150 (MQ=255)
                         agaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAc                                                                                                                              <  1:5323741/150‑1 (MQ=255)
                          gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAATGCACCGGGTTCGGAcc                                                                                                                             >  1:4368085/1‑150 (MQ=255)
                            agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGTCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                           <  2:3290798/150‑1 (MQ=255)
                            agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                           <  2:5134945/150‑1 (MQ=255)
                            agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGTCCGAAccct                                                                                                                           <  2:8271627/150‑3 (MQ=255)
                               tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtgcac                                                                                                                        <  1:645003/150‑5 (MQ=255)
                                ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGGGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGATTGCCCGCTCCCACAGGCATTGCACAGGGTTAAGACCTTGTGc                                                                                                                       <  1:6153532/150‑1 (MQ=255)
                                ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                       <  2:4618850/150‑1 (MQ=255)
                                  aaGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCaa                                                                                                                     <  2:6199610/150‑2 (MQ=255)
                                     cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCCCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                  <  1:2410129/150‑5 (MQ=25)
                                       tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                                <  2:1687258/150‑7 (MQ=21)
                                       tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcct                                                                                                                <  2:237520/150‑7 (MQ=21)
                                       tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                <  1:1449579/150‑1 (MQ=255)
                                        atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCGTGTGTAGGTCCtg                                                                                                               >  2:2031176/1‑150 (MQ=255)
                                         tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCTGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                              <  1:2883506/150‑1 (MQ=255)
                                         tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                              <  2:4596661/150‑1 (MQ=255)
                                           tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                            >  2:7513104/1‑150 (MQ=255)
                                           tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                            >  2:229969/1‑150 (MQ=255)
                                           tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                            <  1:5458152/150‑1 (MQ=255)
                                           tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                            <  2:7898340/150‑1 (MQ=255)
                                            ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAGTGCACCGGGTTCGGACCATGTGCAGGTCCAGTGCg                                                                                                           >  2:2568704/1‑150 (MQ=255)
                                              aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAg                                                                                                         >  2:3541558/1‑150 (MQ=255)
                                              aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAg                                                                                                         >  2:6405710/1‑150 (MQ=255)
                                                 ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                      <  2:3946746/150‑1 (MQ=18)
                                                 ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                      <  1:7244935/150‑1 (MQ=255)
                                                 ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                      <  2:1145383/150‑1 (MQ=255)
                                                   gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGc                                                                                                    <  1:5951702/150‑1 (MQ=18)
                                                               ctgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGAAAGCCCGCGAAt                                                                                        <  2:7189743/150‑1 (MQ=255)
                                                                tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCAAATg                                                                                       <  2:8530094/150‑1 (MQ=17)
                                                                tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                       <  2:5051006/150‑1 (MQ=255)
                                                                tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAATGGCGGGGGTGGCCCCTCCCACAGGCCTTGCAAACGGTTCGGACCATTTGCAGGTACAGTGCGAGCGGGCAAGCCCGCGCATg                                                                                       >  1:689592/1‑150 (MQ=255)
                                                                    ggATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGGTTGCCCGCTCCCCAAGGCAATGCACAGGGTTCGGACCATTTGCAGGTCCTGTGCGAGCGGGCAAGCCCCCGAATGGGcc                                                                                   >  2:3196477/1‑150 (MQ=255)
                                                                           acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGGGTGGCCGCGCCCCCCCGCCATGCCCACGGGTCGGACCGTGTGGCGGGCCGGGGCGAGCGGGCCAGCCCCCGCAGGGGccccccccc                                                                            >  1:7162983/1‑145 (MQ=255)
                                                                                       gtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGAc                                                                <  2:5862313/150‑1 (MQ=255)
                                                                                          cACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCGCCCACAGGCATTTCACAGGGTTCGGACCATTTGCAGGTCCAGTGCGAGCGGGCCACCCCCCCAATGGGGCCCCCACAAATTTACCCGACGGa                                                             >  2:2994438/1‑150 (MQ=255)
                                                                                               cacaCCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCTTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACAtt                                                        <  2:6940620/150‑1 (MQ=255)
                                                                                                    cGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATccc                                                   <  2:1221571/150‑1 (MQ=255)
                                                                                                            gaatgggCCGCAACGCGGCCCCTTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCCCCGGGTTCGGACCTTGGGCAGGTCCGGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGCCGGAACAATATCccccccccca                                           >  2:4924767/5‑148 (MQ=25)
                                                                                                            aCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                           <  2:5633478/150‑1 (MQ=255)
                                                                                                             cgggggCCGGGTGGCGGCCCTTTCGGGGGTTGGCCCGCTCCAACTGGATTTGCACTGGTTTGGGACCTTGTGCCGGCCCTGTGCGAGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAg                                          <  1:4965442/148‑1 (MQ=255)
                                                                                                              atgggCCGCAACGCGGCCCATTCGCGGGCGTGCCCGCTCCCACAGGCAGTGCACAGGGGTCGGTCCATGTGCAGGTCCAGTGCGAGCGGGCAAGCCCGCGGATGGGCCCCGCCCAAATTTACCGGGCAGTACATTATCCCCGCCCCCAgg                                         >  1:8001131/3‑150 (MQ=255)
                                                                                                                  gCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCAAAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAg                                     <  2:6541715/150‑1 (MQ=255)
                                                                                                                    cgcgTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGtt                                   <  2:7140797/150‑1 (MQ=255)
                                                                                                                       gTTGGGGCCCTTTCGGGGGCTTGCCCTCCCCCACAGCCATTGCCCCGGGTTCGCAGCTGGTGCCGGCCGGGGGCGCGCGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                <  1:7503721/150‑1 (MQ=255)
                                                                                                                       aacgCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                         <  1:1953149/58‑1 (MQ=255)
                                                                                                                       aacgCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                         >  2:1953149/4‑61 (MQ=255)
                                                                                                                        ttGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATGGCACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                               <  2:7469350/150‑1 (MQ=255)
                                                                                                                         tGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCGTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                              <  2:3138303/150‑1 (MQ=255)
                                                                                                                         tGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                              <  2:4330367/150‑1 (MQ=255)
                                                                                                                         ggCGGCCCTTCCGCGGGCTGGCCCGCTCCCCCCGGCTTGGCACCGGGTTCGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                              <  1:5757024/149‑1 (MQ=255)
                                                                                                                               ccATTCGCGGGCTTGCCCGCTCCCACGGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACCTTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                        >  2:1109426/1‑150 (MQ=255)
                                                                                                                               ccATTCGCGGGCTTGCCCGCACCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAACCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCCGTTATTCAGGCATc                        >  2:2360235/1‑150 (MQ=255)
                                                                                                                                cttgcgggggCTTGCCCGCTCCCCCCGGCATTCCAGTGGGTTCGGACCTTGTGCCGGCCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                       <  1:4293497/143‑1 (MQ=255)
                                                                                                                                       gggCTTGCCCGCTCCCACAGGCATTGCACAGGGTTGGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCt                <  1:5789446/150‑1 (MQ=255)
                                                                                                                                          ggtGCCCGCGCCCACCGGCTTTGCACCGGGTTCGGACGTTGTGCAGGTCCTGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt             <  1:7911451/148‑1 (MQ=255)
                                                                                                                                          cTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt             >  2:764665/1‑150 (MQ=255)
                                                                                                                                           ttGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg            >  1:5219093/1‑150 (MQ=18)
                                                                                                                                           ttGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg            >  1:4535385/1‑150 (MQ=255)
                                                                                                                                             gcccgcccCCCCGGGCATGGCCGTGGGTTCGGCCCTGGTGCTGGCCGGGTGCGCGGGGCCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGtt          <  1:3834770/150‑1 (MQ=255)
                                                                                                                                             gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGtt          >  2:793279/1‑150 (MQ=255)
                                                                                                                                             gCCCGCACCCACCCGCATTGCACGGGTTCGGCACCTTGTGCCGGGCCTGTGCGCGGGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGtt          <  1:1763152/150‑1 (MQ=255)
                                                                                                                                              cccGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc         <  2:2464958/150‑1 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACCCACCTGGTGTTCg        >  1:855484/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  2:1898715/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  2:2928487/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:7184485/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:6345888/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:364650/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:3552378/1‑150 (MQ=255)
                                                                                                                                                cGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt       >  1:2467665/1‑150 (MQ=18)
                                                                                                                                                cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt       >  1:2593744/1‑150 (MQ=255)
                                                                                                                                                cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGGTCGt       >  2:8402218/1‑150 (MQ=255)
                                                                                                                                                     ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  2:6644897/1‑150 (MQ=255)
                                                                                                                                                     ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  <  2:3430547/150‑1 (MQ=255)
                                                                                                                                                     ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  2:1772040/1‑150 (MQ=255)
                                                                                                                                                     ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  1:1413376/1‑150 (MQ=255)
                                                                                                                                                     ccACAGGCATTGAACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCAGCGAAAGGGCCCCGCACAATTTTACCGGACTGAACATTATCCCCGCCCCAATGTCAGTTATTCAGGCATCACACACCTGGTGTTAGTTGACt  >  2:4292198/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCCCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACT  >  NC_002947/2068028‑2068326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: