Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,909,637 (G)5→4 intergenic (‑133/+14) PA14_RS17840 ← / ← brnQ DUF599 domain‑containing protein/branched‑chain amino acid transport system II carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,909,6330G.99.2% 1895.0 / 6.6 390intergenic (‑129/+18)PA14_RS17840/brnQDUF599 domain‑containing protein/branched‑chain amino acid transport system II carrier protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (189/197);  total (192/198)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 411 total aligned reads displayed.

AATGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCG  >  NC_008463/3909484‑3909778
                                                                                                                                                     |                                                                                                                                                 
aaTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTg                                                                                                                                                   <  1:8452549/150‑1 (MQ=255)
aaTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTg                                                                                                                                                   >  1:438769/1‑150 (MQ=255)
  tGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTggg                                                                                                                                                 >  1:10302558/1‑150 (MQ=255)
   gCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACANGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTgggg                                                                                                                                                <  1:5636934/150‑1 (MQ=255)
   gCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTgggg                                                                                                                                                >  2:3827091/1‑150 (MQ=255)
     cGAGGTCGTCGGTCATGAAATCTCCCTCTTGAGCCGGCCGGGAATAGCCCCGCCGTTGTTTAGCGCCGTTGGGACGCCGCTCGCCGCGGCAAGGACAGGTATTTGGTCGGGTCATTATTATGGACGCCCGGAGTCGGCATTACTggtgtt                                                                                                                                              <  2:352358/150‑5 (MQ=255)
         ttCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGcc                                                                                                                                         <  2:7544211/150‑1 (MQ=255)
         ttCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGcc                                                                                                                                         <  2:4990748/150‑1 (MQ=255)
         ttCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGcc                                                                                                                                         <  1:2632860/150‑1 (MQ=255)
           cggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTg                                                                                                                                       <  2:1800608/150‑1 (MQ=255)
            ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGa                                                                                                                                      >  1:9310546/1‑150 (MQ=255)
            ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGa                                                                                                                                      <  2:7924925/150‑1 (MQ=255)
             gcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGAc                                                                                                                                     <  1:3560517/150‑1 (MQ=255)
              cggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACg                                                                                                                                    >  1:9542274/1‑150 (MQ=255)
               ggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGt                                                                                                                                   >  1:4815111/1‑150 (MQ=255)
                 tCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGtct                                                                                                                                 >  1:1274120/1‑150 (MQ=255)
                 tCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTGGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTGTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGgct                                                                                                                                 >  2:5282715/1‑150 (MQ=255)
                  cATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGtctc                                                                                                                                <  1:2687289/150‑1 (MQ=255)
                  cATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGtctc                                                                                                                                <  1:4631301/150‑1 (MQ=255)
                  cATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGtctc                                                                                                                                >  2:265577/1‑150 (MQ=255)
                     gAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAgg                                                                                                                             <  1:4914892/150‑1 (MQ=255)
                     gAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAgg                                                                                                                             >  1:1461489/1‑150 (MQ=255)
                      aaaaCTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGc                                                                                                                            >  2:4087628/1‑150 (MQ=255)
                      aaaaCTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGc                                                                                                                            <  1:135178/150‑1 (MQ=255)
                            cccTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGc                                                                                                                      <  1:595084/150‑1 (MQ=255)
                             ccTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCt                                                                                                                     <  2:6114180/150‑1 (MQ=255)
                              ctctTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCtt                                                                                                                    <  1:8697882/150‑1 (MQ=255)
                               tctTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCGTc                                                                                                                   >  2:2317093/1‑150 (MQ=255)
                                 ttGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCAGACGTCTCAGGCGACCGCTTcgc                                                                                                                 >  1:6250008/1‑150 (MQ=255)
                                    aGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAATGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTcgcgcg                                                                                                              >  2:5870871/1‑150 (MQ=255)
                                        ggccggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCtt                                                                                                          >  1:1102714/1‑150 (MQ=255)
                                         gccggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTg                                                                                                         >  2:3949533/1‑150 (MQ=255)
                                         gccggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTg                                                                                                         >  1:7744768/1‑150 (MQ=255)
                                            gggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGccc                                                                                                      <  2:1636190/150‑1 (MQ=255)
                                             ggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCca                                                                                                     <  2:2766659/150‑1 (MQ=255)
                                             ggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCca                                                                                                     >  2:10330274/1‑150 (MQ=255)
                                             ggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCGCATGAAATTGGCCGCCCGGAGGCGTCCATCCTGGGGTTCGCCAGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCca                                                                                                     >  2:737108/1‑150 (MQ=255)
                                                  cGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCagcagc                                                                                                <  1:7546159/150‑1 (MQ=255)
                                                   gCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGcc                                                                                               <  2:4794856/150‑1 (MQ=255)
                                                    ccccGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCg                                                                                              >  1:6067773/1‑150 (MQ=255)
                                                     cccGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGa                                                                                             <  1:66503/150‑1 (MQ=255)
                                                     cccGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGa                                                                                             >  2:10191208/1‑150 (MQ=255)
                                                     cccGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGa                                                                                             >  2:7184819/1‑150 (MQ=255)
                                                      ccgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGag                                                                                            >  1:2755860/1‑149 (MQ=255)
                                                      ccgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGAt                                                                                            <  2:4490714/150‑1 (MQ=255)
                                                       cgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATc                                                                                           >  1:452451/1‑150 (MQ=255)
                                                            ttgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCaaa                                                                                      <  2:8593022/150‑1 (MQ=255)
                                                            ttgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCaaa                                                                                      <  1:1542068/150‑1 (MQ=255)
                                                             tgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAAc                                                                                     <  1:8991758/150‑1 (MQ=255)
                                                             tgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACCTTCGTCGGCTCATGAACAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAAc                                                                                     <  1:366821/150‑1 (MQ=255)
                                                               tttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGa                                                                                   >  2:2508611/1‑150 (MQ=255)
                                                                ttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGAc                                                                                  <  2:955368/150‑1 (MQ=255)
                                                                ttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGAc                                                                                  <  1:326959/150‑1 (MQ=255)
                                                                 tCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACt                                                                                 <  2:6656955/150‑1 (MQ=255)
                                                                   gcgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTAATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGc                                                                               <  2:10638677/150‑1 (MQ=255)
                                                                    cgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGcc                                                                              >  1:1591739/1‑150 (MQ=255)
                                                                    cgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGcc                                                                              >  1:68068/1‑150 (MQ=255)
                                                                    cgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGcc                                                                              <  2:2273511/150‑1 (MQ=255)
                                                                      ccGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAg                                                                            >  1:3749876/1‑150 (MQ=255)
                                                                           gggACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCa                                                                       <  2:2909525/150‑1 (MQ=255)
                                                                           gggACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCa                                                                       <  1:3952127/150‑1 (MQ=255)
                                                                           gggACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCa                                                                       <  2:9036613/150‑1 (MQ=255)
                                                                                gccgcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCa                                                                  >  1:9897381/1‑150 (MQ=255)
                                                                                    cTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATgg                                                              <  1:543926/150‑1 (MQ=255)
                                                                                     tCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCCGCGCAATGGa                                                             >  2:506650/1‑150 (MQ=255)
                                                                                            ggCCAGGACCCGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTACGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgc                                                      <  2:8764169/150‑1 (MQ=255)
                                                                                             gCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgca                                                     >  1:6374009/1‑150 (MQ=255)
                                                                                              ccAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgcag                                                    <  1:2558779/150‑1 (MQ=255)
                                                                                                  gACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCagcagc                                                >  1:4996245/1‑150 (MQ=255)
                                                                                                   acacGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGcagcc                                               <  2:5698441/150‑1 (MQ=255)
                                                                                                    cacGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGcagcca                                              <  2:7721549/150‑1 (MQ=255)
                                                                                                     acGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGcagccag                                             >  2:7038888/1‑150 (MQ=255)
                                                                                                       gACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCagccagcc                                           >  2:8221551/1‑150 (MQ=255)
                                                                                                        aCTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAgccagcca                                          <  1:10554453/150‑1 (MQ=255)
                                                                                                         cTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTACTGGGGTTCGACTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCaa                                         >  2:4634430/1‑150 (MQ=255)
                                                                                                           ttCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAgg                                       <  1:5348372/150‑1 (MQ=255)
                                                                                                            tcgtcgGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGc                                      <  2:417128/150‑1 (MQ=255)
                                                                                                                cgGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTg                                  <  2:5500656/150‑1 (MQ=255)
                                                                                                                cgGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTg                                  >  2:9466504/1‑150 (MQ=255)
                                                                                                                 ggCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGa                                 >  2:3393585/1‑150 (MQ=255)
                                                                                                                 ggCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGa                                 <  2:9324095/150‑1 (MQ=255)
                                                                                                                  gCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGAt                                >  1:4117925/1‑150 (MQ=255)
                                                                                                                   cTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATc                               >  2:1140381/1‑150 (MQ=255)
                                                                                                                    tCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCg                              <  1:433448/150‑1 (MQ=255)
                                                                                                                       ttaaaaaGGCCGCCCGGAGGCGTCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCATCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTTCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCg                           >  1:10024991/3‑150 (MQ=255)
                                                                                                                         aaaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAg                         <  1:2045329/150‑1 (MQ=255)
                                                                                                                         aaaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAg                         >  1:1627987/1‑150 (MQ=255)
                                                                                                                          aaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGc                        <  2:6431019/150‑1 (MQ=255)
                                                                                                                           aaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCg                       <  1:3059115/150‑1 (MQ=255)
                                                                                                                             aGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCggg                     <  2:4353021/150‑1 (MQ=255)
                                                                                                                              ggCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGa                    >  1:2404938/1‑150 (MQ=255)
                                                                                                                               gccgccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAg                   >  1:800434/1‑150 (MQ=255)
                                                                                                                                ccgccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGc                  >  2:8076897/1‑150 (MQ=255)
                                                                                                                                 cgccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCt                 <  1:3412652/150‑1 (MQ=255)
                                                                                                                                  gccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCtt                <  1:5167598/150‑1 (MQ=255)
                                                                                                                                       gAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGa           >  1:7289410/1‑150 (MQ=255)
                                                                                                                                         ggcggcCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAg         <  1:8272056/150‑1 (MQ=255)
                                                                                                                                         ggcggcCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACAAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAg         <  1:709470/150‑1 (MQ=255)
                                                                                                                                            ggcCTTCCCGCGGGTCGCCTGACGACACACGCGACCGCATCGCGCGGCATGCCCACCCGCCGAGCGCAAACGACTGCCAGGACCCGCGCAATGGATCCCGGCAGCAGCCAGCCAAGGCTCCCCTCGGCGAGCGGGAGCTTGCCGAATACg      <  2:1882738/150‑1 (MQ=255)
                                                                                                                                                ttCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCg  >  2:2390955/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AATGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCG  >  NC_008463/3909484‑3909778

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: