Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,664,144	A→G	5.1%	N235D (AAC→GAC)	suhB →	inositol monophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,664,144	0	A	G	5.1%	246.2 / 2.5	78	N235D (AAC→GAC)	suhB	inositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base G (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTTCTGCGGCGCTGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAGGCGC  >  NC_000913/2663998‑2664292                                                                                                                                                   |                                                                                                                                                     gTTCTGCGGCGCTGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAAc                                                                                                                                                    >  2:122010/1‑149 (MQ=255)  ttCTGCGGCGCTGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACt                                                                                                                                                   <  4:113310/149‑1 (MQ=255)       cggcgCTGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATg                                                                                                                                              <  3:83149/149‑1 (MQ=255)       cggcgCGGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTCTCTTTGAAATCGGTCTGCGCCCGTGGGACTGCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTAc                                                                                                                                                 <  4:184286/146‑1 (MQ=255)             tGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGAcc                                                                                                                                        >  1:381602/1‑149 (MQ=255)                     gCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGt                                                                                                                                     <  2:381602/144‑1 (MQ=255)                     gCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTaa                                                                                                                                   >  6:110275/1‑146 (MQ=255)                     gCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTaa                                                                                                                                   <  5:110275/146‑1 (MQ=255)                          cGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGtt                                                                                                                            <  7:396940/148‑1 (MQ=255)                          cGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGtt                                                                                                                            <  4:197689/148‑1 (MQ=255)                                   gggTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATg                                                                                                                                              <  7:38876/121‑1 (MQ=255)                                   gggTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATg                                                                                                                                              >  8:38876/1‑121 (MQ=255)                                             gACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCgttgt                                                                                                        >  8:149868/1‑149 (MQ=255)                                                ggTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTa                                                                                                      >  5:313626/1‑148 (MQ=255)                                                 gTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGGGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTaaa                                                                                                    <  6:232054/149‑1 (MQ=255)                                                 gTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTaaa                                                                                                    <  8:318447/149‑1 (MQ=255)                                                      tttGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATGACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCAt                                                                                               >  7:242112/1‑149 (MQ=255)                                                      tttGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATGACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCAt                                                                                               >  7:247885/1‑149 (MQ=255)                                                      tttGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCAt                                                                                               >  3:368506/1‑149 (MQ=255)                                                      tttGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCAt                                                                                               >  7:248604/1‑149 (MQ=255)                                                       ttGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACCTCGTTGCTGGTAACCCCCGCGATGGTAAAGCCATg                                                                                              >  6:180569/1‑149 (MQ=255)                                                       ttGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATg                                                                                              >  4:170943/1‑149 (MQ=255)                                                        tGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGc                                                                                             >  1:33915/1‑149 (MQ=255)                                                         gAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCt                                                                                            >  5:214189/1‑149 (MQ=255)                                                          aaaTCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTg                                                                                           <  1:340475/149‑1 (MQ=255)                                                                gTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGaa                                                                                      >  4:131266/1‑148 (MQ=255)                                                                   tGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATg                                                                                  >  2:269643/1‑149 (MQ=255)                                                                     cgcCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCg                                                                                <  3:170943/149‑1 (MQ=255)                                                                      gcCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGt                                                                               <  5:107808/149‑1 (MQ=255)                                                                        ccGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGa                                                                             >  6:40312/1‑149 (MQ=255)                                                                         cGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGAc                                                                            <  6:313626/149‑1 (MQ=255)                                                                           tGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGa                                                                          <  1:122010/149‑1 (MQ=255)                                                                           tGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGa                                                                          <  1:269643/149‑1 (MQ=255)                                                                             ggACTTCGCCTCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGt                                                                        <  7:111449/149‑1 (MQ=255)                                                                              gACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGtt                                                                       <  4:84979/149‑1 (MQ=255)                                                                              gACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGtt                                                                       <  8:248604/149‑1 (MQ=255)                                                                               aCTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACAt                                                                                                                                >  8:201093/1‑91 (MQ=255)                                                                               aCTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACAt                                                                                                                                <  7:201093/91‑1 (MQ=255)                                                                                   cgccgcAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGATGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGc                                                                   >  4:80380/1‑148 (MQ=255)                                                                                        cAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGct                                                             >  1:46348/1‑149 (MQ=255)                                                                                          ggCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGctct                                                           >  2:305333/1‑149 (MQ=255)                                                                                            cgcgctgccGGTTCCGGAACGGGGCGGCATCGTCACCCTCCTCCCCGGGTGTCATAACTACATGCTGACCGGTAACATCGTTGCTTGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGctct                                                           <  2:68447/144‑1 (MQ=255)                                                                                               gctgctGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATGACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAATCGACGCTCTGAAGc                                                      <  8:242112/149‑1 (MQ=255)                                                                                               gctgctGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATGACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGc                                                      <  8:247885/149‑1 (MQ=255)                                                                                               gctgctGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGCCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGATGc                                                      <  1:214910/149‑1 (MQ=255)                                                                                                 tgctgGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAATGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGt                                                    >  4:247426/1‑149 (MQ=255)                                                                                                  gctgGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGa                                                                          >  8:287981/1‑126 (MQ=255)                                                                                                  gctgGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGa                                                                          <  7:287981/126‑1 (MQ=255)                                                                                                   ctgGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGtt                                                   >  7:43561/1‑148 (MQ=255)                                                                                                    tgGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTaa                                                 <  5:180569/149‑1 (MQ=255)                                                                                                     ggTTCGTGAAGGGGGCGGCATCGTCTGCGACTTCACCGGTGGTCATAACTACCTGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATTCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAAt                                                <  4:21398/149‑1 (MQ=255)                                                                                                      gTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACAGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTaa                                                 <  1:81441/147‑1 (MQ=255)                                                                                                       ttCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGa                                              <  7:149868/149‑1 (MQ=255)                                                                                                          gTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTc                                           >  1:137515/1‑149 (MQ=255)                                                                                                           tGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACt                                            >  6:338575/1‑147 (MQ=255)                                                                                                                cgggcggCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCgg                                     <  8:43561/149‑1 (MQ=255)                                                                                                                 gggcggCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCggg                                    <  3:340112/149‑1 (MQ=255)                                                                                                                 gggcggCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCggg                                    <  4:368506/149‑1 (MQ=255)                                                                                                                  ggcggcATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGt                                   <  1:229836/149‑1 (MQ=255)                                                                                                                    cggcATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGa                                 <  6:214189/149‑1 (MQ=255)                                                                                                                     ggcATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGa                                 >  4:130700/1‑148 (MQ=255)                                                                                                                       cATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAg                                         >  6:24072/1‑138 (MQ=255)                                                                                                                       cATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAg                                         <  5:24072/138‑1 (MQ=255)                                                                                                                       cATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGatat                              >  2:199275/1‑149 (MQ=255)                                                                                                                            tCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTa                                                  <  4:373938/124‑1 (MQ=255)                                                                                                                            tCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTa                                                  >  3:373938/1‑124 (MQ=255)                                                                                                                            tCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATAtcact                          >  3:119134/1‑148 (MQ=255)                                                                                                                            tCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTGAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATAtcact                          >  1:274172/1‑148 (MQ=255)                                                                                                                              aGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATTCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATAtcactca                        <  3:131266/148‑1 (MQ=255)                                                                                                                              aGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTACGCGACGCTCTGAAGCGGTAATGACTCAGGCGGGTGATAtcactcac                       >  3:196463/1‑149 (MQ=255)                                                                                                                              aGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATAtcactcac                       >  3:88708/1‑149 (MQ=255)                                                                                                                               gCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAAAATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACTCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCAcc                      >  2:149926/1‑149 (MQ=255)                                                                                                                                   cTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGcc                  <  4:299189/149‑1 (MQ=255)                                                                                                                                      cACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTc                                           <  7:42032/121‑1 (MQ=255)                                                                                                                                      cACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTc                                           >  8:42032/1‑121 (MQ=255)                                                                                                                                               tCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAATCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAg      >  3:347066/1‑149 (MQ=255)                                                                                                                                                 aaaaCTACATGCTGACCGGTAACATCGTTGCTGGTACCCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAGgc    >  8:391421/3‑149 (MQ=255)                                                                                                                                                 aTAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAgg     <  4:88708/148‑1 (MQ=255)                                                                                                                                                 aTAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAGgc    <  3:247426/149‑1 (MQ=255)                                                                                                                                                   aaCTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGAGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAGgcgc  >  1:367452/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTTCTGCGGCGCTGGATCTGGCTTACGTCGCTGCGGGTCGTGTTGACGGTTTCTTTGAAATCGGTCTGCGCCCGTGGGACTTCGCCGCAGGCGAGCTGCTGGTTCGTGAAGCGGGCGGCATCGTCAGCGACTTCACCGGTGGTCATAACTACATGCTGACCGGTAACATCGTTGCTGGTAACCCGCGCGTTGTTAAAGCCATGCTGGCGAACATGCGTGACGAGTTAAGCGACGCTCTGAAGCGTTAATGACTCAGGCGGGTGATATCACTCACCCGCCCTCGCCTTTCAGGCGC  >  NC_000913/2663998‑2664292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑