Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,184 C→A 13.3% intergenic (+49/‑158) 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,1840CA13.3% 600.3 / 78.3 227intergenic (+49/‑158)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (86/111);  new base A (10/20);  total (96/131)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-02

Only 100 of 228 total aligned reads displayed.

CCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGG  >  NC_002947/2068036‑2068332
                                                                                                                                                    |                                                                                                                                                    
  aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                                   >  1:3162640/1‑150 (MQ=255)
  aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                                   <  1:2293793/150‑1 (MQ=255)
  aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                                   <  1:2527507/150‑1 (MQ=255)
       gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                              >  1:3712841/1‑150 (MQ=255)
        aaaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcacaa                                                                                                                                             <  2:5573485/150‑2 (MQ=255)
          aTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAggg                                                                                                                                           <  2:1003059/150‑1 (MQ=255)
              tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGc                                                                                                                                                 <  2:608730/140‑1 (MQ=255)
              tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCg                                                                                                                                       <  2:1250571/150‑1 (MQ=255)
                aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                                     <  2:3553397/150‑1 (MQ=255)
                  gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCATTGCACAGGGTTCGGGcc                                                                                                                                   >  2:3163109/1‑150 (MQ=255)
                    agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                 <  2:4311065/150‑1 (MQ=255)
                    agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                 <  1:6062977/150‑1 (MQ=255)
                     gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAACCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtg                                                                                                                                <  2:5530214/150‑1 (MQ=255)
                     gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtg                                                                                                                                <  1:4409228/150‑1 (MQ=255)
                     gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtg                                                                                                                                <  1:3459053/150‑1 (MQ=255)
                        ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                             <  2:3917909/150‑1 (MQ=255)
                         gAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCa                                                                                                                            <  2:5065072/150‑1 (MQ=255)
                         gAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGCTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACACGCATTGCACAGGGTTCGGACCTTGTGCa                                                                                                                            <  2:3948677/150‑1 (MQ=255)
                          aaGCGTATATCCAGACCGCGACCACGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAACCCGGGGCCGCGTTGCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAg                                                                                                                           <  2:839087/150‑1 (MQ=255)
                          aaGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAg                                                                                                                           <  1:4952185/150‑1 (MQ=255)
                           aGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCaat                                                                                                                          <  1:6323224/150‑3 (MQ=255)
                             cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                                        <  2:5869970/150‑1 (MQ=255)
                              gTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACCGGGGTCGGACCTTGTGCAGGTcc                                                                                                                       >  2:2585375/1‑150 (MQ=255)
                               tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                      <  1:1495331/150‑1 (MQ=255)
                                atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgcctg                                                                                                                     <  2:5409232/150‑8 (MQ=21)
                                 tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                    <  1:7278009/150‑1 (MQ=255)
                                 tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                    <  1:1977864/150‑1 (MQ=255)
                                  atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                              >  2:2152952/1‑123 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                  <  1:5007164/150‑1 (MQ=255)
                                   tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                  <  1:6107652/150‑1 (MQ=255)
                                    ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCGTGCCCGCGCCCACAGGCATTGCCCACGGTTGGGTCCGTGTGCAGGTGCTGTGCg                                                                                                                 >  1:5423409/1‑150 (MQ=255)
                                     cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                <  1:2648510/150‑1 (MQ=18)
                                     cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGTGCCTGTGGGa                                                                                                                <  1:5473568/150‑1 (MQ=255)
                                      aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCCTTGCACAGGGTTCGGACCGTGTGCAGGGCCGGTGCggg                                                                                                               >  2:372082/1‑148 (MQ=255)
                                      aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCACCCACAGGCCATTCACACGGTTCGGCCCGTGTGCAGGTCCTGTGCGAg                                                                                                               >  2:3746608/1‑150 (MQ=255)
                                      aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCTCCCACAGGCCTTGCACAGGGGTCGGACCTTGTGCAGGTCCGGGGCGAg                                                                                                               >  2:3821911/1‑150 (MQ=255)
                                       gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                              <  2:6827179/150‑1 (MQ=255)
                                       gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCCCGGGCGTGCCCGCGCCCACACGCAGTGCACCGGGGTCGGTCCAAGTGAAGGTCCTGTGCGGGc                                                                                                              >  1:6732281/1‑150 (MQ=255)
                                        aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCg                                                                                                             <  2:1351784/150‑1 (MQ=18)
                                        aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCg                                                                                                             <  2:3914131/150‑1 (MQ=18)
                                         ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                            <  1:1738148/150‑1 (MQ=18)
                                           gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGc                                                                                                          <  2:5468496/150‑1 (MQ=255)
                                             gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCaa                                                                                                        <  1:4635098/150‑1 (MQ=255)
                                                   gagCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGTCCGAAccct                                                                                                                                 >  1:2130928/1‑117 (MQ=255)
                                                     gCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCTTGCCcgctc                                                                                                <  1:2272570/150‑3 (MQ=255)
                                                     gCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGa                                                                                                <  2:7373582/150‑1 (MQ=255)
                                                        tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAGTCGCGGGCTTGCCCCCTCCCACAGGCATTGCACAGGGTTCGGGCCATGTGCAGGGCCTGGGCGAGCCGGCAAGCCcgcggctg                                                                                             >  2:4087430/1‑146 (MQ=255)
                                                             gATAACACCGCCGGCGTGGTGCCCGGCACACCGAAGTAAACGGGGGCCGCGTTGCGGCCCATTCGGGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCGTGTGCGAGCGGGCAAGCCCGCGAATGGGccc                                                                                        <  2:4431915/150‑1 (MQ=38)
                                                             gATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGCATGGGccc                                                                                        >  2:2946535/1‑150 (MQ=255)
                                                              aTAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAACCCGGGGCCGCGTTGCGGCCCATGCGCGGGCTTGCCCGCTCCCACAGGCATGGCACAGGGTTCGGACCTGGTGCAGGTCCTGTGGGAGCGGGCAAGCCCGCGAATGGGcccc                                                                                       <  2:5944817/150‑1 (MQ=38)
                                                              aTAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCTACCACAGGCAATTCACAGGGTTAGGACCATGTGCAGGTCCTGTGCGAGCCGGCAAGCCCCCGAATGGGcccc                                                                                       >  2:1999748/1‑150 (MQ=255)
                                                                 acacCGCGGGCGTGGTGCACGGCACACCGAAGTACACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGca                                                                                    <  2:5297633/150‑1 (MQ=255)
                                                                     cgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCCTTGCACAGGGTTCGGACCTTGTGCCGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACCAt                                                                                >  2:7419802/1‑150 (MQ=255)
                                                                     cgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCGCCCACAGGCCGTGCACACGGTGCGGGCCGTGTGCCGGGCCTGGGCGAGCGGGCCAGCCCGCGACGGGGCCCCCCACCAt                                                                                >  2:657579/1‑150 (MQ=255)
                                                                              ggtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGa                                                                       <  2:5754671/150‑1 (MQ=255)
                                                                                  cACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGGTCGGACCTTGTGCAGGGCCTGTGCGAGCGGGCCAGCCCGCCAAGGGGCCCCGCCCAATTTTCCcggacgga                                                                   >  2:3038296/1‑150 (MQ=38)
                                                                                  cACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAGTGCACCGGGTTCGGACCGTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGCATGGGCCCCGCACAATTTTCCCGGACTGa                                                                   >  2:4521141/1‑150 (MQ=255)
                                                                                                    aCCGGGGCGGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCTTTGAACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                                 <  2:363774/150‑1 (MQ=38)
                                                                                                     ccGGGGCCGCGTTGCGGCCCATTCGCGGGCGTTCCCCCGCCCACACGCATTTCACCGGGGTGGGGCCTTGTGCCGGGCCTGTGGCCGCGGGCCAGCCCGCGCATGGGGCCCGCCCAAATTTGCCGGACTGAACAATATACCCGCCCCCAg                                                >  2:6916591/1‑150 (MQ=255)
                                                                                                         ggCCGGGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACCGGGTTCGGACCTGGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                            <  1:3053018/150‑1 (MQ=255)
                                                                                                         ggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                            <  2:6795335/150‑1 (MQ=255)
                                                                                                           ccGCGTTGCGGCCCATTCGCGGGCTTGCCCGCCCCCCCCGGCTTTGCACCGGGTTGGGACCTTGTGCTGGTCCTGTGCGCGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGt                                          <  1:1685226/150‑1 (MQ=255)
                                                                                                              cgTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAATGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCCCGAAAGGGCCCCGCACAATTTTCCCGGACTGAACAATATCCCCGCCCCAAGGTCAGTTAt                                       >  2:1752157/1‑150 (MQ=255)
                                                                                                               gTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGCACCTTGTGCTGGTCCTGTGCGAGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                      <  2:1562365/150‑1 (MQ=255)
                                                                                                                 ggCGGCCCTTTGGCGGGCTTGCCCGCTCCCCCTGGCTTTGCACTGGGTTGGGACCTTGTGCAGGTCCTGTGCGCGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                    <  2:1533743/149‑1 (MQ=255)
                                                                                                                  gCGGCCCATTGGGGGGCTGGCCCGCTCCCACCGGCATTGCACTGGGTTCGGACCTTGTGCCGGTCGGGTGGGCGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                   <  2:41645/150‑1 (MQ=255)
                                                                                                                  gCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                   <  1:7384531/150‑1 (MQ=255)
                                                                                                                   cGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCAATGCACACGGGTGGGACCGTGTTCCCGTCCGGTGCGAGCCGGCCACCCCCCGCAAGGGCCCCGCACAAATTTACCCGACAGAACAGTATACCCCCCCCCAGGGCCGGTATTAAgg                                  >  2:721580/1‑150 (MQ=255)
                                                                                                                    ggCCCTTTCGCGGGCTGGCCCGCTCCCAAAGGCTTTGCACAGGGTTCGGACCTTGTGATGGTCCTGTGCTAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                                 <  2:2081339/150‑1 (MQ=38)
                                                                                                                    ggCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                                 <  2:2383599/150‑1 (MQ=255)
                                                                                                                       ccATTCGCGGGCTTGCCCGCTCCCACCGGAATTGCACAGGGTTGGGACCTTGTGCTGGTCCTGTGCGAGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                              <  1:4202730/150‑1 (MQ=38)
                                                                                                                         aTTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAAATTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcac                            >  1:6160433/1‑150 (MQ=255)
                                                                                                                            cgcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcacaca                         <  2:6732281/150‑1 (MQ=255)
                                                                                                                             gcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCacacac                        <  2:3712841/150‑1 (MQ=255)
                                                                                                                             gcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCacacac                        >  2:6193923/1‑150 (MQ=255)
                                                                                                                              cgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACAcc                       >  2:4845967/1‑150 (MQ=255)
                                                                                                                                 gattGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                    <  1:7422042/148‑1 (MQ=255)
                                                                                                                                   ttGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg                  >  1:133446/1‑150 (MQ=18)
                                                                                                                                     gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGtt                >  2:4768418/1‑150 (MQ=255)
                                                                                                                                      cccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc               >  2:6358323/1‑150 (MQ=255)
                                                                                                                                      cccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc               >  2:1397728/1‑150 (MQ=255)
                                                                                                                                      cccGCCCCCCCCGGCTTTGCACTGGGGTCGGACGTGGGGCCGGCCGGGGGCGCGGGGGCACCCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc               <  1:5686211/150‑1 (MQ=255)
                                                                                                                                       ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg              >  2:3153837/1‑150 (MQ=18)
                                                                                                                                       ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGCGTTCg              >  1:1567208/1‑150 (MQ=18)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc               <  2:1618958/149‑1 (MQ=255)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg              >  1:5254615/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg              >  1:2234004/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg              >  1:3977613/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg              >  2:6641656/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg              >  2:6254589/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTCCg              >  1:7073715/1‑150 (MQ=255)
                                                                                                                                        cGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGaa                                                                  >  2:3604823/1‑97 (MQ=255)
                                                                                                                                         gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt            >  1:2790389/1‑150 (MQ=21)
                                                                                                                                         gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt            >  1:6865365/1‑150 (MQ=255)
                                                                                                                                             ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt        >  2:2623747/1‑150 (MQ=255)
                                                                                                                                             ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt        >  2:331031/1‑150 (MQ=255)
                                                                                                                                             ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt        >  1:1104957/1‑150 (MQ=255)
                                                                                                                                                 cggCATTGCACCGGGTTCGGACCTTTTGCAGTCCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga    <  2:3736445/149‑1 (MQ=255)
                                                                                                                                                 cggCATTCCACTGGGTCGGGACGGGGGGCTGGCCCTGGGCGCGCGGCCAACCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga    <  1:3821911/149‑1 (MQ=255)
                                                                                                                                                   gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg  >  2:4391977/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGG  >  NC_002947/2068036‑2068332

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: