Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 732,660 G→A intergenic (+114/‑426) tyrS → / → PA14_RS03450 tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463732,6600GA97.3% 268.5 / ‑1.9 75intergenic (+114/‑426)tyrS/PA14_RS03450tyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (50/23);  total (50/25)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGA  >  NC_008463/732512‑732790
                                                                                                                                                    |                                                                                                                                  
agGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagt                                                                                                                                  <  2:8111323/151‑1 (MQ=255)
 ggCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtg                                                                                                                                 >  2:2571005/1‑151 (MQ=255)
 ggCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtg                                                                                                                                 >  1:4729119/1‑151 (MQ=255)
 ggCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtg                                                                                                                                 >  1:2755165/1‑151 (MQ=255)
  gCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                <  2:1791109/151‑1 (MQ=255)
  gCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                <  2:3320775/151‑1 (MQ=255)
   cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               >  1:6915340/1‑151 (MQ=255)
   cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               >  1:11456145/1‑151 (MQ=255)
   cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               >  2:3383186/1‑151 (MQ=255)
   cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               >  2:3391809/1‑151 (MQ=255)
     ttCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGa                                                                                                                             <  2:3394548/151‑1 (MQ=255)
      tCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGaa                                                                                                                            >  1:9626597/1‑151 (MQ=255)
      tCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGaa                                                                                                                            <  2:3477921/151‑1 (MQ=255)
       cgcgcgCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAg                                                                                                                           <  2:5486877/151‑1 (MQ=255)
       cgcgcgCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAg                                                                                                                           <  2:9610919/151‑1 (MQ=255)
        gcgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                          >  1:8270997/1‑151 (MQ=255)
        gcgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                          <  1:11255379/151‑1 (MQ=255)
        gcgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                          >  2:887540/1‑151 (MQ=255)
        gcgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACACCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGAGTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                          >  2:889083/1‑151 (MQ=255)
          gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTAGAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcgc                                                                                                                        >  1:4164921/1‑151 (MQ=255)
          gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCGATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGTAgcgc                                                                                                                        >  1:4159200/1‑151 (MQ=255)
            gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                      >  1:8994780/1‑151 (MQ=255)
            gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                      >  2:9904057/1‑151 (MQ=255)
            gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                      >  2:9904045/1‑151 (MQ=255)
            gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                      >  2:8165143/1‑151 (MQ=255)
             cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  1:10031063/1‑151 (MQ=255)
             cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  2:10531371/1‑151 (MQ=255)
             cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  2:694771/1‑151 (MQ=255)
             cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     <  2:11507989/151‑1 (MQ=255)
             cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCAACACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCATGGATTTACGGATTCTCTCCCAGGAGTGTTGAAGCGCTcc                                                                                                                     <  2:2322682/151‑1 (MQ=255)
              aTCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCg                                                                                                                    >  1:5929654/1‑151 (MQ=255)
               tCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCgg                                                                                                                   >  1:9947858/1‑151 (MQ=255)
               tCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCgg                                                                                                                   <  1:2752859/151‑1 (MQ=255)
                cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  >  1:7020135/1‑151 (MQ=255)
                cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  >  1:1868200/1‑151 (MQ=255)
                cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  <  2:5370053/151‑1 (MQ=255)
                 aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                 >  1:8370263/1‑151 (MQ=255)
                 aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                 <  1:9446176/151‑1 (MQ=255)
                  cGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAg                                                                                                                <  2:11483118/151‑1 (MQ=255)
                   gCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAgg                                                                                                               >  2:3408610/1‑151 (MQ=255)
                    cTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATACTGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgt                                                                                                              <  2:10223010/151‑1 (MQ=255)
                    cTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgt                                                                                                              >  1:3470917/1‑151 (MQ=255)
                    cTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGAAGgt                                                                                                              >  2:3681157/1‑151 (MQ=255)
                     tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                             >  2:9454435/1‑151 (MQ=255)
                     tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                             >  1:5043449/1‑151 (MQ=255)
                     tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                             >  2:1813155/1‑151 (MQ=255)
                      gAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttg                                                                                                            >  1:1247157/1‑151 (MQ=255)
                       aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  1:10299964/1‑151 (MQ=38)
                       aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  1:247032/1‑151 (MQ=38)
                       aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  1:269596/1‑151 (MQ=38)
                       aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  1:269532/1‑151 (MQ=38)
                       aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  1:6060686/1‑151 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          <  2:9154310/151‑1 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          >  2:8201213/1‑151 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          <  1:478798/151‑1 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          >  2:5523676/1‑151 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          >  2:1614571/1‑151 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          >  2:7770557/1‑151 (MQ=38)
                        aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          >  2:7770569/1‑151 (MQ=38)
                         aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCACTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTGGTTc                                                                                                         >  1:7807310/1‑151 (MQ=255)
                          gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcg                                                                                                        >  1:2445016/1‑151 (MQ=25)
                           cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       <  2:8193997/151‑1 (MQ=25)
                           cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  1:8915101/1‑151 (MQ=25)
                           cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  1:3606339/1‑151 (MQ=25)
                           cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  2:3352237/1‑151 (MQ=25)
                           cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  1:7987731/1‑151 (MQ=25)
                           cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  2:1411996/1‑151 (MQ=25)
                           cTGAATAAGTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       <  1:10214936/151‑1 (MQ=18)
                            tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCTCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                      >  2:4652586/1‑151 (MQ=21)
                            tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                      >  2:5239281/1‑151 (MQ=21)
                            tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                      <  1:1199059/151‑1 (MQ=21)
                             gAATAATTCTATCAAATAAGTGTTGACGCCAGTTTCTATCTCCCAATACTGCGCACCACTCCCGGCGACGACGACGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     <  2:11409178/151‑1 (MQ=255)
                              aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                    >  1:1346395/1‑151 (MQ=18)
                              aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                    >  1:1346513/1‑151 (MQ=18)
                               ataatTCTTTCAAATAAGTGTTGACGCCATTTTCTATCTCCCTATACTGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   <  2:4245038/151‑1 (MQ=12)
                               ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATCATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   <  2:824842/151‑1 (MQ=17)
                               ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   >  1:8391318/1‑151 (MQ=17)
                                taatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                  >  1:60413/1‑151 (MQ=14)
                                 aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 >  2:4093450/1‑151 (MQ=14)
                                 aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 <  2:2503756/151‑1 (MQ=14)
                                 aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 <  1:7993105/151‑1 (MQ=14)
                                 aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 <  1:7993106/151‑1 (MQ=14)
                                 aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 >  1:8845024/1‑151 (MQ=14)
                                                  gtagACGCCAGTTTCTATCTCCCTATAATGCGCACCACTACCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAgg                                                                                <  2:8245160/148‑1 (MQ=1)
                                                                                                                                ttcaCGGAATCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGa  <  1:6532655/148‑1 (MQ=1)
                                                                                                                                                    |                                                                                                                                  
AGGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGA  >  NC_008463/732512‑732790

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