Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,112,586	0	C	G	44.1%	76.5 / 24.7	59	A91G (GCG→GGG)	yiiR	DUF805 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (11/22);  new base G (26/0);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 3.54e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAA  >  NC_000913/4112440‑4112731                                                                                                                                                     |                                                                                                                                                        cGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑t                                                                                                                                                   <  2:841624/150‑1 (MQ=255)     atTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTa                                                                                                                                               <  2:431808/150‑1 (MQ=255)     atTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAACA‑GGGG‑G‑G‑tttt                                                                                                                                                >  1:434131/1‑149 (MQ=255)     atTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑TGGC‑G‑G‑ttttt                                                                                                                                               >  1:401738/1‑149 (MQ=255)      ttCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAACT‑TGGG‑G‑G‑TTTTTcc                                                                                                                                             >  2:102101/1‑151 (MQ=255)      ttCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑TGGG‑G‑G‑ttttttc                                                                                                                                             >  1:526548/1‑148 (MQ=255)      ttCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑ttttt                                                                                                                                               >  2:87572/1‑148 (MQ=255)         aGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCggg                                                                                                                                          <  2:411564/151‑1 (MQ=255)              gcggcgTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑t                                                                                                                                    <  1:545088/151‑1 (MQ=255)                      ttGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAACT‑GGGG‑G‑G‑TTTTTCCGGGGGGT‑T‑GGCCATgg                                                                                                                           >  2:515018/1‑151 (MQ=255)                           ttGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTCCCGGGGTC‑T‑GGGAATggggggg                                                                                                                      >  2:733990/1‑150 (MQ=255)                             gtgtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTg                                                                                                                     <  2:597535/150‑1 (MQ=255)                              tgtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTgg                                                                                                                    <  1:123678/150‑1 (MQ=255)                                    tgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGtt                                                                                                              <  1:148933/150‑1 (MQ=255)                                       tctGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAA‑GGGG‑G‑G‑TTTTACCCGGGGTC‑TGGGGAAGGGGGGGGGGGGGGttttt                                                                                                           >  2:181504/1‑150 (MQ=255)                                        ctGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑GAA‑TTTTACGGGGT‑TC‑T‑GGGCAAGGGGGGGGGGGCGGTTTGt                                                                                                          >  2:249203/1‑150 (MQ=255)                                         tGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAACT‑GGGG‑G‑G‑TTTTTCCGGGGGTC‑T‑GGGCATGGGGGGGGGGGGGGTTTGttc                                                                                                        >  2:56506/1‑149 (MQ=255)                                                cAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGtt                                                                                                <  1:651072/151‑1 (MQ=255)                                                cAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAAT‑GGGCGG‑T‑TTTTACCGGG‑GTC‑T‑GGGCATGGGGGGGGGGGGGGTTTTTCCC‑GGCggt                                                                                                >  2:770095/1‑149 (MQ=255)                                                 aGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAC‑GGGG‑G‑G‑TTTTTCCGGGGGTC‑T‑GGCAATGGGGGGGGGGGGGTTTTGTCcc‑gccggt                                                                                                >  2:69649/1‑148 (MQ=255)                                                  gCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGGACC‑TGGG‑G‑A‑TTTTTTCAGGGGTC‑T‑GGCAATGGGGGGGGGGGGGGTTTGTtccca                                                                                                     >  1:702756/1‑143 (MQ=255)                                                         tAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTCCCGGGGTCTT‑GGCATGGGGGGGGGGGGGGTTTTTTCCC‑GCGGTTGTTTTTg                                                                                         >  2:684637/1‑150 (MQ=255)                                                            cTGTTAAGCGCTTGCATGATTGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTCCGGGGGGC‑T‑GGCAATGGGCGGGGGGGGGGTTTGTTCC‑GACGGTGGTTTTgggg                                                                                      >  2:817502/1‑147 (MQ=255)                                                                tAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAACTGGGGG‑G‑T‑TTTTACCGGG‑GGC‑T‑GGGAATGGGCGGGGGGGGGGTTTTTccccgcggtggtttttg                                                                                         >  2:18793/1‑129 (MQ=255)                                                                      cTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTCCGGGGGTC‑T‑GGGAATGGGGGGGGGGGGGTTTTTTTCC‑GACGGTTGTTTTGGTGGTGGTGCtttt                                                                           >  2:81270/1‑148 (MQ=255)                                                                       ttGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAACG‑GGGG‑G‑G‑TTTTACCGGGGGTT‑T‑GGGAA‑GGGGGGGGGGGGGGTTTGTCCCGGACGTTGTTTTTggggtggtg                                                                                 >  1:100000/1‑137 (MQ=255)                                                                         gCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGAt                                                                       <  1:145180/151‑1 (MQ=255)                                                                           aTGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAAT‑GGGG‑G‑G‑TTTTTTCGGGTGGC‑TGGGCATTGGGGGGGGGGGGGTTTT‑TCCC‑CACGCTGGTTTTGGGGGTGGTGCTTATCCATc                                                                      >  2:741825/1‑150 (MQ=255)                                                                             gATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTg                                                                   <  2:777842/151‑1 (MQ=255)                                                                               tCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAAT‑GGGG‑G‑G‑TTTTAACCGGGGTC‑T‑TGGAATTGGGGGGGGGGGGGTTTGTTCC‑GGCGGTGGTTTTGGGGGTGATGCTTTTCGAtttt                                                                    >  1:380194/1‑148 (MQ=255)                                                                                 gcgGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTTCGGGGGGT‑T‑GGGAATGGGGGGGGGGGGGGTTTTTTCC‑CACGTTTTTTCTTGTGGTGGTGCtttttggtttttgtg                                                                >  1:298968/1‑137 (MQ=255)                                                                                 gcgGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTACCGGGGGTC‑T‑GGGAATGGGGGGGGGGGGGGTTTTTCCC‑CACGTTGATTTTGGTGGTGATGGTTATCGATCTTtgtg                                                                >  2:816593/1‑150 (MQ=255)                                                                                     gCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAC‑GGGC‑G‑G‑TTTTTCCGGGGGTT‑TGGGCAAGGGGGGGGGGGGGGTTTTTTCCC‑GCGGTTG‑TTCTTGGGGTGGTGGTTTTCGAATTTTGGGCGttt                                                           >  2:805192/1‑151 (MQ=255)                                                                                         tCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTCCCGGGGGG‑T‑GGGAATGGGGGGGGGGGGGGTTTTTTCC‑GCCGGTTGTTCTTGGGGTGGTGGTTTTTGATCTTTGTGGGGttgt                                                         >  1:543171/1‑149 (MQ=255)                                                                                           cGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGTAAGG‑GGCG‑G‑T‑TTTTACCGGG‑GTG‑T‑GGGAATTGGGGGGGGGGGGGTTTTTTCC‑GGCGGTGGTTTTTGGGGTGGTGGTTTTTGATCTTTGTGGGGTTTTTgg                                                      >  1:779748/1‑150 (MQ=255)                                                                                           cGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑GGGG‑G‑G‑TTTTTTCGGGGGGG‑T‑GGGCATGGGGGGGGGGGGGGTTTTTTCC‑GCCGGTGGTTTTTGGGGTGATGGTTTTTGGTCtttggggggttttt                                                        >  1:484082/1‑136 (MQ=255)                                                                                            ggCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGca                                                    <  2:156319/151‑1 (MQ=255)                                                                                                 aTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑TGGG‑G‑G‑TTTTTCCGGGGGTC‑T‑GGGGATGGGGGGGGGGGGGTTTTTTTCC‑CGCGTTGGTTTTGGGGGTGATGCTTTTTGATCTTTGTGCGTTTTTTGGcacac                                                 >  2:271114/1‑149 (MQ=255)                                                                                                     gCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCg                                           <  2:385983/151‑1 (MQ=255)                                                                                                       aTTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAT‑TGGG‑G‑ATTTTTCCCGGGGGGT‑G‑GGGATTGGGGGGGGGGGGGTTTTTTCCC‑GCGGTTGTTTTTGG‑GGTGGTGCTTTTTGGTATTGGTGGGGTTTTTGGGACACaaaggg                                           >  1:594078/1‑145 (MQ=255)                                                                                                             gatgatTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAaaata                                     <  2:277116/149‑1 (MQ=255)                                                                                                                gatTGTGGCGTGGATGCTGCTGGC‑GGGGAACT‑GGGG‑G‑A‑TTTTACCCGGGGGC‑T‑GGGAATGGGGGGGGGGGCGGTTTTTCCC‑CGCGGTGATTCTGGGGGTGATGCTTTTTGGTTTTTGTGGGGTTTTTGGGACACAAAGGGAaaataaa                                   >  2:587863/1‑148 (MQ=255)                                                                                                                 atTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATg                               <  2:349432/151‑1 (MQ=255)                                                                                                                 atTGTGGCGTGGATGCTGCTGGC‑GGGGGACC‑GGGC‑G‑G‑TTTTTCCTGGGGTC‑T‑GGCGAAGGGGGGGGGGGGGTTTTTTCCC‑GCCGCTGGTTCTGGTGGAGAAGCTTTTTGATATTGGGGCGTTTGTTGGCACACAAGGCGGaaataaa                                   >  2:635212/1‑147 (MQ=255)                                                                                                                     tgGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGGGTC‑T‑GGCAATGGGCGGTGGGGGGTTTTGTCCC‑GACGTTGATTCTGGTGGTGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTa                            >  1:778487/1‑150 (MQ=255)                                                                                                                      ggCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTaaa                          <  1:753587/151‑1 (MQ=255)                                                                                                                      ggCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTaa                           <  2:434183/150‑1 (MQ=255)                                                                                                                       gCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAg                         <  2:606449/151‑1 (MQ=255)                                                                                                                            gATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGGGGGGCGTTTTGTCCC‑GACGTTGATTCTGGGGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACt                    >  2:36883/1‑151 (MQ=255)                                                                                                                                 tgctgGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGGGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAgg                >  2:716468/1‑150 (MQ=255)                                                                                                                                  gctgGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAt              <  1:87556/151‑1 (MQ=255)                                                                                                                                  gctgGC‑GGGTAAAT‑GGGC‑G‑A‑TTTTTCCCGGTGTC‑T‑GGCAATGGGCGGGGGGGGGGTTTGTCCC‑GACGTTGATTCTGGTGGTGATGCTTATCGATCTTGGTGCGTTTGTTTGCACACAAGGCGAAAATAAATATGGTAAAGAAACTCAgg                >  1:513678/1‑149 (MQ=255)                                                                                                                                      gC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgtg           >  2:635399/1‑150 (MQ=255)                                                                                                                                          ggTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAg        >  1:827815/1‑150 (MQ=255)                                                                                                                                           gTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGt       >  1:219497/1‑150 (MQ=255)                                                                                                                                             aaCT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGtata    <  2:660612/151‑1 (MQ=255)                                                                                                                                             aaCT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGtat     <  1:522549/150‑1 (MQ=255)                                                                                                                                                  gggC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATaaa  <  2:790598/149‑1 (MQ=255)                                                                                                                                                  gggC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATaaa  <  1:18793/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                        CGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACT‑GGGC‑G‑A‑TTTTACCGGGTGTC‑T‑GGCAATGGGCGGTGGGGCGTTTTGTCCC‑GACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAA  >  NC_000913/4112440‑4112731

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