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.98.4% 1815.0 / 18.4 381intergenic (+11/+140)PA14_RS25260/PA14_RS25265tRNA‑Met/HigA family addiction module antitoxin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (187/186);  minor base T (5/1);  total (192/189)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 403 total aligned reads displayed.

GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201
                                                                                                                                                     |                                                                                                                                                 
 aaaaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgg                                                                                                                                                  <  2:8556905/150‑1 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                <  1:9245409/150‑1 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                <  2:11080875/150‑1 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGAGTCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                <  2:9235437/150‑1 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               >  1:2156672/1‑149 (MQ=255)
       aTCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggttag                                                                                                                                            >  1:7641007/1‑146 (MQ=255)
        tCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgc                                                                                                                                          >  1:5226638/1‑150 (MQ=255)
        tCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgc                                                                                                                                          >  1:11092804/1‑150 (MQ=255)
        tCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgc                                                                                                                                          >  2:4893880/1‑150 (MQ=255)
          tAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgcgc                                                                                                                                        <  2:654566/150‑1 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATACTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   <  1:9853588/150‑1 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGACACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   >  2:7646645/1‑150 (MQ=255)
                cAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTa                                                                                                                                  <  2:1468603/150‑1 (MQ=255)
                cAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGAGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTa                                                                                                                                  >  1:10818553/1‑150 (MQ=255)
                 aGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTaa                                                                                                                                 >  1:6890097/1‑150 (MQ=255)
                    ggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAccc                                                                                                                              <  2:5321880/150‑1 (MQ=255)
                       tGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCtt                                                                                                                           >  1:6269551/1‑150 (MQ=255)
                         aCAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCtttt                                                                                                                         <  1:5652880/150‑1 (MQ=255)
                           aGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttg                                                                                                                       <  2:9697514/150‑1 (MQ=255)
                            gCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTCGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGAGCAAGCTAACCACTttttgt                                                                                                                      >  2:10260103/1‑150 (MQ=255)
                            gCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgt                                                                                                                      >  2:10311907/1‑150 (MQ=255)
                            gCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgt                                                                                                                      >  2:4780339/1‑150 (MQ=255)
                             ccccGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtt                                                                                                                     >  1:939896/1‑150 (MQ=255)
                                cgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttttg                                                                                                                  <  2:773608/150‑1 (MQ=255)
                                  cccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGt                                                                                                                <  1:3001333/150‑1 (MQ=255)
                                   ccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTc                                                                                                               <  2:4369219/150‑1 (MQ=255)
                                    cgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATGTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCg                                                                                                              >  2:1787498/1‑150 (MQ=255)
                                     gagaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCgg                                                                                                             >  1:8928713/1‑150 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            >  1:4960059/1‑150 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTACAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            <  2:2389065/150‑1 (MQ=255)
                                       gacaCGGCACATGCCCGCCCCTCGGCTACAAAGCTCAGTCGGATAGAGCGCAGCCCACAAACTCCTGAGGTCCCAGGTTCAAGTCCCGGAGTAGCCACCACATTTTACAAGGGGTTAGCGCAAGCTAACCCCATATTGTCTGGGTCGGCa                                                                                                           <  2:9104897/147‑1 (MQ=255)
                                       gaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCa                                                                                                           <  2:6922575/150‑1 (MQ=255)
                                          aCGGACACTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATACTGCTGATGTCCCAGGTTCAAGTCCCTGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        <  1:10307229/150‑1 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        <  2:4742972/150‑1 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        <  1:7447758/150‑1 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        >  1:9515354/1‑150 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        >  1:1313872/1‑150 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        >  1:8423424/1‑150 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        <  1:5562313/150‑1 (MQ=255)
                                            ggAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgacta                                                                                                      >  1:3606269/1‑150 (MQ=255)
                                               aaaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaa                                                                                                   <  1:10002539/150‑1 (MQ=255)
                                               aaaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaa                                                                                                   <  1:9075717/150‑1 (MQ=255)
                                                 aTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaac                                                                                                 <  1:4762018/150‑1 (MQ=255)
                                                 aTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGAGCAAGCTAACCCCTTTTTGTGTTGGTCGGCAgactacaaac                                                                                                 >  2:8173172/1‑150 (MQ=255)
                                                    cgcgcCACACGGCTACATAGCTCAGTCGGTGAGACCGCAGCATTCATAATCCTGATGTCCCAGGTTCAAGTCCCGGTGTAGACACCAGATATTTCCAGGGGTTAGCGCACGCTACCCCCTTTTTGTTTTGGTCGGCAgactacaaaccga                                                                                              <  2:3447838/150‑1 (MQ=255)
                                                     gcgcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGTCAgactacaaaccgac                                                                                             >  1:4828805/1‑150 (MQ=255)
                                                     gcgcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgac                                                                                             <  1:6654819/150‑1 (MQ=255)
                                                     gcgcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgacgacaaaccgac                                                                                             <  1:7226999/150‑1 (MQ=255)
                                                       gcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgacta                                                                                           <  2:1059516/150‑1 (MQ=255)
                                                       gcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgacta                                                                                           <  1:431651/150‑1 (MQ=255)
                                                        ccACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATCATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactac                                                                                          <  2:292744/150‑1 (MQ=255)
                                                        ccACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactac                                                                                          <  2:291816/150‑1 (MQ=255)
                                                         cACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactaca                                                                                         >  1:455495/1‑150 (MQ=255)
                                                         cACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactaca                                                                                         <  2:7222734/150‑1 (MQ=255)
                                                            tCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaac                                                                                      >  2:6801650/1‑150 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   <  2:10904385/150‑1 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   >  1:10948984/1‑150 (MQ=255)
                                                                 tACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCt                                                                                 >  1:2580791/1‑150 (MQ=255)
                                                                 tACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAAGGCTGATGTCCCAGGTTCAAGTCCCGGAGTAGCCACCATATTTTTCAAGGGGTTAGCGCACGCGAACCCCTTTTAGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCt                                                                                 <  1:2759415/150‑1 (MQ=255)
                                                                 aacctAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTATTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCt                                                                                 <  2:10660942/146‑1 (MQ=255)
                                                                       gCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGccc                                                                           >  2:9522948/1‑150 (MQ=255)
                                                                       gCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGccc                                                                           >  2:3387982/1‑150 (MQ=255)
                                                                        cTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGcccc                                                                          <  2:7267377/150‑1 (MQ=255)
                                                                           aGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCtt                                                                       >  1:3658183/1‑150 (MQ=255)
                                                                              cGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTAc                                                                    >  2:3104923/1‑150 (MQ=255)
                                                                               ggTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTAcc                                                                   >  2:8059054/1‑150 (MQ=255)
                                                                                        gcagcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACgg                                                          <  2:249219/150‑1 (MQ=255)
                                                                                        gcagcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACgg                                                          >  2:6107941/1‑150 (MQ=255)
                                                                                         cagcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGt                                                         <  1:738268/150‑1 (MQ=255)
                                                                                          agcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGtt                                                        <  1:971035/150‑1 (MQ=255)
                                                                                          agcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGtt                                                        <  1:4121409/150‑1 (MQ=255)
                                                                                             aTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGc                                                     >  1:8295499/1‑150 (MQ=255)
                                                                                             aTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGc                                                     <  1:5941041/150‑1 (MQ=255)
                                                                                             aTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGc                                                     <  1:1968713/150‑1 (MQ=255)
                                                                                             aTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTATTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGAAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATATACGGTACGc                                                     >  1:8617172/1‑150 (MQ=255)
                                                                                                cATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATg                                                  >  1:2056396/1‑150 (MQ=255)
                                                                                                 aTAATGCTGATGTCCCATTTTCAAGTCCCGGTGTAGCCACAATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTTTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCACCCTTTACCGATTTACGGTTCGCATGa                                                 >  2:1350656/1‑150 (MQ=255)
                                                                                                  tAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaa                                                >  2:5669484/1‑150 (MQ=255)
                                                                                                    aTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAGTTACGGTTCGCATGaaaa                                              >  1:8091658/1‑150 (MQ=255)
                                                                                                          aTGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGcccc                                        >  2:5522227/1‑150 (MQ=255)
                                                                                                            gTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGc                                      >  1:1133189/1‑150 (MQ=255)
                                                                                                               ccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATc                                   >  1:1056687/1‑150 (MQ=255)
                                                                                                               ccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATc                                   <  1:1784414/150‑1 (MQ=255)
                                                                                                                cAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCt                                  >  2:1539634/1‑150 (MQ=255)
                                                                                                                  ggTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGc                                >  2:9324217/1‑150 (MQ=255)
                                                                                                                   gTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCg                               <  2:3783159/150‑1 (MQ=255)
                                                                                                                   gTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCg                               <  2:895755/150‑1 (MQ=255)
                                                                                                                   gTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCg                               <  1:401439/150‑1 (MQ=255)
                                                                                                                     tCAAGTCCCGGTGTAGCCACCATTTTTTTCAAGGGGTTAGCGCAAGCTCACCCCTTTTTGTTTTGGTCGGCAGACTACAATCCGACTACCAACCGTCTACGCCCGCTGTACCGATTTACGGTTCGCATGAAACAGCGCCGCATCTGCggg                             >  2:5607160/1‑150 (MQ=255)
                                                                                                                        aGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGct                          >  2:7466475/1‑150 (MQ=255)
                                                                                                                           cccGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCACTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTa                       <  2:4574344/150‑1 (MQ=255)
                                                                                                                               gtgtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCTGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGag                   >  2:7558652/1‑149 (MQ=255)
                                                                                                                                 gtcgCCACCATAATTATCACGGGGATAGCGCACGCTCACCCCTTTAGGTTTCGGTCGGCAGACTACAAACCGACTACCCACCGTCTACGCCCCCTGTCCCGATTTACGGTTCGCACGAAAAAGCCCCGCCTCTGCGGGGCTCTACGATCa                 <  2:6286671/147‑1 (MQ=255)
                                                                                                                                    gcaaccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACg              >  2:10077250/4‑150 (MQ=255)
                                                                                                                                    gCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACg              <  1:2496096/150‑1 (MQ=255)
                                                                                                                                    gCCACCATATTTTTCATGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACg              >  2:516092/1‑150 (MQ=255)
                                                                                                                                             tttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTCCGATCAACGTCACCGCTa     >  1:10642029/1‑150 (MQ=255)
                                                                                                                                             tttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTAGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTa     >  2:8639220/1‑150 (MQ=255)
                                                                                                                                              ttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCAACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCACGCATCTGCGGGGCTCTACGATCAACGTCACCGCTAt    <  2:4235837/150‑1 (MQ=255)
                                                                                                                                                ttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATcc  >  2:2782971/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 < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: