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% 1833.6 / 5.8 376intergenic (‑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 (3/0);  new base . (181/192);  total (184/192)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 408 total aligned reads displayed.

AAATGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCGG  >  NC_008463/3909483‑3909779
                                                                                                                                                      |                                                                                                                                                  
aaaTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTg                                                                                                                                                    >  2:820768/1‑151 (MQ=255)
aaaTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTg                                                                                                                                                    <  2:5885955/151‑1 (MQ=255)
 aaTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTgg                                                                                                                                                   <  1:10492260/151‑1 (MQ=255)
  aTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGACTCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCGTCCTggg                                                                                                                                                  >  2:513063/1‑151 (MQ=255)
  aTGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTggg                                                                                                                                                  >  2:10378767/1‑151 (MQ=255)
   tGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTgggg                                                                                                                                                 <  1:2687336/151‑1 (MQ=255)
       gAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTggggttcg                                                                                                                                             >  1:659761/1‑147 (MQ=255)
       gAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTggggttcg                                                                                                                                             >  2:10775373/1‑147 (MQ=255)
        aGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGc                                                                                                                                           >  2:1275615/1‑151 (MQ=255)
         gTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGcc                                                                                                                                          >  1:5246732/1‑151 (MQ=255)
          ttCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCt                                                                                                                                         >  2:9387784/1‑151 (MQ=255)
           tCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTg                                                                                                                                        <  1:8202449/151‑1 (MQ=255)
             ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGAc                                                                                                                                      >  1:5653547/1‑151 (MQ=255)
             ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGAc                                                                                                                                      <  1:6967221/151‑1 (MQ=255)
             ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGAc                                                                                                                                      <  1:7701057/151‑1 (MQ=255)
             ggcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGAc                                                                                                                                      >  2:6784897/1‑151 (MQ=255)
              gcggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACg                                                                                                                                     <  1:618864/151‑1 (MQ=255)
               cggTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGt                                                                                                                                    >  1:4165256/1‑151 (MQ=255)
                 gTCAGGACAACTCCCTCTTGAGCCGGCCGGGACTCGCCCCGCCGTTGTTGCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGtct                                                                                                                                  <  2:2731259/151‑1 (MQ=255)
                     tGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAgg                                                                                                                              >  2:9925443/1‑151 (MQ=255)
                       aaaaCTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCg                                                                                                                            <  1:1450473/151‑1 (MQ=255)
                        aaaCTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGa                                                                                                                           >  2:1449527/1‑151 (MQ=255)
                        aaaCTCCCCCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGa                                                                                                                           <  1:117138/151‑1 (MQ=255)
                          aCTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGAcc                                                                                                                         <  1:7614372/151‑1 (MQ=255)
                           cTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCg                                                                                                                        >  2:8514050/1‑151 (MQ=255)
                           cTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCg                                                                                                                        >  1:5792854/1‑151 (MQ=255)
                           cTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCg                                                                                                                        >  1:9705450/1‑151 (MQ=255)
                            tCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGc                                                                                                                       >  1:9977606/1‑151 (MQ=255)
                             cccTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCt                                                                                                                      <  1:7903094/151‑1 (MQ=255)
                              ccTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCtt                                                                                                                     >  2:4069251/1‑151 (MQ=255)
                               ctctTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTc                                                                                                                    <  1:9123092/151‑1 (MQ=255)
                                    gAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTcgcgcg                                                                                                               <  1:8981658/151‑1 (MQ=255)
                                      gccggccgGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGc                                                                                                             <  2:5388288/151‑1 (MQ=255)
                                          gccggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGc                                                                                                         <  1:979888/151‑1 (MQ=255)
                                            cggGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGccc                                                                                                       >  1:7462257/1‑151 (MQ=255)
                                              ggAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcag                                                                                                     >  1:202187/1‑151 (MQ=255)
                                               gAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagc                                                                                                    <  1:5986702/151‑1 (MQ=255)
                                               gAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagc                                                                                                    <  2:372072/151‑1 (MQ=255)
                                               gAATCGCCCCGCCGTTGTGTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagc                                                                                                    <  2:9141670/151‑1 (MQ=255)
                                                aaTCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagca                                                                                                   <  1:4933295/151‑1 (MQ=255)
                                                aaTCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagca                                                                                                   >  1:3019335/1‑151 (MQ=255)
                                                 aTCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCcagcag                                                                                                  <  2:3279138/151‑1 (MQ=255)
                                                    gCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCg                                                                                               >  2:4401951/1‑151 (MQ=255)
                                                      cccGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGAt                                                                                             >  1:469031/1‑151 (MQ=255)
                                                       ccgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATc                                                                                            <  1:6644140/151‑1 (MQ=255)
                                                        cgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCg                                                                                           <  1:8801879/151‑1 (MQ=255)
                                                        cgccgTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCg                                                                                           >  2:338160/1‑151 (MQ=255)
                                                            ggtgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCaaa                                                                                       <  1:10923085/149‑1 (MQ=255)
                                                             ttgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAAc                                                                                      <  1:7833596/151‑1 (MQ=255)
                                                              tgttTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAATAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACg                                                                                     <  2:5483413/151‑1 (MQ=255)
                                                                 ttCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACt                                                                                  <  1:1655125/151‑1 (MQ=255)
                                                                  tCGCGCCGTTGGGACGCCGCTCGCCGCGGCCTGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTg                                                                                 <  2:8078449/151‑1 (MQ=255)
                                                                    gcgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGcc                                                                               >  1:4753477/1‑151 (MQ=255)
                                                                    gcgcCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGcc                                                                               <  1:7388130/151‑1 (MQ=255)
                                                                       ccGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGCGGCGGCCTTCCTGGGGTTCGCCGGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAgg                                                                            <  2:4913046/151‑1 (MQ=255)
                                                                        cGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGa                                                                           <  2:4743220/151‑1 (MQ=255)
                                                                         gTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGAc                                                                          <  1:1183811/151‑1 (MQ=255)
                                                                          ttGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGAcc                                                                         >  2:3914077/1‑151 (MQ=255)
                                                                             ggACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCTGCAGCCGATCGCAAACGACTGCCAGGACCAgc                                                                      <  2:7253466/151‑1 (MQ=255)
                                                                              gACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCAGCACGGAGGCGTCCTTCCTGGGGTTCGGATGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGACGATCGCAAACGACTGCCTGGACCTgcg                                                                     >  2:6279198/1‑151 (MQ=255)
                                                                               aCGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAgcgc                                                                    >  1:10597928/1‑151 (MQ=255)
                                                                                   cgcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGTAAAAGGCCGCCCTGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATg                                                                <  2:29912/151‑1 (MQ=255)
                                                                                   cgcTCGCCGCGGCCAGGACACGACATTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATg                                                                <  2:7770294/151‑1 (MQ=255)
                                                                                    gcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATgg                                                               >  2:6485198/1‑151 (MQ=255)
                                                                                    gcTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATgg                                                               >  1:5493157/1‑151 (MQ=255)
                                                                                     cTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGa                                                              <  2:2368498/151‑1 (MQ=255)
                                                                                        gccgcGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATAc                                                           >  1:1506278/1‑151 (MQ=255)
                                                                                         ccgcGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATAcc                                                          >  1:6812847/1‑151 (MQ=255)
                                                                                             ggCCTGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgca                                                      <  1:2972405/151‑1 (MQ=255)
                                                                                             ggCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgca                                                      >  1:44969/1‑151 (MQ=255)
                                                                                             ggCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgca                                                      >  1:3508619/1‑151 (MQ=255)
                                                                                             ggCCAGGACACGACTTTCGTCGGCTAATGAAACAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGAGGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgca                                                      <  1:844884/151‑1 (MQ=255)
                                                                                              gCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGgcag                                                     >  2:7119003/1‑151 (MQ=255)
                                                                                                cAGGACACGACTTTCCTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCCTCCTGGGGTTCGCCTCACGTCTCAGGCGACCGCTTCGCGCGGCTGGCCCAGCAGCCGATCGCAAACGACGGCCAGGACCAGCGCAATCGATACCGgcagca                                                   <  2:2115034/151‑1 (MQ=255)
                                                                                                 aGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGcagcag                                                  >  1:9269378/1‑151 (MQ=255)
                                                                                                       cGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCACTGGATACCGGCAGCagccagcc                                            >  2:10629071/1‑151 (MQ=255)
                                                                                                         aCTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCaa                                          >  1:2238915/1‑151 (MQ=255)
                                                                                                          cTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAg                                         >  1:10251774/1‑151 (MQ=255)
                                                                                                          cTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAg                                         >  1:2747662/1‑151 (MQ=255)
                                                                                                           tttCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAgg                                        <  2:6656114/151‑1 (MQ=255)
                                                                                                              cgtcgGCTCATGAAACAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGctc                                     <  2:8742924/151‑1 (MQ=255)
                                                                                                                   gCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATc                                >  2:7716728/1‑151 (MQ=255)
                                                                                                                     tCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCgg                              <  2:10089996/151‑1 (MQ=255)
                                                                                                                      cATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGc                             <  2:562018/151‑1 (MQ=255)
                                                                                                                          aaaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAATCGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGc                         >  2:9723942/1‑151 (MQ=255)
                                                                                                                          aaaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGc                         <  1:8954649/151‑1 (MQ=255)
                                                                                                                          aaaaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGc                         >  1:7984145/1‑151 (MQ=255)
                                                                                                                            aaaGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCgg                       >  2:1652982/1‑151 (MQ=255)
                                                                                                                              aGGCCGCCCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGa                     <  2:10043811/151‑1 (MQ=255)
                                                                                                                                 ccgccCGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCt                  >  1:10667742/1‑151 (MQ=255)
                                                                                                                                    cccGGAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCCGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTgg               >  2:11141209/1‑151 (MQ=255)
                                                                                                                                       ggAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGa            >  2:2832604/1‑151 (MQ=255)
                                                                                                                                        gAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGaa           >  2:7353436/1‑151 (MQ=255)
                                                                                                                                        gAGGCGGCCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGaa           <  1:2119735/151‑1 (MQ=255)
                                                                                                                                        gAGGCGGCCTACCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGaa           <  1:9541871/151‑1 (MQ=255)
                                                                                                                                         aGGCGGCCTTCCTGGGGTTCGCCGGACGTCCCAGGCGCCCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAg          <  2:950630/151‑1 (MQ=255)
                                                                                                                                             ggcCTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGt      <  1:10001221/151‑1 (MQ=255)
                                                                                                                                               ccTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTcc    <  2:1841083/151‑1 (MQ=255)
                                                                                                                                                cTTCCTGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCg   <  1:3706325/151‑1 (MQ=255)
                                                                                                                                                cTTCCTGGGGTTCGCCTGACGTCACAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCg   <  2:5760854/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AAATGCCGAGTTCGGCGGTCATGAAAACTCCCTCTTGAGCCGGCCGGGAATCGCCCCGCCGTTGTTTCGCGCCGTTGGGACGCCGCTCGCCGCGGCCAGGACACGACTTTCGTCGGCTCATGAAAAAGGCCGCCCGGAGGCGGCCTTCCTGGGGGTTCGCCTGACGTCTCAGGCGACCGCTTCGCGCGGCTTGCCCAGCAGCCGATCGCAAACGACTGCCAGGACCAGCGCAATGGATACCGGCAGCAGCCAGCCAAGGCTCTGATCGGCGAGCGGGAGCTTGGCGAAGACGTCCGG  >  NC_008463/3909483‑3909779

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: