Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 537,589 A→C 5.5% intergenic (+116/‑37) PP_mr07 → / → rpoB P26 ncRNA/DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947537,5890AC 5.5% 660.1 / 18.4 191intergenic (+116/‑37)PP_mr07/rpoBP26 ncRNA/DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (84/95);  new base C (6/5);  total (91/100)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

AAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTATCGCGAATT  >  NC_002947/537441‑537738
                                                                                                                                                    |                                                                                                                                                     
acaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                                      <  2:1819754/148‑1 (MQ=255)
aaaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                                      <  1:3153282/150‑1 (MQ=255)
   ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGcc                                                                                                                                                   <  1:23122/150‑1 (MQ=255)
    gccgATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGCCa                                                                                                                                                  <  1:618536/147‑1 (MQ=255)
    gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCa                                                                                                                                                  >  1:2456050/1‑150 (MQ=255)
      tGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa                                                                                                                                                >  1:3415540/1‑150 (MQ=255)
      tGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa                                                                                                                                                <  2:896667/150‑1 (MQ=255)
      tGATGGCTGGTGGCTTATGCCACCGGCCTTTGTCCGGTATTGGTGGCCGGTCAGGCCGGTCCCGAAAACGAGCAGCAAGACACCCCGGTGGGAGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGACGCCcaa                                                                                                                                                >  1:1750842/1‑150 (MQ=255)
           gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc                                                                                                                                           >  1:4251825/1‑150 (MQ=255)
           gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc                                                                                                                                           <  2:97030/150‑1 (MQ=255)
            cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa                                                                                                                                          <  1:90773/150‑1 (MQ=255)
            cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa                                                                                                                                          <  1:2031017/150‑1 (MQ=255)
             tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg                                                                                                                                         <  2:3603492/150‑1 (MQ=255)
             tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg                                                                                                                                         >  2:2827148/1‑150 (MQ=255)
                 ggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGa                                                                                                                                     <  2:3475218/150‑1 (MQ=255)
                  gCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAc                                                                                                                                    >  1:399986/1‑150 (MQ=255)
                   cTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAcc                                                                                                                                   <  2:3826381/150‑1 (MQ=255)
                   cTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAcc                                                                                                                                   >  1:2246322/1‑150 (MQ=255)
                        gCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCt                                                                                                                              >  1:384619/1‑150 (MQ=255)
                        gCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCt                                                                                                                              <  1:2811587/150‑1 (MQ=255)
                        gCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCt                                                                                                                              <  2:661350/150‑1 (MQ=255)
                         ccACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTg                                                                                                                             <  2:2356788/150‑1 (MQ=255)
                         ccACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTg                                                                                                                             >  2:517020/1‑150 (MQ=255)
                          cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                            <  1:13715/150‑1 (MQ=255)
                          cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                            >  1:4177260/1‑150 (MQ=255)
                          cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                            <  1:4157932/150‑1 (MQ=255)
                          cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                            <  1:3459207/150‑1 (MQ=255)
                          cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                            <  2:171235/150‑1 (MQ=255)
                                  tttttCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCt                                                                                                                    <  1:1327938/150‑1 (MQ=255)
                                     ttCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGAt                                                                                                                 <  2:1401863/150‑1 (MQ=255)
                                         gTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa                                                                                                                                          <  2:2999081/121‑1 (MQ=255)
                                         gTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa                                                                                                                                          >  1:2999081/1‑121 (MQ=255)
                                         gTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCt                                                                                                             >  2:1833171/1‑150 (MQ=255)
                                         gTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCt                                                                                                             >  2:3653154/1‑150 (MQ=255)
                                          ttATTGGTGGGCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACGGATGTGGGTGGCGCAATCGATGGGGTTTGCTCGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACTAGCAGGTGACCTAGGTGGGGAACGCTGATGGCtt                                                                                                            >  2:1154593/1‑150 (MQ=255)
                                            aTTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCAGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTAc                                                                                                          >  2:3149186/1‑150 (MQ=255)
                                               ggtggCCGGTCAGGCCGGTCCCGATAACGAGCTTCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGACTGTAGAAGCCTAACAAGCATGTGACCAAGCTGTGGAACGCTTATGGATTACTCa                                                                                                       <  2:2258479/150‑1 (MQ=255)
                                               ggtggCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCa                                                                                                       >  1:3922430/1‑150 (MQ=255)
                                                 tggCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATa                                                                                                     <  1:1797489/150‑1 (MQ=255)
                                                     cGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACt                                                                                                 >  1:2578099/1‑150 (MQ=255)
                                                      ggTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTg                                                                                                <  1:748464/150‑1 (MQ=255)
                                                      ggTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTg                                                                                                <  1:4264614/150‑1 (MQ=255)
                                                        tCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTgag                                                                                              <  1:3546505/150‑1 (MQ=255)
                                                               ggTCCCGATATCGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAAc                                                                                       <  2:2936444/150‑1 (MQ=255)
                                                               ggTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAAc                                                                                       <  1:119089/150‑1 (MQ=255)
                                                               ggTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAAc                                                                                       <  1:4442919/150‑1 (MQ=255)
                                                                   ccGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTAt                                                                                   <  2:1978710/150‑1 (MQ=255)
                                                                       tAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGc                                                                               >  2:2286428/1‑150 (MQ=255)
                                                                          cGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAg                                                                            <  2:2163829/150‑1 (MQ=255)
                                                                          cGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAg                                                                            <  1:3130744/150‑1 (MQ=255)
                                                                           gAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGGAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAgg                                                                           >  1:3469550/1‑150 (MQ=255)
                                                                               tGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACtt                                                                       >  2:648126/1‑150 (MQ=255)
                                                                               tGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACtt                                                                       <  1:3180955/150‑1 (MQ=255)
                                                                               ggCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACtt                                                                       <  1:2352410/149‑1 (MQ=255)
                                                                                 cAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTa                                                                     >  1:4259038/1‑150 (MQ=255)
                                                                                 cAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTa                                                                     >  2:764965/1‑150 (MQ=255)
                                                                                       acCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGt                                                               >  1:1391018/1‑150 (MQ=255)
                                                                                       acCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGt                                                               >  1:1561633/1‑150 (MQ=255)
                                                                                         ccATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTg                                                             >  1:1029314/1‑150 (MQ=255)
                                                                                         ccATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTg                                                             >  1:4304543/1‑150 (MQ=255)
                                                                                         ccATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTg                                                             >  2:1766090/1‑150 (MQ=255)
                                                                                         ccATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTg                                                             >  2:1474739/1‑150 (MQ=255)
                                                                                          cATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACTAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTTCTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTATCAAGTTGc                                                            <  1:3806399/150‑1 (MQ=255)
                                                                                           aTGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGcc                                                           <  1:1359254/150‑1 (MQ=255)
                                                                                           aTGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGcc                                                           <  2:3031372/150‑1 (MQ=255)
                                                                                            tgtgGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCg                                                          <  1:3676250/150‑1 (MQ=255)
                                                                                                  tggCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTc                                                    >  2:1919706/1‑150 (MQ=255)
                                                                                                      gcAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATgg                                                >  1:497852/1‑150 (MQ=255)
                                                                                                      gcAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATgg                                                >  1:1306793/1‑150 (MQ=255)
                                                                                                          aCCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGAtgt                                            <  2:2974927/150‑1 (MQ=255)
                                                                                                          aCCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGAtgt                                            <  1:2430761/150‑1 (MQ=255)
                                                                                                           ccAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGAtgtg                                           >  2:4080616/1‑150 (MQ=255)
                                                                                                            cAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGc                                          <  2:1669616/150‑1 (MQ=255)
                                                                                                            cAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGc                                          >  1:2510963/1‑150 (MQ=255)
                                                                                                               ggggTTTGCACGATGGTCGGGCTGCGCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCtt                                       <  2:1750842/150‑1 (MQ=255)
                                                                                                                 ggTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAc                                     <  1:1445037/150‑1 (MQ=255)
                                                                                                                     tGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATacac                                                                                                  <  1:3104154/85‑1 (MQ=255)
                                                                                                                     tGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATacac                                                                                                  >  2:3104154/1‑85 (MQ=255)
                                                                                                                      gCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAcctcct                                >  1:790692/1‑150 (MQ=255)
                                                                                                                       cACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTg                               <  2:2363309/150‑1 (MQ=255)
                                                                                                                       cACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTg                               <  1:3301438/150‑1 (MQ=255)
                                                                                                                        aCGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTgg                              <  1:41085/150‑1 (MQ=255)
                                                                                                                           aTTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCa                           >  2:4259296/1‑150 (MQ=255)
                                                                                                                                 tGGCTTCTCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGc                     >  2:536129/1‑150 (MQ=255)
                                                                                                                                  ggCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACTCTGAGAAATAACGTATCCGCTAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCt                    <  2:2145114/150‑1 (MQ=255)
                                                                                                                                  ggCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCt                    >  2:2317617/1‑150 (MQ=255)
                                                                                                                                   gcttctCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAcctcc                                 <  1:1916857/136‑1 (MQ=255)
                                                                                                                                   gcttctCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAcctcc                                 >  2:1916857/1‑136 (MQ=255)
                                                                                                                                   gcttctCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTg                   >  1:1807868/1‑150 (MQ=255)
                                                                                                                                   gctgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTg                   >  2:4426407/1‑150 (MQ=255)
                                                                                                                                    cttctCCCGGCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTgg                  >  2:2481164/4‑150 (MQ=255)
                                                                                                                                    ctgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTgg                  <  1:274528/150‑1 (MQ=255)
                                                                                                                                     tgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGa                 <  1:2111456/150‑1 (MQ=255)
                                                                                                                                        tCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTc              <  2:4251825/150‑1 (MQ=255)
                                                                                                                                        tCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTc              >  1:2898257/1‑150 (MQ=255)
                                                                                                                                         cccGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCg             <  2:3998714/150‑1 (MQ=255)
                                                                                                                                               gggAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTATcgc       <  2:4177260/150‑1 (MQ=255)
                                                                                                                                                 gagaAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACttt                                                                      <  2:1447803/85‑1 (MQ=255)
                                                                                                                                                 gagaAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTATCGCGa     >  2:3346801/1‑150 (MQ=255)
                                                                                                                                                    aaGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTATCGCGAAtt  >  1:675140/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCGTATCGCGAATT  >  NC_002947/537441‑537738

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