Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,514,060 (G)5→4 intergenic (+15/+136) PA14_RS25260 → / ← PA14_RS25265 tRNA‑Met/HigA family addiction module antitoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,514,0560G.99.7% 1793.5 / 0.6 370intergenic (+11/+140)PA14_RS25260/PA14_RS25265tRNA‑Met/HigA family addiction module antitoxin
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (186/180);  total (188/182)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 407 total aligned reads displayed.

GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201
                                                                                                                                                     |                                                                                                                                                 
 aaacTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAACCTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAAGGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgg                                                                                                                                                  <  2:142439/150‑1 (MQ=255)
  aaaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggg                                                                                                                                                 <  2:2092393/150‑1 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               <  1:561353/150‑2 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               <  2:1929664/150‑2 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               >  2:7627549/1‑149 (MQ=255)
     ttATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtt                                                                                                                                              <  1:1950551/150‑3 (MQ=255)
      tctCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCGACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtta                                                                                                                                             <  2:582844/148‑4 (MQ=255)
      tATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtta                                                                                                                                             <  2:2528938/150‑4 (MQ=255)
       aTCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggttag                                                                                                                                            >  1:5264454/1‑146 (MQ=255)
       aTCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggttag                                                                                                                                            <  1:2892396/150‑5 (MQ=255)
         cTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgcg                                                                                                                                         >  2:8405998/1‑150 (MQ=255)
          tAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgcgc                                                                                                                                        <  1:9059197/150‑1 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      <  2:7396883/150‑1 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      >  2:83023/1‑150 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      <  1:24253/150‑1 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      <  1:7088026/150‑1 (MQ=255)
              cccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaatc                                                                                                                                    <  2:9824911/150‑3 (MQ=255)
              cccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGc                                                                                                                                    <  2:7748256/150‑1 (MQ=255)
              cccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGc                                                                                                                                    <  1:6031973/150‑1 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   <  1:8421566/150‑1 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACAGAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   <  2:4648709/150‑1 (MQ=255)
                        gACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttt                                                                                                                          >  2:9365807/1‑150 (MQ=255)
                             ccccGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCACGCTAACCCCAttttgtt                                                                                                                     >  2:6555087/1‑150 (MQ=255)
                             ccccGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtt                                                                                                                     <  2:8035496/150‑1 (MQ=255)
                              cccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttt                                                                                                                    <  1:5770971/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtttt                                                                                                                   >  1:3224246/1‑150 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACTCGGATACAGAGCTCAGTCGGTTGGAGCGCAGCATTCATACTGCTGAGGTCCCAGGTTGAAGTCCCGGTGTAGCCACCAGATTTTTGAAGGGGTTAGCGCAAGCTAACCCCTgtttggttt                                                                                                                   >  2:2838684/1‑150 (MQ=255)
                                cgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttttg                                                                                                                  >  1:2046672/1‑150 (MQ=255)
                                 gcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTgg                                                                                                                 <  2:10015286/150‑1 (MQ=255)
                                 gcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTgg                                                                                                                 <  1:2531107/150‑1 (MQ=255)
                                  cccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGt                                                                                                                >  1:4729944/1‑150 (MQ=255)
                                   ccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTc                                                                                                               >  1:1638147/1‑150 (MQ=255)
                                    cgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCg                                                                                                              >  1:5366781/1‑150 (MQ=255)
                                    cgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCg                                                                                                              <  1:10079092/150‑1 (MQ=255)
                                       gaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCa                                                                                                           >  2:11118874/1‑150 (MQ=255)
                                       gaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCa                                                                                                           >  2:3166814/1‑150 (MQ=255)
                                           cGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgact                                                                                                       >  1:7995288/1‑150 (MQ=255)
                                            ggAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgacta                                                                                                      <  1:3756814/150‑1 (MQ=255)
                                               aaaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaa                                                                                                   <  2:4722655/150‑1 (MQ=255)
                                               aaaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaa                                                                                                   >  2:7131195/1‑150 (MQ=255)
                                                aaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaa                                                                                                  >  1:1549649/1‑150 (MQ=255)
                                                 aTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaac                                                                                                 <  1:4017215/150‑1 (MQ=255)
                                                 aTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaac                                                                                                 >  2:4494351/1‑150 (MQ=255)
                                                     gcgcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgac                                                                                             <  1:10778177/150‑1 (MQ=255)
                                                      cgcCACTCGGCTACATAGCTCAGTCGGTTTGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgact                                                                                            <  2:4389501/150‑1 (MQ=255)
                                                      cgcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgact                                                                                            <  2:7045198/150‑1 (MQ=255)
                                                        ccACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactac                                                                                          >  1:9247646/1‑150 (MQ=255)
                                                           cTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaa                                                                                       >  2:2931008/1‑150 (MQ=255)
                                                            tCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaac                                                                                      <  1:9781408/150‑1 (MQ=255)
                                                            tCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTGTTTTGGTCGGCAgactacaaaccgactacaaacc                                                                                     <  1:9335760/150‑1 (MQ=255)
                                                              ggCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATACTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGactacaaaccgactacaaaccg                                                                                    <  1:4486576/150‑1 (MQ=255)
                                                              ggCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGactacaaaccgactacaaaccg                                                                                    >  1:7816040/1‑150 (MQ=255)
                                                                  aCATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTa                                                                                >  1:1016438/1‑150 (MQ=255)
                                                                     tAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGAGGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGc                                                                             <  2:8183689/150‑1 (MQ=255)
                                                                      aGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGcc                                                                            >  1:10582572/1‑150 (MQ=255)
                                                                         tCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGccccc                                                                         <  1:8721933/150‑1 (MQ=255)
                                                                           aGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCtt                                                                       >  1:2752828/1‑150 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      >  1:5184034/1‑150 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      >  2:10992836/1‑150 (MQ=255)
                                                                                gTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCg                                                                  >  2:8658259/1‑150 (MQ=255)
                                                                                    gagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAttt                                                              >  1:8537295/1‑150 (MQ=255)
                                                                                    gagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAttt                                                              >  1:1263266/1‑150 (MQ=255)
                                                                                      gcgcAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTAc                                                            <  2:1644531/150‑1 (MQ=255)
                                                                                      gcgcAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTAc                                                            <  1:7374584/150‑1 (MQ=255)
                                                                                      gcgcAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATGTAc                                                            >  2:5120283/1‑150 (MQ=255)
                                                                                         cagcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGt                                                         <  2:6833457/150‑1 (MQ=255)
                                                                                          agcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGtt                                                        >  1:6247087/1‑150 (MQ=255)
                                                                                           gcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTc                                                       <  2:3639017/150‑1 (MQ=255)
                                                                                              ttCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCa                                                    >  1:1401216/1‑150 (MQ=255)
                                                                                                cATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATg                                                  <  2:397515/150‑1 (MQ=255)
                                                                                                 aTAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGa                                                 >  1:6695686/1‑150 (MQ=255)
                                                                                                 aTAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGa                                                 <  1:3168914/150‑1 (MQ=255)
                                                                                                  tAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaa                                                <  1:4682488/150‑1 (MQ=255)
                                                                                                   aaTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaa                                               <  1:8173106/150‑1 (MQ=255)
                                                                                                    aTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaa                                              >  2:277406/1‑150 (MQ=255)
                                                                                                     tGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaaa                                             <  1:9639955/150‑1 (MQ=255)
                                                                                                     tGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaaa                                             >  2:8862443/1‑150 (MQ=255)
                                                                                                     tGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaaa                                             <  1:8339604/150‑1 (MQ=255)
                                                                                                       cTGATGTCCCAGGTTCATGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCAAACCCCTTTTTGTTTTGGTCTGCAGACTACAAACCTACTACAAACCGACAACGCCCCCATTACCGATTTACTTTTCTCATGaaaaatc                                           >  2:9214336/1‑148 (MQ=255)
                                                                                                       cTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGc                                           <  1:4819712/150‑1 (MQ=255)
                                                                                                        tGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGcc                                          <  2:4171284/150‑1 (MQ=255)
                                                                                                        tGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTGTACCGAGTTACGGTTCGCATGAAAAAGcc                                          >  2:6653307/1‑150 (MQ=255)
                                                                                                          aTGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGcccc                                        >  2:10899699/1‑150 (MQ=255)
                                                                                                             tCCCAGGTTCACGTCCGGGGGTAGCCACCATCTTTTTCAAGGGGGTAGCGCAAGCTAACCCGTTTTTGTTTTGGTCGGCAGACTACAAACCGACTAGAAACCGTCTACGCCCCCTGTACCGATTTACGGTTCGCATGAAAAAGCCCCGCa                                     <  2:10467241/150‑1 (MQ=255)
                                                                                                             tCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCa                                     >  2:2138418/1‑150 (MQ=255)
                                                                                                             tCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCGAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCa                                     <  1:6797458/150‑1 (MQ=255)
                                                                                                               ccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATc                                   <  1:8679287/150‑1 (MQ=255)
                                                                                                               ccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCATTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATc                                   <  1:6419138/150‑1 (MQ=255)
                                                                                                                 aGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTg                                 <  2:9000689/150‑1 (MQ=255)
                                                                                                                   gTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACACACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCg                               >  2:2299839/1‑150 (MQ=255)
                                                                                                                   gTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCg                               >  1:4165699/1‑150 (MQ=255)
                                                                                                                    ttCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgg                              <  2:5765150/150‑1 (MQ=255)
                                                                                                                     tCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAGACGACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCggg                             <  2:612594/150‑1 (MQ=255)
                                                                                                                           cccGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTa                       <  2:10682471/150‑1 (MQ=255)
                                                                                                                              ggTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGa                    >  2:5518829/1‑150 (MQ=255)
                                                                                                                               gtgtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGAt                   <  1:10942440/150‑1 (MQ=255)
                                                                                                                                tgtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATc                  >  1:1764726/1‑150 (MQ=255)
                                                                                                                                     ccaccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGt             <  1:3565212/150‑1 (MQ=255)
                                                                                                                                         caTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAATGCCCCGCATCTGCGGGGCTCTACGATCAACGTCAcc         >  1:2994373/1‑150 (MQ=255)
                                                                                                                                               tttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATc   >  2:3730218/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201

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