Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 6.2% intergenic (+70/+17) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,5150TA 6.2% 656.3 / 21.6 196intergenic (+70/+17)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (105/75);  new base A (2/10);  total (109/87)
Fisher's exact test for biased strand distribution p-value = 6.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 202 total aligned reads displayed.

ACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATG  >  NC_000913/3959268‑3959789
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
aCCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGCAAACACGCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                                                                                                                                                          >  2:320877/1‑275 (MQ=255)
                   cAACCGGGCGACCTGGGTAAAGCTCTTCCGGAAGGCGCGGTAGATACCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCCGGTAGGTAAGCAACTGCGCGGCTACATGACAGATATGCAACGTATTGCTGTTGCGGGGTGAGTTCGCGCTGATGCCCTCAACCCGCCCCTCTCCCACAGGGACCGGGAGACAACCCCGCACTCCTTCTCCTGGGCCAGGCCTTGCTCTTACTTCAGGTTCCACACCTTAatgtgttt                                                                                                                                                                                                                                       >  4:239956/1‑270 (MQ=255)
                        gggCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGGGAAAACCCCCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGGGGTACGGCCCCTTAATGATGTGatat                                                                                                                                                                                                                                   >  4:516512/1‑274 (MQ=255)
                            gACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                                                                                                                                                                      <  2:516512/267‑1 (MQ=255)
                                                                     gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                                               >  4:83625/1‑185 (MQ=255)
                                                                                               ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                                                                                                                                                                                                                                                                >  3:90603/1‑158 (MQ=255)
                                                                                               ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGACCTCACCCCGACCCTCCCCCACAGGGTGAGGGAGAAAACACTCAAGGCCTACTCCTGGCGAAGGCCTTGCTAttt                                                                                                                                                                                                                                                                  >  1:3105/1‑171 (MQ=255)
                                                                                                   cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                                                                                                                                                                     <  8:639693/149‑1 (MQ=255)
                                                                                                   cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                                                                                                                                                                     <  7:277118/149‑1 (MQ=255)
                                                                                                   cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                                                                                                                                                            <  1:90603/158‑1 (MQ=255)
                                                                                                    aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                                                                                                                                                    <  8:756475/149‑1 (MQ=255)
                                                                                                     gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTgg                                                                                                                                                                                                                                                                                   >  5:655431/1‑149 (MQ=255)
                                                                                                         tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                                               <  8:539568/149‑1 (MQ=255)
                                                                                                          gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGGGGGAGAACACACTCCAGGCCTTCTCCAGGAGAAg                                                                                                                                                                                                                                                                              <  6:875135/149‑1 (MQ=255)
                                                                                                           cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                                                                                                                                                                                                             <  6:631507/149‑1 (MQ=255)
                                                                                                             aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                                                                                                                                                           <  7:353025/149‑1 (MQ=255)
                                                                                                             aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                                                                                                                                                           >  7:227452/1‑149 (MQ=255)
                                                                                                             aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                                                                                                                                                                                           <  7:702124/149‑1 (MQ=255)
                                                                                                             aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                                                                                                                                                                                           <  8:620911/149‑1 (MQ=255)
                                                                                                                gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                                               >  1:172343/1‑142 (MQ=255)
                                                                                                                gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGAGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCGCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                                                                                                                                                        <  6:114628/149‑1 (MQ=255)
                                                                                                                 aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCTCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTga                                                                                                                                                                                                                                                                       <  8:361353/149‑2 (MQ=255)
                                                                                                                 aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTga                                                                                                                                                                                                                                                                       <  5:600905/149‑2 (MQ=255)
                                                                                                                 aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTga                                                                                                                                                                                                                                                                       <  7:816404/149‑2 (MQ=255)
                                                                                                                  gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                                                                                                                                                                                      <  7:334286/149‑1 (MQ=255)
                                                                                                                     ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                                                                                                                                                                                           >  6:729319/1‑141 (MQ=255)
                                                                                                                      gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                                                                                                                                                                  <  5:938131/149‑1 (MQ=255)
                                                                                                                       taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTAGCTATTAg                                                                                                                                                                                                                                                                 >  6:949660/1‑149 (MQ=255)
                                                                                                                                aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTTCGTTACAGCAACTTAATGATGTGACATCCGAACTGGGTGTGCATCGGGCCGGTTNGCCCCAGAACCGGCCAAGAGAAAACCACTTTAAACAACCCCGGACCCCACTGCACCCAGGGGACATCACCTAAATc                                                                                                                          >  3:151064/1‑274 (MQ=255)
                                                                                                                                    gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                                                                                                                                                    >  4:103532/1‑133 (MQ=255)
                                                                                                                                     cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatg                                                                                                                                                                                                                                            >  1:79961/1‑156 (MQ=255)
                                                                                                                                      gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                                                                                                                                                                  <  8:346066/149‑1 (MQ=255)
                                                                                                                                        ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                                                                                                                                                                                                                                                <  2:103532/133‑1 (MQ=255)
                                                                                                                                        ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                                                                                                                                                                                >  5:662411/1‑149 (MQ=255)
                                                                                                                                         gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                                                                                                                                                                                                                                               <  6:802323/149‑1 (MQ=255)
                                                                                                                                          cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                                                                                                                                                              <  5:458432/149‑1 (MQ=255)
                                                                                                                                           tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                                                                                                                                                                             >  5:1066084/1‑149 (MQ=255)
                                                                                                                                                gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                                                                                                                                                        >  5:758333/1‑117 (MQ=255)
                                                                                                                                                gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                                                        <  6:49081/149‑1 (MQ=255)
                                                                                                                                                 aCAGTTATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                                                                                                                                                                       <  7:99137/149‑1 (MQ=255)
                                                                                                                                                    gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                                                                                                                                                                                                                                 >  3:107971/1‑120 (MQ=255)
                                                                                                                                                     atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                                                                                                                                                                 >  2:177965/1‑151 (MQ=255)
                                                                                                                                                     atatGAAACGTATTGCCGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                                                                                                                                                                                                                                   <  5:115368/149‑1 (MQ=255)
                                                                                                                                                      tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                                                                                                                                                  >  5:225825/1‑149 (MQ=255)
                                                                                                                                                      tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                                                                                                                                                  >  8:1048687/1‑149 (MQ=255)
                                                                                                                                                      taggAACCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAAGCCCTGACCCCGACCCTGGCCGACAGGGAGAGGGAGAAAACACTCAAGGCATTCTCCTGGAGAAGGCCTTGCTATTAGTGGCGGTACAGCACCTTAACGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCGGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATggg                                                                                                    <  4:329946/271‑1 (MQ=255)
                                                                                                                                                        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                                                                                                                                                             <  1:107971/120‑1 (MQ=255)
                                                                                                                                                        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                                                        <  2:585028/141‑1 (MQ=255)
                                                                                                                                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                                                                                                                                                                                                                               >  7:208669/1‑149 (MQ=255)
                                                                                                                                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                                                                                                                                                             <  4:177965/151‑1 (MQ=255)
                                                                                                                                                            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                                                                                                                                                                                     >  8:1077627/1‑140 (MQ=255)
                                                                                                                                                               tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                                                                                                                                                  >  2:41813/1‑92 (MQ=255)
                                                                                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGTGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                                                                                                                                                        <  5:924496/149‑1 (MQ=255)
                                                                                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                                                                                                                                                        <  6:660324/149‑1 (MQ=255)
                                                                                                                                                                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                                                                                                                                                                                                                                 >  2:748023/1‑107 (MQ=255)
                                                                                                                                                                   gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                                                                                                                                                              <  4:41813/92‑1 (MQ=255)
                                                                                                                                                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                                                                                                                                                          >  4:594238/1‑106 (MQ=255)
                                                                                                                                                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                                                                                                                                                          >  3:176773/1‑106 (MQ=255)
                                                                                                                                                                              ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGc                                                                            >  2:239436/1‑274 (MQ=255)
                                                                                                                                                                                  gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                                                                                                                                                                                      >  7:624778/1‑149 (MQ=255)
                                                                                                                                                                                   tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGc                                                                                                                                                                                                     <  7:762338/149‑1 (MQ=255)
                                                                                                                                                                                       cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                                                                                                                                                                            >  1:137160/1‑154 (MQ=255)
                                                                                                                                                                                          tGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                                                                                                                                                              <  2:611914/85‑1 (MQ=255)
                                                                                                                                                                                                 cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                                                                                                                                                                       >  5:537622/1‑149 (MQ=255)
                                                                                                                                                                                                 cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                                                                                                                                                                       >  5:309297/1‑149 (MQ=255)
                                                                                                                                                                                                  tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                                                                                                                                                      >  8:460098/1‑149 (MQ=255)
                                                                                                                                                                                                    aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                                                                                                                                                    >  6:84282/1‑149 (MQ=255)
                                                                                                                                                                                                    aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa                                                                                                                                                                  >  4:705346/1‑167 (MQ=255)
                                                                                                                                                                                                      cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                                                                                                                                                                                     >  4:539741/1‑146 (MQ=255)
                                                                                                                                                                                                      cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                                                                                                                                                  >  6:537639/1‑149 (MQ=255)
                                                                                                                                                                                                      cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                                                                                                                                                  >  5:86178/1‑149 (MQ=255)
                                                                                                                                                                                                       ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                                                                                                                                                 >  7:9035/1‑149 (MQ=255)
                                                                                                                                                                                                       ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGCGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                                                                                                                                                 >  7:831503/1‑149 (MQ=255)
                                                                                                                                                                                                        cGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                                                                                                                                                              <  2:705346/167‑1 (MQ=255)
                                                                                                                                                                                                         gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                                                                                                                                                               <  7:460098/149‑1 (MQ=255)
                                                                                                                                                                                                          aCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                                                                                                                                                 <  2:539741/146‑1 (MQ=255)
                                                                                                                                                                                                          aCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                                                                                                                                                              <  8:772439/149‑1 (MQ=255)
                                                                                                                                                                                                             ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                                                                                                                                                                           <  7:85082/149‑1 (MQ=255)
                                                                                                                                                                                                             ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                                                                                                                                                                           >  6:217084/1‑149 (MQ=255)
                                                                                                                                                                                                               ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                                                                                                                                                                         >  7:955174/1‑149 (MQ=255)
                                                                                                                                                                                                                tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                                                                                                                                                               >  2:688467/1‑158 (MQ=255)
                                                                                                                                                                                                                 cccACAGGTGAAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGAACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCgaga                                        <  2:223975/275‑1 (MQ=255)
                                                                                                                                                                                                                  ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                                                                                                                                                                      >  8:223464/1‑149 (MQ=255)
                                                                                                                                                                                                                   cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                                                                                                                                                     <  8:196677/149‑1 (MQ=255)
                                                                                                                                                                                                                    acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGa                                                                                                                                                           <  4:688467/158‑1 (MQ=255)
                                                                                                                                                                                                                          agagggagaAAACGCGCAAGGCCTTCTCCGGGAGAAGGCCTGGATATAAGTGGGGGTACAGAACATTAATGATGTGATATCCGAACTGGGTGGGCAGCGGGCCGGTCGGCTCCAGTACCGGAAAAGAGAAAACCACTTTATCGAACGCCGGAACAATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGccc                                                            <  1:118283/246‑1 (MQ=255)
                                                                                                                                                                                                                          agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAAACCACTTTATCGAACGCCGGCACCATCTGACCCTGGCGGAATTCACCTAAATCCCCGCCGCGTTTGCCTCATGGGCAAATGGAGTTTTTCTTCGCCAGCTTTCCGACATCGGCCCCCTTCTTCATCTTCTCGAAAACATccc                                 >  4:607072/1‑272 (MQ=255)
                                                                                                                                                                                                                            agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGc                             >  2:266256/1‑275 (MQ=255)
                                                                                                                                                                                                                                        cTCAAGGCCTTTGCCGGGGGAAGGGCTGGCTTTTAGTGGCGGTAAAGAACCTTAATGATGGTATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGACCGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCtctt                  <  4:239436/274‑1 (MQ=255)
                                                                                                                                                                                                                                         tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                                                                                                                                                                  >  2:203668/1‑98 (MQ=255)
                                                                                                                                                                                                                                         tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                                                                                                                                              >  7:301281/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                                                                                                                                              >  7:191480/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                                                                                                                                              >  6:1071309/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                                                                                                                                              >  8:325780/1‑149 (MQ=255)
                                                                                                                                                                                                                                             ggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                                                                                                                                                                                              <  4:203668/98‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa                                                                                                                                      <  8:882998/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ctcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCGTCTGACCNTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCNCCGTTCTTAATCTGCCCGGGAAGATNCAGAGCCATGTTCTCTTCTTTTACaatct      >  3:125631/1‑272 (MQ=255)
                                                                                                                                                                                                                                                    tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTc                                                                                                                                   >  6:1031224/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ccTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCa                                                                                                                                  <  5:834881/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       tGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAcc                                                                                                                                >  8:409001/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
ACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATG  >  NC_000913/3959268‑3959789

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