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.5% 1356.4 / 19.3 363intergenic (+62/+37)alpR/PA14_RS21395transcriptional regulator AlpR/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (8/1);  new base G (133/219);  total (143/220)
Fisher's exact test for biased strand distribution p-value = 2.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 390 total aligned reads displayed.

CACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACG  >  NC_008463/4659674‑4659970
                                                                                                                                                    |                                                                                                                                                    
cACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaa                                                                                                                                                    <  1:6146321/151‑1 (MQ=255)
 aCCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaa                                                                                                                                                   <  1:7719821/151‑1 (MQ=255)
 aCCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaa                                                                                                                                                   >  2:781696/1‑151 (MQ=255)
  cccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaa                                                                                                                                                  <  2:9606979/151‑1 (MQ=255)
   ccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGATAAATCCCGGGCaaagaatat                                                                                                                                                 >  1:4426650/1‑148 (MQ=255)
   ccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaaa                                                                                                                                                 <  1:4811001/151‑1 (MQ=255)
   ccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaaa                                                                                                                                                 <  1:5439944/151‑1 (MQ=255)
    cGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAg                                                                                                                                                <  2:9446242/151‑1 (MQ=255)
     ggACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgg                                                                                                                                               <  1:10011923/151‑1 (MQ=255)
       aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                             <  1:7319980/151‑2 (MQ=255)
         gAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcct                                                                                                                                           <  1:3129747/151‑4 (MQ=255)
         gAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACACCCCCGTGATCGGCCGGGAGTTCTGGTACTCGCTGCTGGGCTGACCCCTTCCCGCCCCGCCTGACGACGGCGCCGCGTTCCGGCGGAAGAACCATCCCGGGCAAAGAACAAgggcct                                                                                                                                           <  2:5936738/151‑4 (MQ=255)
           gcgcTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCAGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGc                                                                                                                                        <  2:6717418/151‑1 (MQ=255)
             gcTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAt                                                                                                                                      >  1:3490526/1‑151 (MQ=255)
             gcTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAt                                                                                                                                      <  2:8462335/151‑1 (MQ=255)
             gcTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAt                                                                                                                                      <  1:9338334/151‑1 (MQ=255)
              cTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAtt                                                                                                                                     <  2:4821927/151‑1 (MQ=255)
               tACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtg                                                                                                                                    <  2:2454680/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgct                                                                                                                                  <  1:990623/151‑1 (MQ=255)
                   aGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgctgc                                                                                                                                >  1:3606507/1‑151 (MQ=255)
                      cccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAgg                                                                                                                             <  1:825500/151‑1 (MQ=255)
                        cAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGcc                                                                                                                           >  2:8555245/1‑151 (MQ=255)
                         aGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGccc                                                                                                                          >  1:6535149/1‑151 (MQ=255)
                           gAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCtt                                                                                                                        <  2:6309313/151‑1 (MQ=255)
                           gAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCtt                                                                                                                        <  1:569710/151‑1 (MQ=255)
                             aGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTc                                                                                                                      <  2:2008520/151‑1 (MQ=255)
                                cgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGt                                                                                                                   <  2:10430001/151‑1 (MQ=255)
                                cgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGt                                                                                                                   <  1:8086007/151‑1 (MQ=255)
                                 gccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTa                                                                                                                  >  2:4074162/1‑151 (MQ=255)
                                 gccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTa                                                                                                                  <  1:6643304/151‑1 (MQ=255)
                                  ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                 >  2:7500530/1‑151 (MQ=255)
                                  ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGACTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                 <  1:802701/151‑1 (MQ=255)
                                   cgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATg                                                                                                                <  2:3233921/151‑1 (MQ=255)
                                    gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTCTgg                                                                                                               >  1:7050706/1‑151 (MQ=255)
                                    gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATgg                                                                                                               <  1:7807708/151‑1 (MQ=255)
                                    gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATgg                                                                                                               <  1:1850903/151‑1 (MQ=255)
                                    gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCTTTCGGTATgg                                                                                                               >  1:5892200/1‑151 (MQ=255)
                                       cacaTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGgcgc                                                                                                            >  1:7144837/1‑151 (MQ=255)
                                                   aTCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGAt                                                                                                >  2:1336882/1‑151 (MQ=255)
                                                   aTCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCAGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGAt                                                                                                <  2:7394882/151‑1 (MQ=255)
                                                     cGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTc                                                                                              <  1:965516/151‑1 (MQ=255)
                                                        ccGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGaa                                                                                           >  1:1634749/1‑151 (MQ=255)
                                                           ggTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCt                                                                                        <  1:3306585/151‑1 (MQ=255)
                                                            gtgtTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTc                                                                                       >  2:7128398/1‑151 (MQ=255)
                                                                     tACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTc                                                                              <  1:2413496/151‑1 (MQ=255)
                                                                      aCTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCg                                                                             <  2:7993773/151‑1 (MQ=255)
                                                                          ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                         <  1:6913309/151‑1 (MQ=255)
                                                                              cTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     <  1:2216947/151‑1 (MQ=255)
                                                                              cTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     <  1:2915103/151‑1 (MQ=255)
                                                                              cTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCATTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     <  2:10648683/151‑1 (MQ=255)
                                                                               tGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGc                                                                    <  1:421093/151‑1 (MQ=255)
                                                                               tGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGc                                                                    <  1:10068883/151‑1 (MQ=255)
                                                                                     gACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGATGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                              <  2:7715421/151‑1 (MQ=38)
                                                                                      aCCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAc                                                             <  1:4629144/151‑1 (MQ=255)
                                                                                           ttCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                        <  2:181263/151‑1 (MQ=255)
                                                                                           ttCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                        >  1:212092/1‑151 (MQ=255)
                                                                                             cccgccccgcCTGACGACGGGGCCGCGTTCCGGCGGAAGACAAATCCCGGGCAATGAAAATGGGCCTGCATTGCTGCAGGCCCTTGCGGTAGGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTCTcc                                                      >  1:9441275/1‑151 (MQ=25)
                                                                                             cccgccccgcCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                      <  2:9759632/151‑1 (MQ=255)
                                                                                              ccgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcct                                                     >  2:5018929/1‑150 (MQ=255)
                                                                                              ccgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATACCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                     >  1:227874/1‑151 (MQ=255)
                                                                                               cgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAGGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                    <  2:1457352/151‑1 (MQ=255)
                                                                                                gccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                   >  2:2137074/1‑151 (MQ=255)
                                                                                                 ccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                  <  1:5096846/151‑1 (MQ=255)
                                                                                                  cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 <  2:3791514/151‑1 (MQ=255)
                                                                                                  cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 <  1:190434/151‑1 (MQ=255)
                                                                                                  cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 >  1:2625033/1‑151 (MQ=255)
                                                                                                    cgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                               >  1:2949236/1‑151 (MQ=255)
                                                                                                        tGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                           >  2:7926136/1‑151 (MQ=255)
                                                                                                        tGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                           >  2:4229935/1‑151 (MQ=255)
                                                                                                         gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATg                                          >  1:7583750/1‑151 (MQ=255)
                                                                                                         gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATg                                          <  2:9223766/151‑1 (MQ=255)
                                                                                                           cgacgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAg                                        >  1:107442/1‑151 (MQ=255)
                                                                                                              cgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGc                                     >  1:4937213/1‑151 (MQ=255)
                                                                                                               ggggCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAATAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCg                                    >  2:10150286/1‑151 (MQ=255)
                                                                                                               ggggCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCCGGCCCTGTCGGTATGGCGCACTCAGGAGGATGCGAACCTCCGACCGCTCGGTGCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCg                                    >  1:55611/1‑151 (MQ=21)
                                                                                                               aggtccGAGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCAGCAATGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCg                                    >  2:6578128/5‑151 (MQ=21)
                                                                                                                 ggCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCgtg                                  <  1:10282724/151‑1 (MQ=255)
                                                                                                                   ccGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTc                                >  1:11165043/1‑151 (MQ=255)
                                                                                                                      cgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGt                             <  2:10888857/151‑1 (MQ=255)
                                                                                                                       gTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGGCGGTg                            >  2:8853373/1‑151 (MQ=255)
                                                                                                                       gTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGGCGGTg                            >  1:8786014/1‑151 (MQ=255)
                                                                                                                        ttCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGc                           <  2:3348586/151‑1 (MQ=255)
                                                                                                                        ttCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGc                           <  1:1716469/151‑1 (MQ=255)
                                                                                                                           cggcggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttt                        >  2:4634656/1‑151 (MQ=255)
                                                                                                                              cggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttttac                     <  1:910602/151‑3 (MQ=255)
                                                                                                                                gaagaaAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCg                   <  1:10415253/151‑1 (MQ=255)
                                                                                                                                  agaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCGGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGc                 <  2:820058/149‑1 (MQ=255)
                                                                                                                                   gaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGcc                <  1:378932/150‑1 (MQ=255)
                                                                                                                                    aaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCa               >  2:5643859/1‑151 (MQ=255)
                                                                                                                                      aaaTCCCGGGCAAAGAAACAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCAcc             <  2:6098875/151‑1 (MQ=255)
                                                                                                                                      aaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCAcc             <  2:524932/151‑1 (MQ=255)
                                                                                                                                        aTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTg           <  2:2760344/151‑1 (MQ=255)
                                                                                                                                        aTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCATTTTCCGGCCACCTg           >  2:1710326/1‑151 (MQ=255)
                                                                                                                                           ccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCg        <  2:10538978/151‑1 (MQ=255)
                                                                                                                                             gggCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTACGGCCACCTGCCGGt      >  2:11171762/1‑151 (MQ=255)
                                                                                                                                                gcaaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTACGTCCACCTGCCGGTaa    >  1:10319979/3‑151 (MQ=255)
                                                                                                                                                 caaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAAc   >  1:10592164/2‑151 (MQ=255)
                                                                                                                                                  aaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGATAGCGTGTCGGTGCTTTTTAAGTCCACAAGCCGGTAACg  >  1:154449/1‑151 (MQ=21)
                                                                                                                                                  aaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACg  <  1:10842313/151‑1 (MQ=255)
                                                                                                                                                  aaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGAATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACg  <  2:3606507/151‑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: