Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,506,038	0	A	C	46.8%	101.4 / 12.9	64	F15V (TTT→GTT)	ydcO	BenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base A (27/6);  new base C (0/29);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 4.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAATGTCAGAGTCAGCGTACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGCTC  >  NC_000913/1505890‑1506187                                                                                                                                                        |                                                                                                                                                      cATAATGTCAGAGTCAGCGTACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAcc                                                                                                                                                     <  2:100528/151‑1 (MQ=255)    aaTGTCAGAGTCAGCGTACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAcccc                                                                                                                                                   <  2:630324/150‑1 (MQ=255)         cAGAGTCAGCGTACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAg                                                                                                                                              >  1:397662/1‑150 (MQ=255)             gTCAGCGTACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCg                                                                                                                                         >  1:629438/1‑151 (MQ=255)                    tACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                                                  >  1:792306/1‑151 (MQ=255)                          cGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAAt                                                                                                                            >  1:996612/1‑151 (MQ=255)                           cgccATTGCCCACCCCCAGCCCCGCCACCCAGCAAGAGT‑TTTGTGCGGTGGTGGCCCGGCCAACATT‑CGCTGCTTCCCAATTTTTTCCCGTTAACCGGGGGTAGCAATTTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTgga                                                                                                                                 <  1:85914/144‑1 (MQ=255)                           cccccttcgccgcccccAGCGCCCCCCCCCGGCCAGATA‑TTTGGGCGGGGGTGGCCCCGGCGACACG‑CGCTGGTTGCCAATTTATGGCCGCTGAGCGGGTGTCGAAAATAAATCCCCCCAAAAACCCCGCCAGTAGCgtgggtg                                                                                                                                   <  1:159717/131‑1 (MQ=255)                             ccATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGAGTGGAGGAAt                                                                                                                            >  2:1130128/1‑148 (MQ=255)                                  g‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                     >  2:271944/1‑150 (MQ=255)                                    ccAGCCCCGGCCCCCTCACCCAG‑CAGAAATTTTGGGCAGTGGGGGC‑CCGGGCACAATCCGCTGTTGGCAAGATTTTTCCCGCTGAAGGGGGGTAGCAAATAAATCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGaaaa                                                                                                                      <  1:852893/148‑1 (MQ=255)                                     cacccccgccgcCGTCATCCAGCCAGAGT‑TTTGGGCTGTGGTGGCTCCGGCGACATT‑CGTTCCTTCCCAAATTATTGCCGCTGAACTGGCGTACCAAATAAAACCCCCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGa                                                                                                                  <  2:480169/142‑1 (MQ=255)                                     cacccccAGCGCCGTCACCCAGCAAGAGT‑TTGGTGCAGTGGGGGCTCCGGCAAAAAT‑CGCTGCTTCCCAGATTATTGCCGCTGAACGGGGGTAGCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACa                                                                                                                    <  2:1055387/146‑1 (MQ=255)                                       cccccccccccttCACCCAGCAAGATT‑TTTGTGCGGTGGTGCCCCCGGCGAAACG‑CGCTGCTTGCCAGTTTATTCCCGCGAAACGGGGGTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                                                                                                 <  1:780318/138‑1 (MQ=255)                                        ccccGGCGCCGTCATCCAGCCAGATT‑TTTGTGAAGTGGGGGCCCCGGCGAAAAT‑CGCTGCTTCCCAGATTATTCCCGCTGAACGGGCGTAGCAATTAAATACCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                                                                                                 <  2:804847/149‑1 (MQ=255)                                        ccccAGCGCCGTCACCCAGCCAGAGT‑TTTGTGAAGTGGGGGCTCGGGCGAAAAC‑CGCTGCTTGCAAGATTATTGCCGTTGACCGGGCGTAGCCAATTAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAc                                                                                                                 <  2:755531/149‑1 (MQ=255)                                           cAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAGCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                                             >  1:205628/1‑150 (MQ=255)                                               gcCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                                         >  2:726374/1‑150 (MQ=255)                                                     aTCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGt                                                                                                  >  2:465660/1‑151 (MQ=255)                                                             agagA‑TTTGTGCGGTGGTGGCCCCGGCGACAAT‑CGCTGCTTGCAAGATTATTGCCGCTGAACGGGGGTAGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                                           <  2:716672/150‑1 (MQ=255)                                                                ga‑TTTGGGCGGTGGTGGCTCGGGCGACAAT‑CGCTGCTTGCCAGATTATGGCCGCTGACCTGGCGTAGCAAATTAATCCCGCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                                             <  2:481817/145‑1 (MQ=255)                                                                   tttGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGt                                                                                      >  2:1121344/1‑150 (MQ=255)                                                                   tttGTGCAGTGGTGGCTCCGGCAAAAAT‑CGCTGCTTCCAAGATTATTGCCGCTGACCTGGCGTAGCAAATAAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                                       <  2:128935/149‑1 (MQ=255)                                                                   tttGGGCGGTGGTGGCCCGGGCGAAAAT‑CGCTGTTTCCAAGTTTTTTGCCGCTGAGCGGGGGTGGCAAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata                                                                                          <  1:1110576/146‑1 (MQ=255)                                                                          aGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                                              >  2:436601/1‑151 (MQ=255)                                                                           ggggtggCCCCGGCAAAAAC‑CGCTGCTGGCAAGATTATGGCCGCTGACCGGGGGAGGCAAATAAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAgg                                                                                 <  1:594647/145‑1 (MQ=255)                                                                             ggtggCCCCGGGAAAAAT‑CGTTGCTTCCAAGATTTTTCCCGCTGATCGGGGGAGGCAAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAcc                                                                             <  2:210818/149‑1 (MQ=255)                                                                                ggCTCCGGCAACAAT‑CGTTGCTTCCCAAATTATCCCCCCGAACCGGGGGTAGCAAATAAACCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAg                                                                         <  2:545215/150‑1 (MQ=255)                                                                                    ccGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAttt                                                                    >  2:647639/1‑151 (MQ=255)                                                                                      ggCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCACGCCAGTAGCGTGGGTGGAGGAATAAAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                                                  >  2:748487/1‑151 (MQ=255)                                                                                      ggCAACAAC‑CGCTGTTTGCCAGATTATTGCCGCGGAACGGGGGAAGCAATTAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAttt                                                                    <  2:136387/149‑1 (MQ=255)                                                                                          caaattcgttttttcccAAATTATCGCCGCTAAAGGGGCGTCACTAATAAACCCCGCGAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCTCCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                                                <  2:860570/135‑1 (MQ=255)                                                                                           cAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                             >  1:824184/1‑151 (MQ=255)                                                                                            aaT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                             >  1:679928/1‑150 (MQ=255)                                                                                                  ctgctTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACg                                                       >  2:1099724/1‑151 (MQ=255)                                                                                                      ttccaaGATTTTTGCCGCTGACTGGGCGTAGCAAATAAACCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAAc                                                        <  1:82440/141‑1 (MQ=255)                                                                                                      ttGCCATTTTATTCCCGCTAAACGGGGGAACCCATTAAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtat                                                           <  2:692115/143‑1 (MQ=255)                                                                                                       tGCCATTTTATCGCCGATGACGGGGCCTAGCCATTTAATACCGCCAAACCCACCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACg                                                       <  1:521338/146‑1 (MQ=255)                                                                                                       tGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                                                   >  2:440878/1‑150 (MQ=255)                                                                                                         ccAGATTATTGCCGCTGACCGGGGGTAGCCAATTAACCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                                                   <  2:437178/148‑1 (MQ=255)                                                                                                          cAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                                               >  2:938439/1‑151 (MQ=255)                                                                                                           aGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                                               >  2:475400/1‑150 (MQ=255)                                                                                                           aGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                                               >  2:507653/1‑150 (MQ=255)                                                                                                              ttattGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATa                                            >  1:1008988/1‑150 (MQ=255)                                                                                                              ttattGCCGCTGAACGGGCGTGGCAAATTAAACCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGa                                              <  2:482488/148‑1 (MQ=255)                                                                                                              ttagtccggCTGACCTGGCGTAGACTATTACTGCCTCAAGAAGCCCCGCCAGTAGCGTGGGTGGAGGAATGGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACaa                                                  <  1:748230/136‑1 (MQ=255)                                                                                                                 ttccccttGAACGGGCGTCCCAAATTAAAACCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                                         <  1:38454/143‑1 (MQ=255)                                                                                                                 tggCCGCTGACCGGGCGTACCAAATTAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                           <  2:708292/146‑1 (MQ=255)                                                                                                                  tGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                       >  1:799176/1‑151 (MQ=255)                                                                                                                    ccGCTGATCTGGCTTACCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTcc                                     <  2:694904/151‑1 (MQ=255)                                                                                                                    ccGCTGAGGTGGCGTAGCCAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                       <  2:1028386/149‑1 (MQ=255)                                                                                                                      gCTGACCGGGCGTAGCAAATTAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGt                                   <  2:996736/151‑1 (MQ=255)                                                                                                                         aaaCGGGGGTAGCAAATAAATGCCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                                         <  2:288660/141‑1 (MQ=255)                                                                                                                          aaCTGGCGTAGCAAATAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGttt                                 <  2:847438/149‑1 (MQ=255)                                                                                                                            cTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGcc                             >  2:340816/1‑151 (MQ=255)                                                                                                                             tGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg                            >  2:166992/1‑151 (MQ=255)                                                                                                                               gCGTAGCCAATTAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga                          <  2:322752/151‑1 (MQ=255)                                                                                                                                cGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgg                           >  2:688194/1‑149 (MQ=255)                                                                                                                                 gTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                        >  2:472158/1‑151 (MQ=255)                                                                                                                                      cAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCg                   >  1:126930/1‑151 (MQ=255)                                                                                                                                           aaTCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCAt                <  2:584288/149‑1 (MQ=255)                                                                                                                                                 gCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGGAAAtc        >  1:623014/1‑151 (MQ=255)                                                                                                                                                  ccAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtc        <  2:979763/150‑1 (MQ=255)                                                                                                                                                    aGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAt         <  2:503562/147‑1 (MQ=255)                                                                                                                                                       aaCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcgctc  >  2:180144/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CATAATGTCAGAGTCAGCGTACTGACGCCCATTG‑CCAGCCCCAGCGCCGTCATCCAGCCAGAGA‑TTTGTGCAGTGGTGGCTCCGGCGACAAT‑CGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGCTC  >  NC_000913/1505890‑1506187

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