Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	837,162	G→A	5.9%	C271Y (TGC→TAC)	ybiC →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	837,162	0	G	A	5.9%	211.7 / 3.0	68	C271Y (TGC→TAC)	ybiC	putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (35/33)
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.11e-01

GCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTGGGAAGTGAACACGCG  >  NC_000913/837019‑837309                                                                                                                                                |                                                                                                                                                    gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGa                                                                                                                                                >  3:310188/1‑149 (MQ=255)  ccGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGAc                                                                                                                                               <  4:310188/149‑1 (MQ=255)   cGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGAc                                                                                                                                               >  8:5774/1‑148 (MQ=255)         gTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTat                                                                                                                                         >  7:348883/1‑148 (MQ=255)             gCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatca                                                                                                                                    >  5:255239/1‑149 (MQ=255)             gCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatc                                                                                                                                     <  7:303597/148‑1 (MQ=255)              cTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                                                                                                                   <  7:191677/149‑1 (MQ=255)               tGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTTCATGACCACTatcatcat                                                                                                                                   >  2:447851/1‑148 (MQ=255)                   taCCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGCGTCAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAAcc                                                                                                                              <  2:392681/148‑1 (MQ=255)                   gACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAACCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAAc                                                                                                                               <  7:299338/148‑1 (MQ=255)                      cTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCgg                                                                                                                           >  1:303109/1‑149 (MQ=255)                      cTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCgg                                                                                                                           <  7:5774/149‑1 (MQ=255)                       tttGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatca                                                                                                                                    <  2:234584/139‑1 (MQ=255)                       tttGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatca                                                                                                                                    >  1:234584/1‑139 (MQ=255)                         tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGaa                                                                                                                         <  6:255239/148‑1 (MQ=255)                            cGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAActct                                                                                                                     >  3:398624/1‑149 (MQ=255)                             gAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCtt                                                                                                                    >  1:268234/1‑149 (MQ=255)                               aCATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCg                                                                                                                  <  8:348883/149‑1 (MQ=255)                               aCATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCg                                                                                                                  <  4:407145/149‑1 (MQ=255)                                 aTAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGgc                                                                                                                >  2:18400/1‑149 (MQ=255)                                  tAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGgcg                                                                                                               >  6:454519/1‑149 (MQ=255)                                     aGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTACATGACCACTATCATCATCAACCCGGAACTCTTCGgcgcgc                                                                                                            >  2:332093/1‑149 (MQ=255)                                       gCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCg                                                                                                          <  2:303109/149‑1 (MQ=255)                                            gCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTg                                                                                                     <  7:292853/149‑1 (MQ=255)                                              ccTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                                                                                                    >  1:35685/1‑148 (MQ=255)                                                  gCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAc                                                                                                 <  2:153931/147‑1 (MQ=255)                                                   cAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg                                                                                              <  7:71665/149‑1 (MQ=255)                                                   cAGCGATGGGTGAAATTCTTGGCGCGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgc                                                                                               <  7:426484/148‑1 (MQ=255)                                                    aGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg                                                                                              >  1:225745/1‑148 (MQ=255)                                                      cGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg                                                                                              >  7:147492/1‑146 (MQ=255)                                                      cGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg                                                                                              >  7:396422/1‑146 (MQ=255)                                                             tgAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatc                                                                                                                                     <  7:43880/100‑1 (MQ=255)                                                             tgAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatc                                                                                                                                     >  8:43880/1‑100 (MQ=255)                                                             tgAAATTCTTGGCGGGCCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTACAGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCTGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAg                                                                                    >  2:276872/1‑149 (MQ=255)                                                               aaaTTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGcc                                                                                  >  1:17774/1‑149 (MQ=255)                                                                  ttCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCAGCAGGAAACGTTACAAACCAGTCCCGAGGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCttt                                                                               <  6:377277/149‑1 (MQ=255)                                                                   tCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGCCGCGCCTGAATGTAACGCGCAGACCGAAGCCTTTg                                                                              >  8:211408/1‑149 (MQ=255)                                                                     ttGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGcc                                                                            >  6:134359/1‑149 (MQ=255)                                                                         cGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGt                                                                        <  2:35685/149‑1 (MQ=255)                                                                             gcgcTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGt                                                                    >  5:451025/1‑149 (MQ=255)                                                                                   tCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGaaa                                                                <  2:268234/147‑1 (MQ=255)                                                                                   tCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTACATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGaaa                                                                <  1:332093/147‑1 (MQ=255)                                                                                      ggcggTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTc                                                           <  6:451025/149‑1 (MQ=255)                                                                                       gcggTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGaa                                                                                     >  5:386416/1‑122 (MQ=255)                                                                                       gcggTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGaa                                                                                     <  6:386416/122‑1 (MQ=255)                                                                                             aaaCGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCat                                                    <  1:18400/149‑1 (MQ=255)                                                                                             aaaCGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCat                                                    <  1:447851/149‑1 (MQ=255)                                                                                              aaCGACGCATCAGGTAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGa                                                                          >  4:252479/1‑126 (MQ=255)                                                                                              aaCGACGCATCAGGTAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGa                                                                          <  3:252479/126‑1 (MQ=255)                                                                                                  acgCATCAGGAAACGTTACATACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCatgatga                                               <  4:398624/149‑1 (MQ=255)                                                                                                  acgCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCatgatga                                               <  5:134359/149‑1 (MQ=255)                                                                                                       tCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATgatgatgat                                           >  6:456304/1‑148 (MQ=255)                                                                                                          ggAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGc                                       <  5:454519/149‑1 (MQ=255)                                                                                                                 ttACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGAtttt                                <  2:17774/149‑1 (MQ=255)                                                                                                                   aCAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGc                              >  7:267586/1‑149 (MQ=255)                                                                                                                     aaaCCAGTCCCGATGCCATTCTTAACTACATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTa                            >  3:240062/1‑149 (MQ=255)                                                                                                                      aaCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTAc                           <  8:147492/149‑1 (MQ=255)                                                                                                                      aaCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTAc                           <  8:396422/149‑1 (MQ=255)                                                                                                                       aCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTAcc                          >  7:345926/1‑149 (MQ=255)                                                                                                                                   gCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTggg               >  1:420670/1‑148 (MQ=255)                                                                                                                                   gCCATTCTTAACTGCATGACAACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTggg               >  5:76295/1‑148 (MQ=255)                                                                                                                                    ccATTCTTAACTACATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTGGGaa             <  4:240062/149‑1 (MQ=255)                                                                                                                                     cATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTGGGAAg            >  3:39835/1‑149 (MQ=255)                                                                                                                                       ttCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGc                      <  8:76248/137‑1 (MQ=255)                                                                                                                                       ttCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGc                      >  7:76248/1‑137 (MQ=255)                                                                                                                                          ttAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTGGGAAGTGaa        >  7:467382/1‑148 (MQ=255)                                                                                                                                           tAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTGCCGGGCGAGTGGGAAGTGAaca      >  3:259777/1‑149 (MQ=255)                                                                                                                                               tGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTGGGAAGTGAACAcgcg  >  4:408652/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCGAGTGGGTGAAAGCCTCGCCGCATGATGATGATAAGCCGATTTTGCTACCGGGCGAGTGGGAAGTGAACACGCG  >  NC_000913/837019‑837309

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: ‑