Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,056,057 T→A 5.5% intergenic (‑63/+23) PP_1825 ← / ← PP_1826 hypothetical protein/isochorismatase family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,056,0570TA 5.5% 1178.0 / 39.4 386intergenic (‑63/+23)PP_1825/PP_1826hypothetical protein/isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (168/187);  new base A (15/7);  total (192/194)
Fisher's exact test for biased strand distribution p-value = 7.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 390 total aligned reads displayed.

CGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCG  >  NC_002947/2055908‑2056206
                                                                                                                                                     |                                                                                                                                                     
      gccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCggg                                                                                                                                                 >  1:6400773/1‑150 (MQ=255)
      gccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCggg                                                                                                                                                 <  1:2423733/150‑1 (MQ=255)
       ccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGGCGGGc                                                                                                                                                >  1:5564958/1‑150 (MQ=255)
       ccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGc                                                                                                                                                <  1:5379880/150‑1 (MQ=255)
       ccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGc                                                                                                                                                <  2:5014669/150‑1 (MQ=255)
         ggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCgg                                                                                                                                              <  1:2595698/150‑1 (MQ=255)
           ccgccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                            <  2:6871816/150‑3 (MQ=255)
             gccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCt                                                                                                                                          >  2:130031/1‑150 (MQ=255)
               cAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGc                                                                                                                                        <  2:3359565/150‑1 (MQ=255)
                 ggcgggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCaa                                                                                                                                      <  2:2348685/150‑1 (MQ=255)
                    gggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGcc                                                                                                                                   >  2:3553313/1‑150 (MQ=255)
                        gcCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCCGCTGCCGGGCGGTCCGGCAAGccccgc                                                                                                                               >  1:5854347/1‑146 (MQ=255)
                         ccccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCa                                                                                                                              <  1:6582311/150‑1 (MQ=255)
                          cccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCCGGCGGTCCTGCAAGCCGCTCac                                                                                                                             >  1:1850080/1‑149 (MQ=255)
                             aGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATc                                                                                                                          <  2:2682105/150‑1 (MQ=255)
                               ttACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAg                                                                                                                        <  1:2346750/150‑1 (MQ=255)
                                tACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCCAGCCGCTCAGATCAgg                                                                                                                       >  1:4097728/1‑150 (MQ=255)
                                 aCGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGc                                                                                                                      >  2:4543289/1‑150 (MQ=255)
                                   gcgGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTg                                                                                                                    <  2:3988658/150‑1 (MQ=255)
                                   gcgGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTg                                                                                                                    >  2:1251040/1‑150 (MQ=255)
                                       ttCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCgg                                                                                                                >  2:6134007/1‑150 (MQ=255)
                                        tCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGt                                                                                                               <  2:6831789/150‑1 (MQ=255)
                                        tCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGt                                                                                                               <  2:1580061/150‑1 (MQ=255)
                                          gcgcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgg                                                                                                             <  2:3829585/150‑1 (MQ=255)
                                           cgcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggg                                                                                                            <  1:6236499/150‑1 (MQ=255)
                                            gcgGGCGCAGTCGATTGGCGTCTAAACTCCCTTTTAGCTTCATTCCGCTCTCATGGTTGGATATACGGAATTCGAGGTGCTGCTGCCGATCGGGTTCCGCCCGTCTGCAGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGATGTgggg                                                                                                           <  2:1670221/150‑1 (MQ=255)
                                            gcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgggg                                                                                                           <  2:2016159/150‑1 (MQ=255)
                                                gcgcAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCa                                                                                                       <  1:1340217/150‑1 (MQ=255)
                                                   cAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGc                                                                                                    <  1:2968762/150‑1 (MQ=255)
                                                     ggCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCaa                                                                                                  <  1:3672631/150‑1 (MQ=255)
                                                      gcgcTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAg                                                                                                 >  2:5964255/1‑150 (MQ=255)
                                                          ttGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                             <  1:6034439/150‑1 (MQ=255)
                                                          ttGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                             <  2:83414/150‑1 (MQ=255)
                                                          ttGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                             <  1:4052101/150‑1 (MQ=255)
                                                              cGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTg                                                                                         >  1:6706988/1‑150 (MQ=255)
                                                               gTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGt                                                                                        <  1:2479134/150‑1 (MQ=255)
                                                               gTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGt                                                                                        <  2:6018857/150‑1 (MQ=255)
                                                                tCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTc                                                                                       <  2:793780/150‑1 (MQ=255)
                                                                 cGAACCTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCa                                                                                      <  2:1139210/150‑1 (MQ=255)
                                                                 cGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCCGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCa                                                                                      <  1:4908887/150‑1 (MQ=255)
                                                                  gAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAg                                                                                     <  2:5759364/150‑1 (MQ=255)
                                                                   aaaCTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGc                                                                                    <  1:1803733/150‑1 (MQ=255)
                                                                    aaCTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGcc                                                                                   >  2:5526836/1‑150 (MQ=255)
                                                                       tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATc                                                                                <  1:6080201/150‑1 (MQ=255)
                                                                       tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATc                                                                                <  1:1251040/150‑1 (MQ=255)
                                                                       tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATc                                                                                >  1:1051360/1‑150 (MQ=255)
                                                                           tttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                            >  2:4368826/1‑150 (MQ=255)
                                                                           tttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                            >  2:1387062/1‑150 (MQ=255)
                                                                           tttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCCGCTGCCCGGCGGTCCTGCAAGCAGCACACATCAGGCTGGCGGGGGGGGGCACGCAAGCGAAGATGTAACCCATCACgg                                                                            >  1:3287067/1‑150 (MQ=255)
                                                                             tAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCCTCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGcc                                                                          <  1:1003019/150‑1 (MQ=255)
                                                                                gCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCGGCCGGGCGGGCCTGCAAGCCGCGCAGATCAGGCTGGCGGGGGGGGCCACGCAAGCGAAGTTGTCAGCCCTCACGGCCCtt                                                                       >  1:6185251/1‑150 (MQ=255)
                                                                                  ttCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTc                                                                     <  1:2015531/150‑1 (MQ=255)
                                                                                   tCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCg                                                                    <  1:2733587/150‑1 (MQ=255)
                                                                                     aTTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCa                                                                  <  1:5624168/150‑1 (MQ=255)
                                                                                        ccGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTg                                                               <  1:1089564/150‑1 (MQ=255)
                                                                                         cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTgg                                                              >  2:2774977/1‑150 (MQ=255)
                                                                                            tctcCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGa                                                           <  1:2208240/150‑1 (MQ=255)
                                                                                            tctcCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGa                                                           <  2:3168370/150‑1 (MQ=255)
                                                                                            tctcCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGa                                                           <  2:3011101/150‑1 (MQ=255)
                                                                                                 tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCGGCCGGGCGGTCCTGCAAGCCGCTCCGATCAGGCTGGCGGTGGGGGCGACGCAAGCGCAGTTGTCAGCCCTCACGGCCAATTAGGACAGGTGGATCTg                                                      >  1:3566037/1‑150 (MQ=255)
                                                                                                  cGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTGCGCCCCGCTGCCGGGGGGTCCGGCAAGCCCCTAAGATAAAGCTGGCGGTGGGGGGCACCCAAGCCAAGTTTTCAGCCAACACGGCCCATTCGAACAGGTTAATATgc                                                     >  1:310757/1‑150 (MQ=255)
                                                                                                  cGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgc                                                     <  1:96318/150‑1 (MQ=255)
                                                                                                   ggTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgcg                                                    >  2:5090457/1‑150 (MQ=255)
                                                                                                    gTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgcgc                                                   <  2:528142/150‑1 (MQ=255)
                                                                                                    gTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgcgc                                                   <  2:6728126/150‑1 (MQ=255)
                                                                                                         tatCTGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCg                                              >  2:1779970/1‑150 (MQ=255)
                                                                                                         tatCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCGGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCg                                              >  2:306585/1‑150 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGa                                            <  2:6637170/150‑1 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACAGGTGAATCTGCGCGGCCGGa                                            >  2:5181932/1‑150 (MQ=255)
                                                                                                             cGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAAt                                          <  2:2424808/150‑1 (MQ=255)
                                                                                                             cGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAATTAGCACTGGTGAATCTGCGCGGCCGGAAt                                          >  1:4638631/1‑150 (MQ=255)
                                                                                                              ggAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATg                                         >  2:253698/1‑150 (MQ=255)
                                                                                                                  ttCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCCAGCGAAGTTGTCAGCCCTCACGGCCCTTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGc                                     >  1:2818277/1‑150 (MQ=255)
                                                                                                                  ttCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGc                                     <  2:6225247/150‑1 (MQ=255)
                                                                                                                     gACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGt                                  >  2:6535482/1‑150 (MQ=255)
                                                                                                                     gACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGt                                  >  1:4231863/1‑150 (MQ=255)
                                                                                                                       ccTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTcc                                <  1:6830690/150‑1 (MQ=255)
                                                                                                                          gctgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCCAGCCGCTCAGATCAGGCGGGCGGTGGGGGCCCCGCCAGCGAAGTTGTCAGCCATAACCGCCCTTTCGCACAGGTGCAACAGCGCGGCCGGAAAGACGCCGGCCTTg                             >  1:3769207/1‑150 (MQ=255)
                                                                                                                          gctgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCCGGCCTGCAACCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACAGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTg                             >  1:1935753/1‑150 (MQ=255)
                                                                                                                            ggctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGaa                           <  2:3728293/149‑1 (MQ=255)
                                                                                                                             gctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAAt                          <  2:4724489/150‑1 (MQ=255)
                                                                                                                             gctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAAt                          <  2:5925175/150‑1 (MQ=255)
                                                                                                                             gctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCGGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAAACTGCGCGGCCGGAATGACGCCGTCCGTGAAt                          >  1:5280135/1‑150 (MQ=255)
                                                                                                                                 cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTTTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                      >  2:2184627/1‑150 (MQ=255)
                                                                                                                                 cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                      <  1:4436486/150‑1 (MQ=255)
                                                                                                                                  ccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAg                     <  1:1464345/150‑1 (MQ=255)
                                                                                                                                    aTCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGt                   <  2:2074613/150‑1 (MQ=255)
                                                                                                                                       gggTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCAc                >  2:5580739/1‑150 (MQ=255)
                                                                                                                                       caggaccGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCAc                >  1:325338/6‑150 (MQ=255)
                                                                                                                                        aggaccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACg               >  2:4878588/5‑150 (MQ=255)
                                                                                                                                         gTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg              <  2:4097728/150‑1 (MQ=255)
                                                                                                                                         gTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg              <  1:837301/150‑1 (MQ=255)
                                                                                                                                         gTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg              <  2:107667/150‑1 (MQ=255)
                                                                                                                                          ttCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAGGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACggg             <  2:5400110/150‑1 (MQ=255)
                                                                                                                                          gaccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACggg             >  2:4917320/3‑150 (MQ=255)
                                                                                                                                           accGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGt            >  2:1915323/2‑150 (MQ=255)
                                                                                                                                            ccGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTg           <  1:4661921/150‑1 (MQ=255)
                                                                                                                                               cccGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGcc        >  2:443886/1‑150 (MQ=255)
                                                                                                                                                 cGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGa      >  1:1711987/1‑150 (MQ=255)
                                                                                                                                                    cTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAcgc   <  1:3374449/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCG  >  NC_002947/2055908‑2056206

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: