Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,185 A→C 22.7% intergenic (+50/‑157) 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,1850AC22.7% 322.1 / 95.5 166intergenic (+50/‑157)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (51/71);  new base C (19/19);  total (70/96)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

Only 100 of 169 total aligned reads displayed.

CCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGT  >  NC_002947/2068036‑2068334
                                                                                                                                                     |                                                                                                                                                     
ccAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCa                                                                                                                                                       <  1:252930/150‑1 (MQ=255)
 cAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAt                                                                                                                                                      >  2:4095360/1‑150 (MQ=255)
 cAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCAt                                                                                                                                                      >  1:2767432/1‑150 (MQ=255)
  aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                                     <  2:6323160/150‑1 (MQ=255)
   gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTg                                                                                                                                                    <  1:1056972/150‑1 (MQ=255)
   gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTg                                                                                                                                                    <  1:3268745/150‑1 (MQ=255)
    ccGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGc                                                                                                                                                   <  1:377509/150‑1 (MQ=255)
     cGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGca                                                                                                                                                  >  2:4701286/1‑150 (MQ=255)
     cGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGca                                                                                                                                                  <  1:6779789/150‑1 (MQ=255)
      ggAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcac                                                                                                                                                 >  2:1101861/1‑150 (MQ=255)
       gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                                <  2:6615046/150‑1 (MQ=255)
        aaaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAg                                                                                                                                               <  2:2058776/150‑1 (MQ=255)
          aTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAggg                                                                                                                                             <  2:1051349/150‑1 (MQ=255)
             cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTc                                                                                                                                          >  2:310193/1‑150 (MQ=255)
             cTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                                          <  2:3797395/150‑1 (MQ=255)
              tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCg                                                                                                                                         >  2:412542/1‑150 (MQ=255)
               cAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCCATGCACAGGGTTCgg                                                                                                                                        >  1:5138269/1‑150 (MQ=255)
                aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                                       >  2:1532148/1‑150 (MQ=255)
                 agaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAc                                                                                                                                      <  1:6261221/150‑1 (MQ=255)
                  gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAcc                                                                                                                                     <  2:2554120/150‑1 (MQ=255)
                  gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAcc                                                                                                                                     >  2:5914944/1‑150 (MQ=255)
                    agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAAccct                                                                                                                                   <  2:3694803/150‑3 (MQ=37)
                    agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAAccct                                                                                                                                   <  1:4365487/150‑3 (MQ=37)
                    agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCCCCCACAGGCATTGCACAGGGGTCGGACCtt                                                                                                                                   >  1:1069959/1‑150 (MQ=255)
                      cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGAGTGGTGCACGGCACACCGAAGTATACCGGTGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCtgt                                                                                                                                 <  2:333629/150‑1 (MQ=25)
                       tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCtgtg                                                                                                                                <  1:2257563/150‑1 (MQ=37)
                       tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                <  2:2410359/150‑1 (MQ=255)
                       tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCTCCCACAGGCCTTGCACCGGGGTCGGACCGtgtg                                                                                                                                >  1:1658274/1‑150 (MQ=255)
                          aaGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCACTGCACAAGGTCCGAACCCTGTGCaa                                                                                                                             <  1:552988/150‑2 (MQ=38)
                           aGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAgg                                                                                                                            <  2:1825595/150‑1 (MQ=255)
                             cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                          <  2:2348951/150‑5 (MQ=25)
                             cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                          <  2:1318076/150‑5 (MQ=25)
                              gTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCtgtg                                                                                                                                <  1:3212166/143‑1 (MQ=25)
                              gTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCATGTGCAGGTcc                                                                                                                         >  2:6545129/1‑150 (MQ=255)
                                atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGGGGTGAACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGACCATTCGCGGGCTTTCCCGCACCCCCACGCATTGGACAGGGGTAGGACCAATTGCAGGTcccg                                                                                                                       >  2:4851575/1‑148 (MQ=255)
                                 tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                      <  2:534666/150‑1 (MQ=255)
                                 tagccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                      <  2:3049337/147‑1 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACCCCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCCGTCGCGGGCGTGCCCGCGCCCCCACGCCGTGCCCCGGGGTCGGCCCGTgggccgggccggggc                                                                                                                    >  2:4760061/1‑145 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                    <  2:3153240/150‑2 (MQ=21)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                    <  2:3144143/150‑1 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                    >  2:4528071/1‑150 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                    <  2:4213882/150‑1 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                    <  2:3816964/150‑1 (MQ=255)
                                    ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCg                                                                                                                   >  2:3385067/1‑150 (MQ=255)
                                     cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                                  <  2:3214463/150‑1 (MQ=255)
                                     cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                                  <  2:366555/150‑1 (MQ=255)
                                      aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAg                                                                                                                 >  2:2685153/1‑150 (MQ=255)
                                       gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                <  1:6614017/150‑1 (MQ=18)
                                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                              <  1:745682/150‑1 (MQ=18)
                                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                              <  2:6266396/150‑1 (MQ=18)
                                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                              <  2:2897510/150‑1 (MQ=255)
                                             gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCaa                                                                                                          <  1:1313551/150‑1 (MQ=255)
                                               ccAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGc                                                                                                        <  2:6382076/150‑1 (MQ=255)
                                                       ctgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATGGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                                <  2:4426778/150‑1 (MQ=255)
                                                         gctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                                          >  2:2307395/1‑106 (MQ=255)
                                                           tgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAATCGCGGGGGTGCCCCCGCCCCCACGCCATGCACACGGTTCGGACCTTTTGCAGGGCCGGTGCGAGCGGGCAAGCCCCCGCAAGGGc                                                                                            >  1:4961850/1‑150 (MQ=255)
                                                                aaCACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                  <  2:3623192/123‑1 (MQ=14)
                                                                aaCACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                  >  1:3623192/1‑123 (MQ=14)
                                                                  cacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCCTTGCACAGGGTTCGGACCGTGTGCCGGGCCGGTGCGAGCGGGCAAGCCCCCCCATGGGcccccccc                                                                                     >  2:2490643/1‑146 (MQ=37)
                                                                   acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCCATTCCCACGGTTCGGACCATGTGCAGGTCCGGTGCGAGCCGGCAAGCCCCCGCAAGGGcccccccca                                                                                    >  2:3764060/1‑145 (MQ=255)
                                                                                   aCGGCACACCGACGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGaa                                                                    <  2:3366593/150‑1 (MQ=255)
                                                                                                 tAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCTCCCACAGGCAATTCCCACGGTTGGGGCCGTGTGCAGGTCCAGTGGCTGCGGGCCACCCCCCGAATGGGCCCCCCCCAATTTTACCGGACTGAACAGTATAACCGccc                                                      >  2:921618/1‑150 (MQ=255)
                                                                                                    acgggggCCGGGTGGCGGCCCATTCGCGGGCTGGCCCGCTCCCCCCGGCTTTGCACCGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                                   <  2:1822868/147‑1 (MQ=38)
                                                                                                    acgggggCCGCGTTGCGGCCCTTTCGCGGGCTTGCCCGCTCCCACCGGCTTTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                                   <  2:1822875/147‑1 (MQ=38)
                                                                                                     cgggggCCGCGTTGGGGCCCTTTCGCGGGCTTGCCCGCTCCCACCGGATTTGCACTGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAg                                                  <  1:2238633/148‑1 (MQ=37)
                                                                                                     ccGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCCCAGGGTTCGGACCTTGTGCAGGTCCTGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACAATATCCCCGCCCCAAg                                                  >  2:2183062/1‑150 (MQ=255)
                                                                                                         ggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGTTCCTGTGCGAGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                              <  2:6429213/150‑1 (MQ=255)
                                                                                                                 ggcgcCCCATGCGCGGGCTGCCCCGCTCCCACAGCCATACCACAGGTTGCGGACCTTGTGCAGGCCCGGTCCGAGCGGCCACCCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                      <  1:6563770/149‑1 (MQ=255)
                                                                                                                 ggCGGCCCATTGGCGGGCTTGCCCGCTCCCACCGACATTGCACTGGTTTGGGACCTTTTGCCGTTCCTGCGGGTGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                      <  1:5400943/149‑1 (MQ=255)
                                                                                                                  gCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                     <  1:6257621/150‑1 (MQ=255)
                                                                                                                  gCGGCCACTCCGCGGGCTGCCCCCCCCCCCCCGGCTTTCCCCTGGGTGCGCCCCTTGTGCCGGCCGGGTGCGCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                     <  1:1821019/150‑1 (MQ=255)
                                                                                                                           tCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcacac                            >  2:6155404/1‑150 (MQ=255)
                                                                                                                               gggCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCt                        >  1:4118715/1‑150 (MQ=255)
                                                                                                                                gggttGCCCGCTCCCACAGGCTTTGCACAGGTTTCGGACGTTGTGCTGGTCCTGTGGCCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg                       <  1:3764060/147‑1 (MQ=255)
                                                                                                                                      cccGCTCCCACCGGCATTGCACGGGGTTCGGACCTTGTGCTGGTCCTGTGCCCGCGGCCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                 <  1:4871250/150‑1 (MQ=255)
                                                                                                                                      cccGCGCCCACCGGCATGGACCTGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                 <  1:1527709/150‑1 (MQ=255)
                                                                                                                                       ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                >  1:2341044/1‑150 (MQ=18)
                                                                                                                                       ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                >  1:2493877/1‑150 (MQ=18)
                                                                                                                                       ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                >  1:6747661/1‑150 (MQ=18)
                                                                                                                                       ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                >  1:4578625/1‑150 (MQ=18)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                >  2:4257150/1‑150 (MQ=255)
                                                                                                                                        cGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt               >  2:1226682/1‑150 (MQ=18)
                                                                                                                                        cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt               >  2:1793444/1‑150 (MQ=255)
                                                                                                                                         gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAgg                                    >  2:1509812/1‑128 (MQ=14)
                                                                                                                                         gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAgg                                    <  1:1509812/128‑1 (MQ=14)
                                                                                                                                         gCTCCCACAGGCATTGCACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt              <  2:2663207/150‑1 (MQ=255)
                                                                                                                                         gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt              <  2:5138269/150‑1 (MQ=255)
                                                                                                                                           tCCCACAGGCATTGCCCTGGGTGCGGACCTTGTGCCGGCCCTGTGCGCGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGa            <  1:2077894/150‑1 (MQ=255)
                                                                                                                                           tCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGa            >  1:6401612/1‑150 (MQ=255)
                                                                                                                                             ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt          >  1:941622/1‑150 (MQ=255)
                                                                                                                                               acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga        <  2:4300408/150‑1 (MQ=255)
                                                                                                                                                 cggCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga      <  2:3810961/149‑1 (MQ=255)
                                                                                                                                                 aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGc                                                        <  1:2684476/100‑1 (MQ=255)
                                                                                                                                                 aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga      <  2:4546436/150‑1 (MQ=255)
                                                                                                                                                   gCATTGCACTGGGGGCGGACCTTGTGCCGGCCCTGTGCTCGGGGCCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg    <  1:1808285/150‑1 (MQ=255)
                                                                                                                                                   gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGa                                                                         <  2:5782383/81‑1 (MQ=255)
                                                                                                                                                   gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGa                                                                         >  1:5782383/1‑81 (MQ=255)
                                                                                                                                                   gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg    >  2:2375552/1‑150 (MQ=255)
                                                                                                                                                   gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCACACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg    >  1:564510/1‑150 (MQ=255)
                                                                                                                                                     aTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGt  >  1:3926710/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGT  >  NC_002947/2068036‑2068334

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: