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.98.9% 1829.9 / 10.0 378intergenic (‑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/2);  new base . (196/178);  total (198/180)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 403 total aligned reads displayed.

AATGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCG  >  NC_008463/3909484‑3909778
                                                                                                                                                     |                                                                                                                                                 
aaTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTg                                                                                                                                                   >  1:6120922/1‑150 (MQ=255)
 aTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTgg                                                                                                                                                  <  1:7101745/150‑1 (MQ=255)
  tcccGAGCTCGGCGGTCCAGAAACCTCCCACCCGAGCCGGCCGGGAAGCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTggg                                                                                                                                                 <  2:9044779/148‑1 (MQ=255)
         ttCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGcc                                                                                                                                         >  2:5751385/1‑150 (MQ=255)
           cggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTg                                                                                                                                       <  1:10088019/150‑1 (MQ=255)
           cggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTg                                                                                                                                       >  1:5399722/1‑150 (MQ=255)
            ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGa                                                                                                                                      >  2:8557240/1‑150 (MQ=255)
              cggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACg                                                                                                                                    >  2:2321844/1‑150 (MQ=255)
               ggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGt                                                                                                                                   <  1:9994549/150‑1 (MQ=255)
               ggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGt                                                                                                                                   >  2:6058604/1‑150 (MQ=255)
                  cATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGtctc                                                                                                                                <  1:6384267/150‑1 (MQ=255)
                   aTGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCa                                                                                                                               <  2:5897371/150‑1 (MQ=255)
                   aTGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCa                                                                                                                               >  2:4527402/1‑150 (MQ=255)
                   aTGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCa                                                                                                                               >  1:8160620/1‑150 (MQ=255)
                   aTGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCa                                                                                                                               >  1:1835976/1‑150 (MQ=255)
                    tGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAg                                                                                                                              >  2:2045678/1‑150 (MQ=255)
                    agAACACTCCCTCCCGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAg                                                                                                                              <  2:8975387/149‑1 (MQ=255)
                     gAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAAGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAgg                                                                                                                             <  2:8633021/150‑1 (MQ=255)
                       aaaCTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCg                                                                                                                           <  2:10991044/150‑1 (MQ=255)
                          cTCCCTCTTGAGCCGGCCGGGACTCGCCCCGCCGGTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGAcc                                                                                                                        >  2:6496376/1‑150 (MQ=255)
                                  tGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTcgcg                                                                                                                <  1:3564087/150‑1 (MQ=255)
                                   gAGCCGGCCGGCAATCGCCCCGCCGGCGTGTCGCGCCGATCGGACGCCGCTCGCCGCGGCCAGGACACGACTATCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTcgcgc                                                                                                               <  2:7581079/150‑1 (MQ=255)
                                     gccggccgGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCgg                                                                                                             >  1:3902009/1‑150 (MQ=255)
                                      ccggccggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGGCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGc                                                                                                            >  2:299200/1‑150 (MQ=255)
                                           cggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGcc                                                                                                       >  2:5629532/1‑150 (MQ=255)
                                           cggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGcc                                                                                                       <  2:83934/150‑1 (MQ=255)
                                            gggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGccc                                                                                                      <  1:1282859/150‑1 (MQ=255)
                                            gggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGccc                                                                                                      <  1:10691442/150‑1 (MQ=255)
                                            gggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGccc                                                                                                      >  1:223791/1‑150 (MQ=255)
                                            gggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTAGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGccc                                                                                                      <  1:6016337/150‑1 (MQ=255)
                                                 tCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagcag                                                                                                 <  1:2136495/150‑1 (MQ=255)
                                                 tCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagcag                                                                                                 <  2:1836515/150‑1 (MQ=255)
                                                  cGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGACAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCagcagc                                                                                                <  2:4361155/150‑1 (MQ=255)
                                                   gCCCCGCCGTTGTTTCGCGCCGTTGGGTCGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCAAGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGcc                                                                                               <  1:4572786/150‑1 (MQ=255)
                                                    ccccGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAATGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCg                                                                                              >  1:5062350/1‑150 (MQ=255)
                                                      ccgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGAt                                                                                            >  1:8742098/1‑150 (MQ=255)
                                                      ccgccgTTGTTTCGCGCCGTGGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGAt                                                                                            <  2:9377287/150‑1 (MQ=255)
                                                       cgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATc                                                                                           <  2:4835813/150‑1 (MQ=255)
                                                        gccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCg                                                                                          >  2:1660772/1‑150 (MQ=255)
                                                          cgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCa                                                                                        >  2:2566935/1‑150 (MQ=255)
                                                            ttgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCaaa                                                                                      <  2:10370202/150‑1 (MQ=255)
                                                               tttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGa                                                                                   >  2:1556023/1‑150 (MQ=255)
                                                               tttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCAAGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGa                                                                                   <  2:7343593/150‑1 (MQ=255)
                                                               ttgcgcgCCGTCGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGa                                                                                   <  2:6916212/147‑1 (MQ=255)
                                                                ttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGAc                                                                                  <  1:4846908/150‑1 (MQ=255)
                                                                ttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGACCAGCAGCCGATCGCAAACGAc                                                                                  <  2:5239832/150‑1 (MQ=255)
                                                                 tCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACt                                                                                 >  2:3491344/1‑150 (MQ=255)
                                                                 tCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACt                                                                                 <  1:2978981/150‑1 (MQ=255)
                                                                   gcgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGc                                                                               >  1:3338184/1‑150 (MQ=255)
                                                                    cgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGcc                                                                              >  2:6400767/1‑150 (MQ=255)
                                                                     gcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCa                                                                             <  2:9685694/150‑1 (MQ=255)
                                                                      ccGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCCGCAGCCGATCGCAAACGACTGCCAg                                                                            >  2:10733028/1‑150 (MQ=255)
                                                                      ccGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAg                                                                            >  2:10574036/1‑150 (MQ=255)
                                                                      ccGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAg                                                                            >  1:561089/1‑150 (MQ=255)
                                                                         ttGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGAc                                                                         >  2:2392103/1‑150 (MQ=255)
                                                                              aCGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAgcg                                                                    >  1:6787190/1‑150 (MQ=255)
                                                                               cgccgcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAATGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAgcgc                                                                   >  2:3742323/1‑150 (MQ=255)
                                                                                gccgcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCa                                                                  >  2:6180721/1‑150 (MQ=255)
                                                                                 ccgcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCaa                                                                 <  2:10100957/150‑1 (MQ=255)
                                                                                   gcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATg                                                               <  2:3339434/150‑1 (MQ=255)
                                                                                    cTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATgg                                                              >  1:11137115/1‑150 (MQ=255)
                                                                                       gccgcGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATa                                                           <  1:2032021/150‑1 (MQ=255)
                                                                                       gccgcGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATa                                                           <  2:8163139/150‑1 (MQ=255)
                                                                                        ccgcGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATAc                                                          >  2:5100882/1‑150 (MQ=255)
                                                                                         cgcgGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATAcc                                                         <  2:9196765/150‑1 (MQ=255)
                                                                                         cgcgGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATAcc                                                         >  1:8227026/1‑150 (MQ=255)
                                                                                          gcgGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCg                                                        <  1:9816553/150‑1 (MQ=255)
                                                                                            ggCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATCCCGgc                                                      >  2:4689131/1‑150 (MQ=255)
                                                                                               cagtaaacGAATTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgcagc                                                   <  1:3058670/144‑1 (MQ=255)
                                                                                                   acacGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGcagcc                                               <  2:10466965/150‑1 (MQ=255)
                                                                                                    cacGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGcagcca                                              >  1:2623167/1‑150 (MQ=255)
                                                                                                       gACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCagccagcc                                           >  1:1685453/1‑150 (MQ=255)
                                                                                                         cTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCACACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCaa                                         >  2:10243117/1‑150 (MQ=255)
                                                                                                              gtcgGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGctc                                    >  2:8790940/1‑150 (MQ=255)
                                                                                                              gtcgGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGctc                                    >  1:3176661/1‑150 (MQ=255)
                                                                                                                  gCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGAt                                <  2:3907530/150‑1 (MQ=255)
                                                                                                                   cTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATc                               >  1:10880453/1‑150 (MQ=255)
                                                                                                                   cTCATCAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGAACGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATc                               <  1:1882261/150‑1 (MQ=255)
                                                                                                                     cATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCgg                             >  1:4251575/1‑150 (MQ=255)
                                                                                                                      aTGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGc                            >  1:8488162/1‑150 (MQ=255)
                                                                                                                      aTGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGc                            >  1:7686626/1‑150 (MQ=255)
                                                                                                                        gAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGa                          <  1:8577038/150‑1 (MQ=255)
                                                                                                                          aaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGc                        >  1:4760632/1‑150 (MQ=255)
                                                                                                                            aaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCgg                      >  1:9210435/1‑150 (MQ=255)
                                                                                                                            aaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCgg                      >  2:131413/1‑150 (MQ=255)
                                                                                                                             aGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCggg                     <  1:1036365/150‑1 (MQ=255)
                                                                                                                             aGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCggg                     >  1:9095700/1‑150 (MQ=255)
                                                                                                                             aGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCggg                     >  2:2769376/1‑150 (MQ=255)
                                                                                                                              ggCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGa                    >  1:366636/1‑150 (MQ=255)
                                                                                                                               gccgccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCCGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAg                   >  1:7464763/1‑150 (MQ=255)
                                                                                                                                ccgccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGc                  <  2:7893406/150‑1 (MQ=255)
                                                                                                                                    ccGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTgg              >  1:2165476/1‑150 (MQ=255)
                                                                                                                                     cGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGc             <  1:8029820/150‑1 (MQ=255)
                                                                                                                                     cGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGc             <  1:9663228/150‑1 (MQ=255)
                                                                                                                                        aGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGaa          >  2:3698291/1‑150 (MQ=255)
                                                                                                                                         ggcggcCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAg         >  1:10458530/1‑150 (MQ=255)
                                                                                                                                          gcggcCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGa        <  1:949103/150‑1 (MQ=255)
                                                                                                                                            ggcCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACg      >  2:2635097/1‑150 (MQ=255)
                                                                                                                                               cTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTcc   <  1:7317826/150‑1 (MQ=255)
                                                                                                                                               cTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTcc   >  1:1582041/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: