Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,745,494	0	C	T	39.7%	72.4 / 26.6	63	R231C (CGC→TGC)	yiaN	2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (11/27);  new base T (25/0);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 4.86e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTC  >  NC_000913/3745346‑3745644                                                                                                                                                      |                                                                                                                                                           gggCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTc                                                                                                                                                            <  2:927879/149‑1 (MQ=255)       aaCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTCGCTcc                                                                                                                                                       >  2:393507/1‑148 (MQ=255)        aCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGcc                                                                                                                                                  <  1:436021/151‑1 (MQ=255)        aCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCC‑GGcc                                                                                                                                                  >  1:108187/1‑151 (MQ=255)           cTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTg                                                                                                                                                <  1:1279425/150‑1 (MQ=255)                gCTTCCCTGGTGGTGGCGGGCCAGTCGCCTTAATTCACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑tt                                                                                                                                             <  1:760712/147‑1 (MQ=255)                  ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGGTCC‑GGGCTG‑TTTTcccc                                                                                                                                       >  2:1024501/1‑151 (MQ=255)                     ccTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCt                                                                                                                                      <  1:131378/149‑1 (MQ=255)                               gCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCggg                                                                                                                          <  1:143942/151‑1 (MQ=255)                                    ccAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTTCTCC‑GGCCCG‑TTTAtccccaacccagcccggg                                                                                                                         >  2:428525/1‑133 (MQ=255)                                    ccAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGt                                                                                                                     <  1:1150149/151‑1 (MQ=255)                                     cAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGtt                                                                                                                    <  1:161347/151‑1 (MQ=255)                                     cAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGtt                                                                                                                    <  2:1158211/151‑1 (MQ=255)                                         cGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGtt                                                                                                                    <  1:367936/147‑1 (MQ=255)                                         cGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑GGCCTG‑TTTTCCCCTCCCCAAGCCgggg                                                                                                                         >  2:764696/1‑142 (MQ=255)                                         cGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTtgc                                                                                                                  <  1:375916/149‑1 (MQ=255)                                                 tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCt                                                                                                          <  2:81387/149‑1 (MQ=255)                                                 tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCC‑GGCCTG‑TTTACCCCCACCGAAGCCGGGGGCGGTGCTGCTTTCttt                                                                                                        >  2:1102964/1‑149 (MQ=255)                                                 tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTCGCTCC‑GGCCCG‑TTTACCCCCACCCAAACCGGGGCCGTTGCTGCTTTCT‑At                                                                                                       >  1:919020/1‑151 (MQ=255)                                                  ttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTTCTCC‑GGCCTG‑TTTTCCCCCACCCAAACCGGGGCCGTTGCTGCttttttt                                                                                                        >  2:62134/1‑146 (MQ=255)                                                    aCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTTCCCCTCCCGAAGCCCGGGCCGTTGCTGCTTTCT‑TTgc                                                                                                     >  1:391887/1‑150 (MQ=255)                                                    aCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCC‑GGCCTG‑TTTTCCCCTACCCCAGCCGGGGCCCTTTCTGCTTTTT‑At                                                                                                       >  1:113879/1‑148 (MQ=255)                                                     ccccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCCT‑TTTTCCCCTTCCGCAGCCCGGGCCGGTGCTGCTTTCtttccg                                                                                                     >  1:977613/1‑145 (MQ=255)                                                       ccGCCAGCAAAAAGCAACGATGCAGGAAATCTGCCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCACCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATgcg                                                                                                    <  1:718378/148‑1 (MQ=255)                                                             gCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTc                                                                                           <  1:393335/151‑1 (MQ=255)                                                               aaaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCC‑GGCCCG‑TTTTTCCCTACCCAAACCGGGGCCGTTGCTGCTTTTT‑TTGCGCTGttt                                                                                              >  2:143478/1‑146 (MQ=255)                                                                aaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCC‑GGGCTG‑TTTTCCCCCACCGAAGCCGGGGCCCTTGCCGCTTTTT‑TTGCGCCGTTTTTCGcc                                                                                        >  1:609891/1‑151 (MQ=255)                                                                 aaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑GGCCTG‑TTTTCCCCTACCCAAGCCCGGGCCGTTGCTGCTTTCT‑TTTCGCTGTTTTTCGcc                                                                                        >  1:1278162/1‑150 (MQ=255)                                                                    gCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑CGCCCG‑TTTTCCCCTACCTAAACCGGGGGCGGTGTTGCTTTCT‑ATGCGCTGTTTTTAGCcac                                                                                      >  1:44924/1‑149 (MQ=255)                                                                    gCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCTCTCC‑CGCCTG‑TTTACCCCTACCGAAACCCGGGCCGTTGCTGCTTTTT‑ATGCGCCGTTTGTTGCcac                                                                                      >  1:916174/1‑149 (MQ=255)                                                                     cAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGt                                                                                   <  1:696995/151‑1 (MQ=255)                                                                         gATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑GGCCTG‑TTTTCCCCTACCGAAGCCCGGGCCGTTGCTGCTTTCT‑TTGCGCTGTTTGTCGCCACAGTTA‑tt                                                                              >  1:178443/1‑151 (MQ=255)                                                                              aGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCC‑GGGCCG‑TTTTGCCCTACCCGAACCCGGGCGGCTGTTGTTTTCT‑TTGCCCTCTTTTTGGCCCCAGGTA‑ttt                                                                             >  1:557486/1‑147 (MQ=255)                                                                                 aaaTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGa                                                                      <  1:1154947/151‑1 (MQ=255)                                                                                    tCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑GGCCCG‑TTTACCCCTACCCAAGCCGGGGGCGTTGCTGCTTTTT‑ATGTGCTGTTTTTCGCCCCCGTTA‑TTTTCCGTGaaa                                                                    >  2:520464/1‑150 (MQ=255)                                                                                     cTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTCCTCC‑TGCCTC‑TTTTCCCCCACCCAAGACGCGGCGGGTGCTGTTTTCTGTTGCGTCGTTT‑TCGCCCCAGCTA‑TTTTGCGTGGAAt                                                                   >  2:499810/1‑150 (MQ=255)                                                                                        gCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTTCTCC‑GGCCCG‑TTTTCCCCTACCCAAGCCGGGGCCGGTGCTGCTTTTT‑TTGCGCTGTTTTTCGCCACACTTA‑TTTTCACTGAAAAGAcc                                                               >  2:773529/1‑151 (MQ=255)                                                                                         cACGCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCTGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCt                                                              <  2:978111/151‑1 (MQ=255)                                                                                          aCTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑GGGCCG‑TTTTTCCCTACCCAAACCGGGGGCGTTGCTGCTTTCT‑TTGCGCCGTTTGTCGCCACAATTT‑TTTACCGTGAAATGAACtt                                                             >  1:1042092/1‑151 (MQ=255)                                                                                             ctTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTc                                                            <  1:644019/149‑1 (MQ=255)                                                                                             ctTTGGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTc                                                            <  1:1119284/149‑1 (MQ=255)                                                                                               tttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCg                                                           <  1:689674/148‑1 (MQ=255)                                                                                                 tGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCcac                                                       <  2:1297004/150‑1 (MQ=255)                                                                                                 tGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTCGCCCC‑GGCCCG‑TTTTACCCTACCCAAACCGGGGCCGTTGCCTCTTTTT‑ATTCGCGGTTTTTCGCCCAAGTTT‑TTTACCCGGAAAAGACCTTTGcc                                                         >  1:124208/1‑148 (MQ=255)                                                                                                  gTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑CGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCCCAGTTA‑TTTAACGTGAAATGACCTTCGCcac                                                       >  1:42909/1‑149 (MQ=255)                                                                                                   tctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTCCTCC‑GCCCCGTTTTTCCCCCACCAAACCCGGGG‑CGGTGTTGCCTTTT‑TTTCGCCGTTTTTCGCCCCACATA‑TTTTCCCGGAAAAGaacctcccccc                                                       >  2:1235997/1‑138 (MQ=255)                                                                                                      ccGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCTCCCTTTg                                                 <  2:1151665/151‑1 (MQ=255)                                                                                                       cGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑CGCCTG‑TTTACCCCCACCCAAGCCGGGGGCGGTGCTGCTTTTT‑ATTCGCTGTTTTTTGCCCCAGTTT‑TTTTCCGTGGAATGACCCTCGCCACACttttg                                                >  1:1262736/1‑149 (MQ=255)                                                                                                              tGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCCGGGCGTT‑TTTCCCCCTCCCGAAG‑CGGGGGCGTTGGCGCTTTTT‑TTTCGTCGTTTTTCGGCCCAGTTT‑TTTTACGTGGAATTACCTTCTCCCCACCTTTGCATGTAc                                         >  1:500205/1‑151 (MQ=255)                                                                                                              tGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTTC‑CGCCCG‑TTTTCCCCTCCCCAAACCCGGGGCGGTTCGGCCTTTT‑TTTCGCTGTTTTTTCCCCCAGTTA‑TTTTACGCGCAAAGACCCTCGCCCCAATTTTGCATGtgc                                         >  2:1187354/1‑149 (MQ=255)                                                                                                              tGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTCGCTCC‑GGCCTG‑TTTTCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTTTTGCCACAGTTA‑TTTTCCGTGAAATGACCTTCGCCCCACTTTGGCATGTa                                          >  1:379685/1‑150 (MQ=255)                                                                                                               gggCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTAc                                         <  1:833057/150‑1 (MQ=255)                                                                                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCCCC‑GGCCCG‑TTTTCCCCCACCCAACCCGGGGGCGGTGGTGCTTTCT‑TTGCGGCGTTTTTTCCCCCCGTTT‑TTTTTCCTGGAATGACCTTCGTCACACCTTTGCATGTACTc                                       >  2:114243/1‑150 (MQ=255)                                                                                                                   gcTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑CGCCTG‑TTTATCCCCACCGAACCCGGGGCCGTTGGTGCTTTTT‑ATGGGCCGTTTTTCGCCCCAGGTA‑TTTACCGCGAAATGACCTTCGCCACACTTTTGCCTGTAACCAtt                                    >  1:1065114/1‑151 (MQ=255)                                                                                                                     tGTTTCTTCCAGTCATCATTATTGGCGGGGTTTTCCCCC‑GGCAAT‑TTTACCACTACAGAGGCC‑GGGACGGGGTGGGTTTTT‑CATCGCTGTTTTTTGCCCCCGTTATTTTACCG‑GAAAATACCTTCGCCACACTTTTGCaatttctcctttg                                  >  1:179968/1‑139 (MQ=255)                                                                                                                     tGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCC‑CGCCCG‑TTTTTCCCCCCCCAAACCGGGGCCGGTGCTGCTTTTT‑ATGCGCTGTTTGTCTCCCCAGTTA‑TTTTCAGGGAAATGAACTTCGCCCCACCTTTGCATGTACTCAtt                                    >  2:978072/1‑149 (MQ=255)                                                                                                                      gTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTg                                   <  1:1152518/149‑1 (MQ=255)                                                                                                                            tcCAGTCATCATTATTGGC‑GGGCTTTGCTCC‑CGCCCG‑TTTTCCCCCACCCAAACCGGGGGCGGTGGTGCTTTTT‑TTTCGCTGTTTTTTGCCCCAGTTA‑TTTTTCGTGAAAAGACCTTCCCCCCCCTTTTGCATTTACTCATTTGCgcgggg                           >  2:1220685/1‑149 (MQ=255)                                                                                                                               aGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGaaa                        <  2:1263439/151‑1 (MQ=255)                                                                                                                                  catcatTATTGGC‑GGGTTTCGCTCC‑CGCCTG‑TTTACCCCCACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTTGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTTGCATGTACTCATTGGCGCGGCGaaaa                       >  1:3100/1‑149 (MQ=255)                                                                                                                                        tattGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTc                <  2:905965/150‑1 (MQ=255)                                                                                                                                         attGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAg              <  1:201671/151‑1 (MQ=255)                                                                                                                                          ttGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAg              >  2:207003/1‑150 (MQ=255)                                                                                                                                            ggC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTgg           <  2:568930/151‑1 (MQ=255)                                                                                                                                                      cGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTc  >  2:207182/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                           GGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCC‑GGCCTG‑TTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCT‑ATGCGCTGTTTGTCGCCACAGTTA‑TTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTC  >  NC_000913/3745346‑3745644

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