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,6600GA93.2% 290.2 / 9.4 88intergenic (+114/‑426)tyrS/PA14_RS03450tyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (47/35);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 8.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGA  >  NC_008463/732513‑732790
                                                                                                                                                   |                                                                                                                                  
ggCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtg                                                                                                                                 >  1:4175901/1‑151 (MQ=255)
 gCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                <  2:7934483/151‑1 (MQ=255)
 gCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                >  1:2522484/1‑151 (MQ=255)
  cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               >  2:8653744/1‑151 (MQ=255)
  cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               >  2:5385794/1‑151 (MQ=255)
  cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtggt                                                                                                                               <  2:4926102/151‑1 (MQ=255)
   tttCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAgtggtg                                                                                                                              <  1:8331997/151‑1 (MQ=255)
   tttCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAgtggtg                                                                                                                              >  1:8336958/1‑151 (MQ=255)
   tttCGCGCGCATCACGCTGAAAGCTGAATAATTCGTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCCCTCCCGGCGACGACGAAGACTACATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAgtggtg                                                                                                                              <  2:9086559/151‑1 (MQ=255)
    ttCGCGCGCATCACGCTGAAAGATGAATATTTCTTTAAAATAAGTGTTGACTCCAGTTTCTATCTCCCTATAATGCGCAACACTCCCGGCGACGACGAAGAGTATATGCTTGAAAATCAAGGTTTTACGGATTATCTCCCTTGAGTGGTGa                                                                                                                             <  2:1529471/151‑1 (MQ=17)
     tCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGaa                                                                                                                            <  1:5600616/151‑1 (MQ=255)
     tCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGaa                                                                                                                            >  2:3372010/1‑151 (MQ=255)
     tCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGaa                                                                                                                            <  1:5849941/151‑1 (MQ=255)
       gcgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                          <  2:9484401/151‑1 (MQ=255)
        cgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcg                                                                                                                         >  2:10806865/1‑151 (MQ=255)
        cgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcg                                                                                                                         >  1:7028257/1‑151 (MQ=255)
         gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcgc                                                                                                                        >  1:6110162/1‑151 (MQ=255)
         gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCGATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcgc                                                                                                                        >  2:3248942/1‑151 (MQ=255)
          cgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCt                                                                                                                       <  2:3151679/151‑1 (MQ=255)
          cgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCt                                                                                                                       >  1:1212989/1‑151 (MQ=255)
           gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                      <  2:10743112/151‑1 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     <  2:5538161/151‑1 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     <  2:368195/151‑1 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  2:3028063/1‑151 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  2:3026668/1‑151 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  1:8424146/1‑151 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                     >  1:751252/1‑151 (MQ=255)
             aTCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCg                                                                                                                    <  1:5492417/151‑1 (MQ=255)
             aTCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCg                                                                                                                    <  1:5122492/151‑1 (MQ=255)
             aTCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCg                                                                                                                    <  1:3285081/151‑1 (MQ=255)
              tCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCgg                                                                                                                   >  2:6755667/1‑151 (MQ=255)
              tCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCgg                                                                                                                   <  1:4876750/151‑1 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTCTAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGAGTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  >  2:1082405/1‑151 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTATATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  <  2:9713297/151‑1 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  <  2:4496733/151‑1 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  <  1:214064/151‑1 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  >  1:3888091/1‑151 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAACTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                  <  2:9713136/151‑1 (MQ=255)
                aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                 >  1:3922243/1‑151 (MQ=255)
                aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                 <  2:9159294/151‑1 (MQ=255)
                aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                 <  1:3562964/151‑1 (MQ=255)
                aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                 >  2:11204394/1‑151 (MQ=255)
                  gCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAgg                                                                                                               >  1:4219322/1‑151 (MQ=255)
                   cTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATGTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgt                                                                                                              >  2:295507/1‑151 (MQ=255)
                    tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAGATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                             >  1:3637067/1‑151 (MQ=255)
                    tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                             >  1:6180674/1‑151 (MQ=255)
                    tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                             >  2:10620745/1‑151 (MQ=255)
                     gAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttg                                                                                                            <  2:8350383/151‑1 (MQ=255)
                      aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  2:3268355/1‑151 (MQ=38)
                      aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  2:3268339/1‑151 (MQ=38)
                      aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           <  2:1970035/151‑1 (MQ=38)
                      aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  1:7161343/1‑151 (MQ=38)
                      aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                           >  2:10627497/1‑151 (MQ=38)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          <  1:6894494/151‑1 (MQ=38)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                          >  1:9207458/1‑151 (MQ=38)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                         >  1:2153157/1‑151 (MQ=37)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                         <  1:5010404/151‑1 (MQ=37)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                         >  1:3442679/1‑151 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcg                                                                                                        <  2:1987794/151‑1 (MQ=25)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  1:3436680/1‑151 (MQ=25)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  2:11303911/1‑151 (MQ=25)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  2:10113225/1‑151 (MQ=25)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                       >  2:832167/1‑151 (MQ=25)
                           tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                      >  1:7152546/1‑151 (MQ=21)
                           tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                      >  2:5240558/1‑151 (MQ=21)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     >  1:6523622/1‑151 (MQ=18)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     >  1:6516994/1‑151 (MQ=18)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     >  2:3926202/1‑151 (MQ=18)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     >  2:3926373/1‑151 (MQ=18)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     >  2:8065848/1‑151 (MQ=18)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                     <  1:1644860/151‑1 (MQ=18)
                             aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                    >  2:5647981/1‑151 (MQ=18)
                             aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                    >  2:607351/1‑151 (MQ=18)
                             aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                    >  2:4771398/1‑151 (MQ=18)
                             aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                    <  1:6813463/151‑1 (MQ=18)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   >  1:7836935/1‑151 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   <  2:9707619/151‑1 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   <  1:7856551/151‑1 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   >  2:5088235/1‑151 (MQ=17)
                              ataatTCTTTCAAAGAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                   <  2:8711047/151‑1 (MQ=17)
                               taatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                  <  1:11343482/151‑1 (MQ=14)
                               taatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                  >  1:944027/1‑151 (MQ=14)
                               taatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                  <  2:8537881/151‑1 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 <  1:5552884/151‑1 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 <  1:4578638/151‑1 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 >  2:11251276/1‑151 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 <  1:2953614/151‑1 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGc                                                                                                 >  2:5501398/1‑151 (MQ=14)
                                     ttacAAAGAAGTGATGACGCCAGGTTCTAGCTCCATAGAATGCGCACCACGCCCGGCGACGACGAAGAGTAAAAGCTTGACAATCAAGGATTTACGGATTCTCGCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTc                                                                                            <  2:5010872/148‑1 (MQ=1)
                                           atcaGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTccc                                                                                      <  2:1251921/148‑1 (MQ=1)
                                                                tacctcCCGATCATGCGCACCACTCCCGGCGACGACGAAGCGTACATGCTTGCAAATCCAGGATTTACGGAGTCTCTCCCAGGGGTGGTGCAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGc                                                                 <  2:11107615/148‑1 (MQ=1)
                                                                          tactGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAgg                                                       <  1:8855335/148‑1 (MQ=1)
                                                                                                   cggagagTAAATGCTTGAAAATCAAGGATTTACGGATTCGCTCCCCGGGGTGGGGAAGCGCTCCGGCCGGTTGTTCGCCGCAGCGGTTCGGTCCCGAATGGGGCTCGAAACGAAGCTTCAGCGAGGTGCTTGACAGCGTATGTGATCGCTg                              >  2:2804710/4‑151 (MQ=1)
                                                                                                                               ttgaCGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGa  <  1:2941663/148‑1 (MQ=1)
                                                                                                                                                   |                                                                                                                                  
GGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGA  >  NC_008463/732513‑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: