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.3% 1218.0 / 18.5 340intergenic (+62/+37)alpR/PA14_RS21395transcriptional regulator AlpR/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (8/1);  new base G (126/201);  total (136/204)
Fisher's exact test for biased strand distribution p-value = 3.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 358 total aligned reads displayed.

CACCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACG  >  NC_008463/4659674‑4659970
                                                                                                                                                    |                                                                                                                                                    
 aCCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaa                                                                                                                                                   >  1:1122854/1‑151 (MQ=255)
 aCCCGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaa                                                                                                                                                   <  1:1709375/151‑1 (MQ=255)
   ccGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCaaagaaaaa                                                                                                                                                 <  2:6880152/151‑1 (MQ=255)
    cGGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAg                                                                                                                                                <  2:4932120/151‑1 (MQ=255)
      gACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAggg                                                                                                                                              <  1:6936547/151‑1 (MQ=255)
       aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                             <  1:2758833/151‑2 (MQ=255)
       aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                             >  1:3072886/1‑150 (MQ=255)
       aCGAGCGCTAAGATCCCCAGGAAGCCGCCGAGCACAACAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACACCTTCCCGCCCCGCCTGACGACGGGGCAGCGTTCCGGCTGAAGAAAAATCCCGGGCAAAGAAATAgggc                                                                                                                                             >  2:331457/1‑150 (MQ=255)
        cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                            >  2:6493676/1‑149 (MQ=255)
        cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                            <  1:3494144/151‑3 (MQ=255)
           gcgcTACGAGCCCCAGGAAGCTGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGc                                                                                                                                        <  2:3626872/151‑1 (MQ=255)
           gcgcTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGc                                                                                                                                        <  2:7878147/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgct                                                                                                                                  <  2:7397702/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgct                                                                                                                                  >  2:7646156/1‑151 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgct                                                                                                                                  <  1:1972469/151‑1 (MQ=255)
                 cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCGGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCCAAGAAAAAGGGCCTGCAGtgct                                                                                                                                  >  2:2097352/1‑151 (MQ=255)
                  gAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgctg                                                                                                                                 >  2:10730331/1‑151 (MQ=255)
                     ccccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAg                                                                                                                              <  1:9428880/151‑1 (MQ=255)
                       ccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGc                                                                                                                            <  1:3157338/151‑1 (MQ=255)
                        cAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGcc                                                                                                                           >  1:7057076/1‑151 (MQ=255)
                        cAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAAGCCCGGGCAAGGAAAATGGGCCTGCCGGGCTGCAGGcc                                                                                                                           >  2:9095902/1‑151 (MQ=255)
                         aGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGccc                                                                                                                          <  2:1507942/151‑1 (MQ=255)
                         aGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGccc                                                                                                                          <  1:6182809/151‑1 (MQ=255)
                          ggAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCt                                                                                                                         <  2:8091405/151‑1 (MQ=255)
                           gAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCtt                                                                                                                        <  1:5410038/151‑1 (MQ=255)
                           gAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCtt                                                                                                                        >  2:3262313/1‑151 (MQ=255)
                             aGCCGCCGAGCACATCACCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTc                                                                                                                      <  1:5796675/151‑1 (MQ=255)
                             aGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCGGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGTCATATCCCGGGAAATGAACATGGGCGTGCATTGCTGCAGGCCGtttg                                                                                                                      >  2:4991080/1‑150 (MQ=255)
                              gccgccGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCg                                                                                                                     <  1:6738415/151‑1 (MQ=255)
                              gccgccGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCg                                                                                                                     <  2:9193423/151‑1 (MQ=255)
                               ccgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCgg                                                                                                                    <  1:8381073/151‑1 (MQ=255)
                               ccgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCgg                                                                                                                    <  2:281382/151‑1 (MQ=255)
                                cgccgAGCACATCAACGTGCTCGGCCGGGTGGTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAACGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGt                                                                                                                   <  2:4361935/151‑1 (MQ=255)
                                cgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGt                                                                                                                   <  1:9312477/151‑1 (MQ=255)
                                  ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                 >  2:9794365/1‑151 (MQ=255)
                                     aGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGgc                                                                                                              <  1:8633443/151‑1 (MQ=255)
                                      gCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGgcg                                                                                                             <  1:3646019/151‑1 (MQ=255)
                                         caTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCAc                                                                                                          <  1:7544291/151‑1 (MQ=255)
                                          aTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACt                                                                                                         >  1:10820637/1‑151 (MQ=255)
                                            cAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCa                                                                                                       <  2:2297854/151‑1 (MQ=255)
                                            cAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCa                                                                                                       <  2:9529294/151‑1 (MQ=255)
                                               cGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagga                                                                                                    <  2:3958655/151‑1 (MQ=255)
                                                  gATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAggagga                                                                                                 <  1:960130/151‑1 (MQ=255)
                                                   aTCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGAt                                                                                                >  2:7488632/1‑151 (MQ=255)
                                                     cGGCCGGGTGTTCTGGTACTCGGTTCTGGTCTGACCCCATCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAACAAGCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTc                                                                                              <  2:9640659/151‑1 (MQ=255)
                                                       gCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTATCGGTATGGCGCACTCAGGAGGATTCGa                                                                                            <  1:1474023/151‑1 (MQ=255)
                                                        ccGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGAGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTGTCGGTATGGCGCACTCAGGAGGATTCGaa                                                                                           >  2:27582/1‑151 (MQ=255)
                                                         cGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAATAAAAACGGCCTGCAATGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAAc                                                                                          >  2:8243240/1‑151 (MQ=255)
                                                         cGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAAc                                                                                          >  1:9694372/1‑151 (MQ=255)
                                                            gtgtTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTc                                                                                       >  2:4533833/1‑151 (MQ=255)
                                                              gtTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCg                                                                                     <  2:9303529/151‑1 (MQ=255)
                                                              gtTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCg                                                                                     <  1:6527436/151‑1 (MQ=255)
                                                                 cTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGCACCTCCGAcc                                                                                  <  2:7031948/151‑1 (MQ=255)
                                                                 cTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGAcc                                                                                  <  1:8554075/151‑1 (MQ=255)
                                                                     tACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCCAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTGTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTc                                                                              >  2:6061928/1‑151 (MQ=39)
                                                                        tCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAACAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGt                                                                           <  1:1078186/151‑1 (MQ=255)
                                                                        acGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTACCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGt                                                                           <  1:10331925/150‑1 (MQ=255)
                                                                         cGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                          >  1:5592466/1‑151 (MQ=255)
                                                                          ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                         >  1:8056616/1‑151 (MQ=255)
                                                                           gTGCTGGTCTGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTAAGGAGGATTCGAACCTCCGACCGCTCGGTTCg                                                                        <  2:5234530/151‑1 (MQ=39)
                                                                             gCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                      <  1:8847174/151‑1 (MQ=255)
                                                                              cTGGTCTGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                     <  2:2872805/151‑1 (MQ=255)
                                                                                ggTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                   <  2:4629266/151‑1 (MQ=255)
                                                                                 gTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCg                                                                  >  1:5748273/1‑151 (MQ=255)
                                                                                    tGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGt                                                               <  1:4377957/151‑1 (MQ=255)
                                                                                    tGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGt                                                               <  2:2346013/151‑1 (MQ=255)
                                                                                     gACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                              <  1:1764494/151‑1 (MQ=255)
                                                                                        cccTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActc                                                           <  1:3553497/151‑1 (MQ=255)
                                                                                        cccTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActc                                                           >  2:3468976/1‑151 (MQ=255)
                                                                                         ccTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTCTGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActct                                                          >  2:5889382/1‑151 (MQ=255)
                                                                                           tacccgccccgcCTGACGCCGGGGCCGCGTTCCGGCGGAAGAACAATCCCGGGAAAAGAAAAAGGGCCTGCATTGCTGCAGGCCATTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                        <  2:11056338/149‑1 (MQ=25)
                                                                                             cccgccccgcCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                      <  2:7091090/151‑1 (MQ=255)
                                                                                               cgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCCAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                    >  1:7996840/1‑151 (MQ=255)
                                                                                                gccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAATGAAATAGGGCCTGCATTGCTGCTGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCGCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                   >  1:1608142/1‑151 (MQ=25)
                                                                                                 ccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                  <  1:8271212/151‑1 (MQ=255)
                                                                                                 ccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                  <  2:8542565/151‑1 (MQ=255)
                                                                                                  cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                 >  1:480687/1‑151 (MQ=255)
                                                                                                   ccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                <  1:10984082/151‑1 (MQ=255)
                                                                                                     gccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                              <  2:658016/151‑1 (MQ=255)
                                                                                                       cTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTa                                            <  1:3997857/151‑1 (MQ=255)
                                                                                                           cgacgGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAg                                        <  1:7690054/151‑1 (MQ=255)
                                                                                                              cgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGc                                     <  1:761580/151‑1 (MQ=255)
                                                                                                              cgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGc                                     >  2:8366899/1‑151 (MQ=255)
                                                                                                                    cgcgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCg                               <  1:7318748/151‑1 (MQ=255)
                                                                                                                     gcgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCgg                              >  2:1962612/1‑151 (MQ=255)
                                                                                                                      cgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGt                             >  1:125444/1‑151 (MQ=255)
                                                                                                                        ttCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGc                           >  2:8635524/1‑151 (MQ=255)
                                                                                                                             gcggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttttt                      <  2:7542953/151‑1 (MQ=255)
                                                                                                                                 aagaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCgg                  >  1:7859486/4‑151 (MQ=255)
                                                                                                                                  agaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGc                 >  1:4936995/3‑151 (MQ=255)
                                                                                                                                   gaaaacTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCGATGAGTGCGTGTCGGTGCTTTTTCCGGcc                <  2:10004985/150‑1 (MQ=255)
                                                                                                                                    acaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCa               <  1:10071410/149‑1 (MQ=255)
                                                                                                                                    aaacaTCCCGGGCAACCAAAAAGCGCCTGCATTGCTGCAGGCCCTATCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCACGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGACGGTGCATTTTCCGGCCa               <  2:4770585/151‑1 (MQ=18)
                                                                                                                                      aaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCAcc             <  1:21556/151‑1 (MQ=255)
                                                                                                                                          cccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGcc         <  2:10219528/151‑1 (MQ=255)
                                                                                                                                            cGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCgg       >  1:5176937/1‑151 (MQ=255)
                                                                                                                                            cGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCgg       <  2:3007326/151‑1 (MQ=255)
                                                                                                                                                gcaaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTaa    >  2:6334242/3‑151 (MQ=255)
                                                                                                                                                gcaaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCAGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTAAGGCCTCCTGCCGTTaa    >  2:6933285/3‑151 (MQ=21)
                                                                                                                                                 caaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAAc   >  2:5456997/2‑151 (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 < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: