Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,112,578	0	T	G	41.4%	100.1 / 10.9	58	G88G (GGT→GGG)	yiiR	DUF805 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (5/29);  new base G (24/0);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATCTACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTAT  >  NC_000913/4112430‑4112728                                                                                                                                                      |                                                                                                                                                            aGAATCTACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTaa                                                                                                                                                           <  1:707654/151‑1 (MQ=255)   aaTCTACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑t                                                                                                                                                        <  1:300884/151‑1 (MQ=255)   aaTCTACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGaa                                                                                                                                                           >  2:246073/1‑149 (MQ=255)         cTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑g                                                                                                                                                 <  2:787123/150‑1 (MQ=255)          tCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGaaattgggggt                                                                                                                                                  >  2:442320/1‑142 (MQ=255)            gATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGTAAC‑g‑gggggttt                                                                                                                                               >  2:242297/1‑146 (MQ=255)              tatTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAC‑T‑GGGC‑GG‑ttttt                                                                                                                                          >  1:129690/1‑150 (MQ=255)                             gTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAC‑T‑TGGG‑GAGTTTTACTggggttt                                                                                                                                 >  1:521845/1‑142 (MQ=255)                               tttGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑TGGGGG‑GA‑TTTTTCCggggtcttgg                                                                                                                              >  1:2603/1‑140 (MQ=255)                                ttGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGGAAT‑G‑GGGG‑GT‑TTTTACCGGG‑GTCTgggcaagg                                                                                                                        >  2:570660/1‑145 (MQ=255)                                         gtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAC‑T‑GGGC‑GG‑TTTTTTCCGGGGGGTGGGAATGggg‑ggggggg                                                                                                              >  2:488438/1‑151 (MQ=255)                                          tgcttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑TCGGGG‑GG‑TTTTACCGGGGGCTTGGCAATGggg‑ggggggg                                                                                                              >  1:27698/1‑151 (MQ=255)                                                ctctGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGG‑GG‑TTTTTTCGGGTGTCTGGCCATGGGC‑GGGGGGGGGgttt                                                                                                        >  1:99499/1‑150 (MQ=255)                                                       cAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGACGATTGTGGCGTGGATGCTGCTGGC‑GGGGTAC‑T‑GGGG‑GG‑TTTTTCCGGGGGGGTGGGCATGGGG‑GGGGGGGGGgttttttcc                                                                                                   >  1:119848/1‑143 (MQ=255)                                                        aaCAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑TGGGGG‑GT‑TTTTACCGGG‑GGCTGGCAATGGGG‑GGGGGGGGGGTTTGTccccg                                                                                                 >  2:66401/1‑147 (MQ=255)                                                        aaCAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGG‑GG‑TTTTTCCGGGGGTCTGGGAATGGGG‑GGGGGGGGGgttttttcc                                                                                                   >  2:782251/1‑142 (MQ=255)                                                          cAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGtt                                                                                             <  2:655456/151‑1 (MQ=255)                                                          cAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAC‑T‑GGGG‑GG‑TTTTACCTGGGGGCTGGGAAATGGG‑GGGGGGGGGTTTTTTccccgcggt                                                                                             >  1:277855/1‑145 (MQ=255)                                                                cAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGGTGGC‑GGGGAAA‑T‑GGGC‑GG‑GTTTTCCGGGGGGGTGGGGATGGGG‑GGGGGGGGGttttttgcccgcgt                                                                                              >  1:182569/1‑134 (MQ=255)                                                                 aGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCt                                                                                       <  2:441918/150‑1 (MQ=255)                                                                  gTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTgg                                                                                     <  2:630032/151‑1 (MQ=255)                                                                         ttAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAC‑A‑GGGG‑GG‑TTTTTCCCGGGGGTTGGGAATGGGG‑GGGGGGGGGGTTTTTGCCCCCGTTGGttt                                                                                        >  1:782214/1‑140 (MQ=255)                                                                          tAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGG‑GG‑TTTTTACGGGGGTCTGGCAATGGGG‑GGGGGGGGGGTTTGTTCCGACggttgtttttgtgat                                                                                 >  1:603671/1‑133 (MQ=255)                                                                            aGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGatgatg                                                                             <  1:557535/149‑1 (MQ=255)                                                                                  tGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCg                                                                      <  1:444129/150‑1 (MQ=255)                                                                                     aTGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGTCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATcg                                                                  <  1:703442/151‑2 (MQ=255)                                                                                      tGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑G‑TGGG‑GG‑TTTTTTCGGGTGGTTGGGCAGGGGCGGGGGGGGGGGTTTTTCCCG‑CGTTTGTTCTGGGGGTGATGGTTGTCgggtt                                                                  >  1:405762/1‑146 (MQ=255)                                                                                       gATCGCGGGCGTTCCTGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCtt                                                                 <  1:276690/150‑1 (MQ=255)                                                                                           gcgGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTg                                                             <  1:638303/150‑1 (MQ=255)                                                                                           gcgGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGc                                                            <  1:182533/151‑1 (MQ=255)                                                                                               gCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGttt                                                        <  2:565896/151‑1 (MQ=255)                                                                                                cGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTttg                                                       <  1:352944/151‑1 (MQ=255)                                                                                                    ccGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAACAG‑GGGG‑GT‑TTTTACCGGG‑GGCTGGGCATGGGG‑GGGGGGGGGGTTTTTTCCGGCGTTTGTTTTGGGGGTGGTGCTTATCGGTCTTTGTGGGtttttttg                                                   >  1:50019/1‑149 (MQ=255)                                                                                                         gcATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGcacac                                              <  1:149725/151‑1 (MQ=255)                                                                                                         gcATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGCGGT‑TTTTACCGGGGGT‑TGGGCATGGGG‑GGGGGGGGGGTTTGTTCCGGCGGTGGTTTTGGGGGTGATGGTTTTTGGTCTTGGTggggttttttggaca                                               >  1:604966/1‑136 (MQ=255)                                                                                                           aTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACaa                                            <  1:246954/151‑1 (MQ=255)                                                                                                              ggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGc                                         <  2:319018/151‑1 (MQ=255)                                                                                                              ggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGC‑GG‑TTTTTCCGGGGGTCTGGGAATGGGG‑GGGGGGGGGTTTTTTTCCCACGTTGGTTTTTGGGGTGATGCTTATCGATCTTTGGGGGTTTGTTTGCACCCaaag                                          >  2:97798/1‑148 (MQ=255)                                                                                                               gCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGTAAT‑G‑GGGG‑GT‑TTTTACCGGGGGTTGGGCATGGGGG‑GGGGGGGGTTTTTTTCCCGGGTTTGTTTTTGGGGTTGGT‑CTTTTTGGTCTTGGTGGGGTTGTTTGcaacaaaag                                          >  1:427070/1‑142 (MQ=255)                                                                                                                cATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGa                                       <  1:347640/151‑1 (MQ=255)                                                                                                                cATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑C‑GGGC‑GG‑TTTTTCCCGGGGGCTGGCCATGGGG‑GGGGGGGCGGTTTTTTCCGACGTTGGTTTTGGGGATGATGCTTTTCGATCTTTGGGCGGTTTTTTGGACACAAGGCg                                        >  2:176095/1‑150 (MQ=255)                                                                                                                 aTTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGG‑GG‑TTTTATCGGGGGTCTTGCAATGGGG‑GGGGGGGGGGTTTTTTCCGGCGTTGGTTTTGGGGGTGATGGTTTTTGATCTTTGTGCGGTTTTTTTGACACaaagggg                                       >  1:420829/1‑145 (MQ=255)                                                                                                                   ttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGaaa                                     <  2:287671/150‑1 (MQ=255)                                                                                                                           atTGTGGCGTGGATGCTGCTGGC‑GGGGAAA‑T‑GGGCGGA‑TTTTCCCGGG‑GTCTGGGTATGGGG‑GGTGGGGGGTTTTTTTCCCACGGTGGTTTTGGGGGTGGTGGTTTTTGGTCTTGGTGGGGTTGTTTGCACACAAGGCGGAAATAAaaat                             >  1:675925/1‑150 (MQ=255)                                                                                                                            ttGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATg                            <  2:405793/150‑1 (MQ=255)                                                                                                                             tgtgGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATgg                           <  2:745732/150‑1 (MQ=255)                                                                                                                              gtgGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTa                         <  2:318033/151‑1 (MQ=255)                                                                                                                               tgGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTTCCGGGGGTTTGGCAATGGGG‑GGGGGGGCGTTTTGTCCCGACGTTGATTCTTGGGATGATGCTTATCGATCTTGGTGCGGTTGTTGGCACACAAGGCGAAAATAAATATGGTa                         >  2:89769/1‑150 (MQ=255)                                                                                                                               tgGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGGGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTTTTGGCACACAAGGCGAAaataaata                              >  1:611039/1‑145 (MQ=255)                                                                                                                                ggCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTaa                        <  2:368130/150‑1 (MQ=255)                                                                                                                                        tgctgcTGGC‑GGGGAAA‑T‑GGGC‑GG‑TTTTACCGGGGGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTTGTGGTGATGCTTATCGATCTTGGTGCGTTTTTTGGAAAAAAAGGCGAAAATAAATATGGTAAAGATACTc                >  1:716305/1‑150 (MQ=255)                                                                                                                                            gctgGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGa            <  1:451560/150‑1 (MQ=255)                                                                                                                                            gctgGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAt           <  1:176074/151‑1 (MQ=255)                                                                                                                                             ctgGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtg          >  2:297650/1‑151 (MQ=255)                                                                                                                                               ggC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgtg        <  1:616106/151‑1 (MQ=255)                                                                                                                                               ggC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgt         <  1:427465/150‑1 (MQ=255)                                                                                                                                                 c‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGaa      <  2:217180/151‑1 (MQ=255)                                                                                                                                                     gTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGta   <  2:217777/151‑1 (MQ=255)                                                                                                                                                      tAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGtat  >  2:381558/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                            AGAATCTACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAC‑T‑GGGC‑GA‑TTTTACCGGGTGTCTGGCAATGGGC‑GGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTAT  >  NC_000913/4112430‑4112728

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