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,8220AG96.3% 1285.1 / 24.6 353intergenic (+62/+37)alpR/PA14_RS21395transcriptional regulator AlpR/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (12/1);  new base G (144/195);  total (156/197)
Fisher's exact test for biased strand distribution p-value = 3.56e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 373 total aligned reads displayed.

CACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACG  >  NC_008463/4659674‑4659970
                                                                                                                                                    |                                                                                                                                                    
 aCCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaa                                                                                                                                                   <  1:5017468/151‑1 (MQ=255)
  cccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACTGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaa                                                                                                                                                  >  2:2485973/1‑151 (MQ=255)
   ccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaaa                                                                                                                                                 <  1:10761862/151‑1 (MQ=255)
   ccGGACGAGCGATACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaaa                                                                                                                                                 <  2:1222908/151‑1 (MQ=255)
    cGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAg                                                                                                                                                <  2:3462379/151‑1 (MQ=255)
    cGGACGAGCGCAACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAg                                                                                                                                                <  2:3280404/151‑1 (MQ=255)
     ggccGAGCGCTACGCGCCCCCGCAAGCCGCCGAGCACAACAACGTGATCGGCCGGGTGTTCTGGTCCTCTGTGCTGGTCTGACCCCTTCCCGCCCCGCCAGACGACGGGGCCGCGCTCCGGCGGAAGACACAGCCCGGGCACAGAAAAAgg                                                                                                                                               <  2:9018423/148‑1 (MQ=255)
     ggACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgg                                                                                                                                               <  1:5831382/151‑1 (MQ=255)
     ggACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgg                                                                                                                                               >  1:8112510/1‑151 (MQ=255)
      gACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAggg                                                                                                                                              >  1:7522194/1‑151 (MQ=255)
      gACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAggg                                                                                                                                              <  2:11655459/151‑1 (MQ=255)
        cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                            <  1:11222533/151‑3 (MQ=255)
         gAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcct                                                                                                                                           <  2:11501362/151‑4 (MQ=255)
          aGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcctg                                                                                                                                          >  1:7134872/1‑147 (MQ=255)
              cTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAtt                                                                                                                                     >  1:4223965/1‑151 (MQ=255)
                aCGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgc                                                                                                                                   <  1:1909340/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAACGAAAAAGGGCCTGCATtgct                                                                                                                                  <  2:7795699/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgct                                                                                                                                  <  2:4028419/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCANtgct                                                                                                                                  <  2:7699783/151‑1 (MQ=255)
                  gAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgctg                                                                                                                                 >  2:1336938/1‑151 (MQ=255)
                     ccccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAg                                                                                                                              >  1:6179359/1‑151 (MQ=255)
                      cccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGCCGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAgg                                                                                                                             <  1:905535/151‑1 (MQ=255)
                      cccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAgg                                                                                                                             >  1:289383/1‑151 (MQ=255)
                       ccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGc                                                                                                                            <  2:5534843/151‑1 (MQ=255)
                           gAAGCCGCCGAGCACCTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCtt                                                                                                                        <  2:1489334/151‑1 (MQ=255)
                           gAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCtt                                                                                                                        >  1:3681259/1‑151 (MQ=255)
                               ccgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCgg                                                                                                                    >  1:2388934/1‑151 (MQ=255)
                               ccgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCgg                                                                                                                    <  2:6967983/151‑1 (MQ=255)
                                cgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGt                                                                                                                   <  2:5598920/151‑1 (MQ=255)
                                cgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGt                                                                                                                   <  2:11349627/151‑1 (MQ=255)
                                  ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCACAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                 >  2:3227177/1‑151 (MQ=255)
                                  ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                 >  2:7050311/1‑151 (MQ=255)
                                  ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                 <  1:2053666/151‑1 (MQ=255)
                                   cgAGCACATCACCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATg                                                                                                                <  2:8674889/151‑1 (MQ=255)
                                    gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATgg                                                                                                               >  2:5101426/1‑151 (MQ=255)
                                    gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATgg                                                                                                               >  1:9921378/1‑151 (MQ=255)
                                      gCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCACAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGgcg                                                                                                             >  2:6185849/1‑151 (MQ=255)
                                      gCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGgcg                                                                                                             >  1:1401926/1‑151 (MQ=255)
                                      gCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGgcg                                                                                                             <  2:102346/151‑1 (MQ=255)
                                          aTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACt                                                                                                         <  1:8244011/151‑1 (MQ=255)
                                           tCAACGTGATCGGCCGGGTGTTCTGGTCCTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTc                                                                                                        <  1:2249763/151‑1 (MQ=255)
                                           tCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTGTCGGTATGGCGCACTc                                                                                                        >  2:8806394/1‑151 (MQ=255)
                                            cAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCa                                                                                                       <  2:8927184/151‑1 (MQ=255)
                                             aaCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCag                                                                                                      <  2:4728381/151‑1 (MQ=255)
                                              aCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagg                                                                                                     <  2:9781687/151‑1 (MQ=255)
                                              aCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagg                                                                                                     <  1:649746/151‑1 (MQ=255)
                                               cGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagga                                                                                                    >  1:3048084/1‑151 (MQ=255)
                                                gTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCaggag                                                                                                   <  1:3438252/151‑1 (MQ=255)
                                                      ggCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCg                                                                                             >  1:11589926/1‑151 (MQ=255)
                                                      ggCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCg                                                                                             <  2:9620312/151‑1 (MQ=255)
                                                      ggCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCGTTCGGTATGGCGCACTCAGGAGGATTCg                                                                                             >  2:6607383/1‑151 (MQ=255)
                                                          gggTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAAcc                                                                                         >  1:4929273/1‑151 (MQ=255)
                                                          gggTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAAcc                                                                                         >  2:2332113/1‑151 (MQ=255)
                                                             tgtTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTcc                                                                                      >  1:8552030/1‑151 (MQ=255)
                                                              gtTCTGGTACTCGGTGCGGGTCTGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAACGACAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCg                                                                                     <  2:2980914/151‑1 (MQ=255)
                                                                    gTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCt                                                                               <  2:6046460/151‑1 (MQ=255)
                                                                     tACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTc                                                                              >  1:3884515/1‑151 (MQ=255)
                                                                      aCTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCg                                                                             <  1:10398855/151‑1 (MQ=255)
                                                                       cTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCgg                                                                            <  1:5549302/151‑1 (MQ=255)
                                                                        tCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGCCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGt                                                                           <  2:6407819/151‑1 (MQ=255)
                                                                         cGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                          <  1:4306041/151‑1 (MQ=255)
                                                                         cGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                          <  2:1711237/151‑1 (MQ=255)
                                                                         cGGTGCTGGTCTGACCCCTTCCCGCCCCGCCCGACGACGGGGCCGCGATCCGGCGGAAGAAAAATCCCCGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATCGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtc                                                                          <  1:5597861/151‑2 (MQ=38)
                                                                          ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                         <  2:4605583/151‑1 (MQ=255)
                                                                            tGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGt                                                                       >  1:7963895/1‑151 (MQ=255)
                                                                              cTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     <  2:7448919/151‑1 (MQ=255)
                                                                                ggTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                   <  1:4086604/151‑1 (MQ=255)
                                                                                  tCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                 <  2:10097055/151‑1 (MQ=255)
                                                                                  tCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                 <  1:10068159/151‑1 (MQ=255)
                                                                                  tCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                 <  1:7802862/151‑1 (MQ=255)
                                                                                   cTGACCCCTGCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAg                                                                <  1:9297898/151‑1 (MQ=255)
                                                                                     gACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGAGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                              <  2:10892503/151‑1 (MQ=255)
                                                                                         ccgtCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGACGACATATCCCGGGCATAGACACTGGGCCCGCATGGCTGCAGGCCCTTGCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActct                                                          <  2:8552436/148‑1 (MQ=21)
                                                                                           ttCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                        <  1:9691994/151‑1 (MQ=255)
                                                                                               cgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                    <  2:10244186/151‑1 (MQ=255)
                                                                                                  cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 <  1:10003838/151‑1 (MQ=255)
                                                                                                  cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 >  2:8126717/1‑151 (MQ=255)
                                                                                                   ccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                >  2:10409631/1‑151 (MQ=255)
                                                                                                    cgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                               <  2:4241509/151‑1 (MQ=255)
                                                                                                         gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATg                                          >  1:739301/1‑151 (MQ=255)
                                                                                                         gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAGTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGTTTCGAACCTCCGACCGCTCGGGTCGTAGCCGAGTACTCTATCCAGCTGAGCTCTg                                          >  1:6635357/1‑151 (MQ=25)
                                                                                                            gacgGGGCCGCGTTCCGGCGGACGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGGTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGt                                       >  2:4323814/1‑151 (MQ=25)
                                                                                                                gggCCGCGTTCCGGCGGAAGAAACATCCCGGGCAAAGACAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCgt                                   <  1:4442092/151‑1 (MQ=21)
                                                                                                                gggCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCgt                                   <  1:1764459/151‑1 (MQ=255)
                                                                                                                   ccGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTc                                >  2:8373750/1‑151 (MQ=255)
                                                                                                                     gcgTTCCGGCGGAAGAAATATCCCGGGCAAAGAAAAAGGGCCTGCAGTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCGCGGTTCGTAGCCGCGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCgg                              >  2:5412672/1‑151 (MQ=21)
                                                                                                                      cgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGt                             >  2:3886334/1‑151 (MQ=255)
                                                                                                                         tCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCt                          >  1:1550825/1‑151 (MQ=255)
                                                                                                                          ccGGCGGAAGAAACATCCCGCGCAAAGCAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGAATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCtt                         <  2:1878752/151‑1 (MQ=21)
                                                                                                                           cggcggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCCGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttt                        <  1:3575167/151‑1 (MQ=255)
                                                                                                                              cggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTc                     <  1:10585832/151‑1 (MQ=255)
                                                                                                                               ggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTcc                    <  1:6921062/151‑1 (MQ=255)
                                                                                                                               ggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTcc                    <  2:8265071/151‑1 (MQ=255)
                                                                                                                                 aagaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCgg                  <  1:8649455/148‑1 (MQ=255)
                                                                                                                                  agaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActct                                                          >  1:8806309/3‑110 (MQ=255)
                                                                                                                                  agaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGc                 <  1:2664112/149‑1 (MQ=255)
                                                                                                                                     aaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTCTCCAGCTGAGCTATGAGTGCGTGTCGGTGCATATTCCGACCAc              >  2:6591925/1‑151 (MQ=21)
                                                                                                                                         tCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGc          <  1:4624605/151‑1 (MQ=255)
                                                                                                                                          cccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTATACGGCCCCCTggc         >  2:7307729/1‑149 (MQ=21)
                                                                                                                                                 caaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAAc   <  1:475421/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACG  >  NC_008463/4659674‑4659970

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