Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,659,822 A→G intergenic (+62/+37) alpR → / ← PA14_RS21395 transcriptional regulator AlpR/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,659,8220AG97.8% 1329.3 / 10.1 372intergenic (+62/+37)alpR/PA14_RS21395transcriptional regulator AlpR/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base G (119/244);  total (126/246)
Fisher's exact test for biased strand distribution p-value = 2.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 391 total aligned reads displayed.

ACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACG  >  NC_008463/4659675‑4659970
                                                                                                                                                   |                                                                                                                                                    
   cGGACGAGCGCTACGAGCCCCAGGCCGCCGCCGAGCACATCCACGGGCTCGGCCGGGTGTTCTGGGACTCGGCGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAg                                                                                                                                                <  2:8147580/151‑1 (MQ=255)
   cGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCCAAGAAAAAg                                                                                                                                                >  2:8726685/1‑151 (MQ=255)
    ggACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGACGAAAAATCCCGGGCAAAGAAAAAgg                                                                                                                                               >  2:4943597/1‑151 (MQ=255)
      aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                             <  2:5539621/151‑2 (MQ=255)
      aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                             <  1:2253972/151‑2 (MQ=255)
      aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                             <  1:3760535/151‑2 (MQ=255)
       cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                            >  2:9138798/1‑149 (MQ=255)
       cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                            <  1:1635801/151‑3 (MQ=255)
       cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                            <  2:8974092/151‑3 (MQ=255)
        gAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcct                                                                                                                                           <  1:4612805/151‑4 (MQ=255)
         aGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcctg                                                                                                                                          <  1:10316619/151‑5 (MQ=255)
         aGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcctg                                                                                                                                          <  1:6068326/151‑5 (MQ=255)
           cgcTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCa                                                                                                                                       <  2:1300768/151‑1 (MQ=255)
              tACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtg                                                                                                                                    <  2:5947667/151‑1 (MQ=255)
                 gAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgctg                                                                                                                                 <  1:6187007/151‑1 (MQ=255)
                  aGCCCCAGGAAGCCGCCGAGCACCTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAATAGGGCCTGCATtgctgc                                                                                                                                >  2:6512735/1‑151 (MQ=255)
                    ccccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAg                                                                                                                              >  1:1287952/1‑151 (MQ=255)
                     cccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAgg                                                                                                                             <  2:4337579/151‑1 (MQ=255)
                     cccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAgg                                                                                                                             >  1:8253458/1‑151 (MQ=255)
                      ccAGGAGGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGc                                                                                                                            >  1:8056334/1‑151 (MQ=255)
                      ccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGc                                                                                                                            >  2:8790294/1‑151 (MQ=255)
                      ccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGc                                                                                                                            <  1:6021359/151‑1 (MQ=255)
                         ggAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCt                                                                                                                         <  1:7257636/151‑1 (MQ=255)
                             gccgccGAGCACAACAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGACTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCg                                                                                                                     <  2:6400944/151‑1 (MQ=255)
                              ccgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCgg                                                                                                                    >  1:8904725/1‑151 (MQ=255)
                                gccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTa                                                                                                                  <  1:10634132/151‑1 (MQ=255)
                                  cgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATg                                                                                                                <  2:4745137/151‑1 (MQ=255)
                                  cgAGCACATCAACGAGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATg                                                                                                                >  1:8960628/1‑151 (MQ=255)
                                   gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATgg                                                                                                               >  2:10914681/1‑151 (MQ=255)
                                        caTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCAc                                                                                                          <  2:207979/151‑1 (MQ=255)
                                         aTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACt                                                                                                         >  2:7694653/1‑151 (MQ=255)
                                         aTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACt                                                                                                         >  1:8781301/1‑151 (MQ=255)
                                          tCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTc                                                                                                        <  2:5351603/151‑1 (MQ=255)
                                              cGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagga                                                                                                    >  1:4966491/1‑151 (MQ=255)
                                              cGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagga                                                                                                    <  2:4114633/151‑1 (MQ=255)
                                               gTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCaggag                                                                                                   <  2:3515789/151‑1 (MQ=255)
                                                 gATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAggagga                                                                                                 >  2:4231361/1‑151 (MQ=255)
                                                  aTCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAACAATCCCGGGCAAAGAAATAGGGCCTGCATTGCTGCAGGCCCTTGCGGTATGGCGCACTCAggaggag                                                                                                >  2:11122752/1‑150 (MQ=255)
                                                  aTCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGAt                                                                                                >  1:5759412/1‑151 (MQ=255)
                                                    cGGCCGGGTGTTCTGCTACTCGGTGCTGGTCTGACCCCTGCCCGCCCCGCCTGACGACGGGGCCGCGATCCGGCGGAACACCAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTc                                                                                              <  2:10364421/151‑1 (MQ=255)
                                                           gtgtTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTc                                                                                       <  2:1464148/151‑1 (MQ=255)
                                                            tgttgtGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACTACGGGGCCGCGTTCCGGCGTAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTcc                                                                                      <  1:8852266/151‑1 (MQ=255)
                                                                cTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGAcc                                                                                  <  1:2940908/151‑1 (MQ=255)
                                                                  ggTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGc                                                                                >  1:7679599/1‑151 (MQ=255)
                                                                  ggTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGc                                                                                >  1:687694/1‑151 (MQ=255)
                                                                    tACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTc                                                                              >  2:619285/1‑151 (MQ=255)
                                                                        cGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                          <  1:2674594/151‑1 (MQ=255)
                                                                        cGGTGCGGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                          <  2:6227341/151‑1 (MQ=255)
                                                                         ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                         <  1:1096150/151‑1 (MQ=255)
                                                                         ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                         <  2:4655472/151‑1 (MQ=255)
                                                                          gTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCg                                                                        <  2:6768190/151‑1 (MQ=255)
                                                                           tGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGt                                                                       <  1:5429023/151‑1 (MQ=255)
                                                                             cTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     >  1:3865017/1‑151 (MQ=255)
                                                                             cTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     >  1:4249467/1‑151 (MQ=255)
                                                                              tGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGc                                                                    <  1:3437161/151‑1 (MQ=255)
                                                                               ggTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                   <  2:36229/151‑1 (MQ=255)
                                                                                   tGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGt                                                               >  1:2486620/1‑151 (MQ=255)
                                                                                    gACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                              <  2:5822316/151‑1 (MQ=255)
                                                                                        ccTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActct                                                          >  1:3105104/1‑151 (MQ=255)
                                                                                         ctccccgccccgcCTGACGCCGGGGCCGCGTTCCGGCGGAAGACCAATCCCGGGCAAAGAACAAGGGCCTGCCTTGCTGCAGGCCCTATCGGTAAGGCGCCCTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                         <  2:2839650/148‑1 (MQ=25)
                                                                                         cgtCCCCCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAATATCCCGGGCAAAGACAAAGCGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                         <  1:5670882/149‑1 (MQ=25)
                                                                                         cTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                         <  2:2629234/151‑1 (MQ=255)
                                                                                          ttCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                        >  1:7100367/1‑151 (MQ=255)
                                                                                            cccgccccgcCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                      >  2:668373/1‑151 (MQ=255)
                                                                                             ccgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                     <  2:9374643/151‑1 (MQ=255)
                                                                                               gccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                   <  1:10141455/151‑1 (MQ=255)
                                                                                               gccccgacTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                   <  1:9430951/151‑1 (MQ=25)
                                                                                                ccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                  <  2:3958924/151‑1 (MQ=255)
                                                                                                 cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCCTTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTCCTCTATCCAGCTg                                                 >  2:946472/1‑151 (MQ=37)
                                                                                                 cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 <  1:8593845/151‑1 (MQ=255)
                                                                                                    gccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                              <  1:5186554/151‑1 (MQ=255)
                                                                                                    gccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                              <  1:7297097/151‑1 (MQ=255)
                                                                                                     ccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                             >  1:1981945/1‑151 (MQ=255)
                                                                                                       tGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                           <  1:8385400/151‑1 (MQ=255)
                                                                                                       tGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                           >  2:5013980/1‑151 (MQ=255)
                                                                                                         acgacgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAATGTAAAAGGGCCTGCGTGGCTGCAGGCGCTTGCGGGTTGGCGCACTCTGGAGGTTTCGAACCTCCGACCGCTGGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTCTGa                                         >  2:2293296/1‑151 (MQ=18)
                                                                                                         acgacgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                         >  2:2392870/1‑151 (MQ=255)
                                                                                                                ggCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCGATGAGTGCgtg                                  >  2:161491/1‑151 (MQ=255)
                                                                                                                   cgcgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCg                               <  2:10603547/151‑1 (MQ=255)
                                                                                                                     cgTTCCGGCGGAAGAAAAATCCCTGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCggg                             >  2:10048435/1‑150 (MQ=255)
                                                                                                                     cgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGt                             <  1:10072370/151‑1 (MQ=255)
                                                                                                                      ggtCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTg                            <  2:7163101/149‑1 (MQ=255)
                                                                                                                            gcggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttttt                      <  2:8518585/151‑1 (MQ=255)
                                                                                                                              ggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTcc                    >  2:9784636/1‑151 (MQ=255)
                                                                                                                                aagaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCgg                  >  1:10028046/4‑151 (MQ=255)
                                                                                                                                    caaaTCCCGGGCACAGAACCAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATACGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCAc              <  2:3343499/150‑1 (MQ=21)
                                                                                                                                    aaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCAc              <  1:1471707/151‑1 (MQ=255)
                                                                                                                                      aaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCt            <  1:9152688/151‑1 (MQ=255)
                                                                                                                                      aaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCAGGCCACCt            >  1:10796474/1‑151 (MQ=255)
                                                                                                                                        tCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGc          <  1:4070269/151‑1 (MQ=255)
                                                                                                                                        acccGCGCAAAGCAAAATGGACTGCATTTCTGCAGGCCCTTTCGGTATGTCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGATGAGCAATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGc          <  2:5158338/150‑1 (MQ=18)
                                                                                                                                         cccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGcc         >  1:6740818/1‑151 (MQ=255)
                                                                                                                                          ccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCg        <  1:2161145/151‑1 (MQ=255)
                                                                                                                                           cGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCgg       <  1:9696072/151‑1 (MQ=255)
                                                                                                                                            gggCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGGGCGTGTTACGGCATCCTGcctga      >  2:9903193/1‑148 (MQ=18)
                                                                                                                                             ggCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTa     <  1:7843284/151‑1 (MQ=255)
                                                                                                                                                caaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAAc   >  1:6569146/2‑151 (MQ=255)
                                                                                                                                                 aaagaaaaTGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCGTTTTACGGCCACCTGCCGGAAACg  >  1:479796/1‑151 (MQ=255)
                                                                                                                                                 aaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACg  >  1:1846951/1‑151 (MQ=255)
                                                                                                                                                 aaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTGTTCCGGCCAACTGCCCGTAACg  >  2:8195151/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACG  >  NC_008463/4659675‑4659970

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: