Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	359,089	T→C	5.4%	G97G (GGT→GGC)	cynT →	carbonic anhydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	359,089	0	T	C	5.4%	215.7 / 2.5	75	G97G (GGT→GGC)	cynT	carbonic anhydrase
Reads supporting (aligned to +/- strand):  ref base T (28/41);  new base C (2/2);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGATTTACCGTCAAAAGCTGCG  >  NC_000913/358941‑359236                                                                                                                                                     |                                                                                                                                                    tGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGt                                                                                                                                                     >  7:357740/1‑149 (MQ=255)         cAACGGGAGCCTGGCGATCGGGTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCaa                                                                                                                                             <  3:59694/149‑1 (MQ=255)          aaCGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAAc                                                                                                                                            <  8:357740/149‑1 (MQ=255)          aaCGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGAGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAAc                                                                                                                                            >  2:158465/1‑149 (MQ=255)              tGAGCCTGGCGATCTGTTCGTTATTCGCAACGGGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTTTCTGACATTGTGATTTGTGGTCATTCCAACtgtg                                                                                                                                        <  6:368366/149‑1 (MQ=255)              tGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACtgtg                                                                                                                                        <  5:193744/149‑1 (MQ=255)              tGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGGCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACtgtg                                                                                                                                        <  2:207105/149‑1 (MQ=255)                aGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGgc                                                                                                                                      <  7:347735/149‑1 (MQ=255)                 gCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGGGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGgcg                                                                                                                                     <  2:229944/149‑1 (MQ=255)                    tGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGa                                                                                                                                  >  3:49345/1‑149 (MQ=255)                     ggCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGAt                                                                                                                                 >  3:183130/1‑149 (MQ=255)                      gCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGAt                                                                                                                                 <  1:158125/148‑1 (MQ=255)                      gCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATg                                                                                                                                >  4:1920/1‑149 (MQ=255)                            tGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGc                                                                                                                           >  8:369388/1‑148 (MQ=255)                             gTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCa                                                                                                                         >  6:279762/1‑149 (MQ=255)                              ttCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGAc                                                                                                                              <  7:21700/143‑1 (MQ=255)                              ttCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTGCGGTGGAGTATGCCGGCGATGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGAcc                                                                                                                             <  6:29661/144‑1 (MQ=255)                                cGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCAtt                                                                                                                       <  2:123050/148‑1 (MQ=255)                                 gTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGACCCCGGGGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCGGCGCGGCGGGGATCTGACATTGTGATTTGTGGCCAGTCCACCTGTGGCGCGATGACCGCCCTTGc                                                                                                                     >  4:157498/1‑149 (MQ=255)                                 gTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGCCATTCCAACTGTGGCGCGATGACCGCCATTGc                                                                                                                     >  8:50528/1‑149 (MQ=255)                                  ttATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGGGGAGTAGGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGcc                                                                                                                    <  5:205379/149‑1 (MQ=255)                                     ttCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGc                                                                                                                 >  1:327044/1‑149 (MQ=255)                                       cGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCt                                                                                                                <  7:357608/148‑1 (MQ=255)                                       cGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTg                                                                                                               >  5:385496/1‑149 (MQ=255)                                        gCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCAtt                                                                                                                                                 <  5:84987/114‑1 (MQ=255)                                        gCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCAtt                                                                                                                                                 >  6:84987/1‑114 (MQ=255)                                        gCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGt                                                                                                              >  2:170782/1‑149 (MQ=255)                                        gCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGt                                                                                                              <  3:83469/149‑1 (MQ=255)                                          aaCGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCa                                                                                                            <  7:259450/149‑1 (MQ=255)                                            cgcgGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGt                                                                                                          <  1:170782/149‑1 (MQ=255)                                            cgcgGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGt                                                                                                          <  5:254466/149‑1 (MQ=255)                                                ggCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCAt                                                                                                      <  2:398154/149‑1 (MQ=255)                                                     tatCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGCCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGAcc                                                                                                 <  7:50528/149‑1 (MQ=255)                                                     tatCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTGCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTCTCTGACATTGTGATTTGTGGCCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGAcc                                                                                                 <  3:157498/149‑1 (MQ=255)                                                          tCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATAtg                                                                                            >  8:333915/1‑149 (MQ=255)                                                            ccTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTCGAGTATGCCGTCGCTGCTCTTCGGGTATATGACATTGTGAATTGTGGTCAATCCAACTTAGGCGCGATGACAGCCATTGCCAGCTGTCAGTGCATGGACCATAtgc                                                                                           >  2:86889/1‑148 (MQ=255)                                                             cTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATAtgcct                                                                                         >  7:207627/1‑149 (MQ=255)                                                                 cTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCCGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCg                                                                                     >  2:138574/1‑149 (MQ=255)                                                                      ggCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTcc                                                                                <  1:53857/149‑1 (MQ=255)                                                                       gCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTccc                                                                               <  5:279762/149‑1 (MQ=255)                                                                       gCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTcc                                                                                >  2:337256/1‑148 (MQ=255)                                                                        ccGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCa                                                                              >  3:115597/1‑149 (MQ=255)                                                                        ccGGAACCCGGAGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCa                                                                              <  1:158465/149‑1 (MQ=255)                                                                         cGGAACCGGGTGTCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCTGTTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCTCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCAc                                                                             <  5:68712/149‑1 (MQ=255)                                                                          ggAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCa                                                                              <  6:26613/147‑1 (MQ=255)                                                                          ggAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCa                                                                              >  5:26613/1‑147 (MQ=255)                                                                          ggAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACt                                                                            <  2:386727/149‑1 (MQ=255)                                                                           gAACCCGGTGGCGTTTCTGCTTCTGTGGAGTATCCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGGCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTg                                                                           <  3:85271/149‑1 (MQ=255)                                                                           gAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTg                                                                           <  2:327044/149‑1 (MQ=255)                                                                            aaCCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCAATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTgg                                                                          <  7:26628/149‑1 (MQ=255)                                                                              cccGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCt                                                                        <  4:183130/149‑1 (MQ=255)                                                                               ccGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTg                                                                       <  6:385496/149‑1 (MQ=255)                                                                                   tGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGtt                                                                   <  7:333915/149‑1 (MQ=255)                                                                                               ttCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGCACATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAg                                                       >  1:297249/1‑149 (MQ=255)                                                                                                 cGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAg                                                       <  7:369388/147‑1 (MQ=255)                                                                                                         tATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGcgtcgg                                             >  6:194499/1‑148 (MQ=255)                                                                                                           tGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGtt                                            <  5:21271/148‑1 (MQ=255)                                                                                                             ccGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAAt                                         >  6:168678/1‑149 (MQ=255)                                                                                                             ccGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAAt                                         <  4:115597/149‑1 (MQ=255)                                                                                                                  gCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGgc                                    <  1:57176/149‑1 (MQ=255)                                                                                                                   ctccgagtCGGGTAGCTGACATTGTGATTTGTGGTCATTCCAACTGCGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGCCCATATGCCTGCCGTCCCCCCCTCGCTGCGTTATGCCGATTCAGCCCGCGTCGTAAATGAGgcg                                   <  4:49345/142‑1 (MQ=255)                                                                                                                    tGCGCTTCGGGTATCTTACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGgcgc                                  <  4:6231/149‑1 (MQ=255)                                                                                                                    tGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTc                                                         >  7:231958/1‑126 (MQ=255)                                                                                                                    tGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTc                                                         <  8:231958/126‑1 (MQ=255)                                                                                                                      cgcTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATTAGgcgcgc                                >  4:217917/1‑149 (MQ=255)                                                                                                                       gcTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGcc                               >  5:192749/1‑149 (MQ=255)                                                                                                                        cTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGAGccc                              >  3:66584/1‑149 (MQ=255)                                                                                                                               tATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTcc                       <  1:30102/149‑1 (MQ=255)                                                                                                                               tATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTcc                       >  7:85475/1‑149 (MQ=255)                                                                                                                                 tCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGa                     >  6:45158/1‑149 (MQ=255)                                                                                                                                    gACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGAtt                   >  2:321808/1‑148 (MQ=255)                                                                                                                                     aCATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGAt                    <  8:207627/146‑1 (MQ=255)                                                                                                                                     aCATTGTGAGTTGTGGTCATTCCATCTGTGTCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGAt                    <  2:297249/146‑1 (MQ=255)                                                                                                                                          gtgATTTGTGGGCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGATTTACCGt             <  1:377426/148‑1 (MQ=255)                                                                                                                                           tgATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGATTTACCGTc            <  1:138574/148‑1 (MQ=255)                                                                                                                                                    gTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGACGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGATATACCGTCAAAAGCTgcg  >  1:336397/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTTCTGCTTCGGTGGAGTATGCCGTCGCTGCGCTTCGGGTATCTGACATTGTGATTTGTGGTCATTCCAACTGTGGCGCGATGACCGCCATTGCCAGCTGTCAGTGCATGGACCATATGCCTGCCGTCTCCCACTGGCTGCGTTATGCCGATTCAGCCCGCGTCGTTAATGAGGCGCGCCCGCATTCCGATTTACCGTCAAAAGCTGCG  >  NC_000913/358941‑359236

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