Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	38.6%	96.2 / 25.8	58	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (4/31);  new base T (22/0);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGA  >  NC_000913/1413586‑1413882                                                                                                                                                      |                                                                                                                                                       aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtat                                                                                                                                                       <  2:880724/151‑1 (MQ=255)     gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAt                                                                                                                                                   <  1:216709/151‑1 (MQ=255)       tttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCGCCTTCCTTTTCAATAGTGGCGGTAATTTTATACATg                                                                                                                                                  <  1:916461/150‑1 (MQ=255)         tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑aa                                                                                                                                               <  2:1086351/150‑1 (MQ=255)           cTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAGTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATG‑Aaca                                                                                                                                             >  1:226835/1‑150 (MQ=255)            tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtattatata                                                                                                                                                 >  2:1325325/1‑140 (MQ=255)              tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCt                                                                                                                                          <  2:800705/150‑1 (MQ=255)                tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATT‑AAAACcccc                                                                                                                                        >  1:612457/1‑150 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCa                                                                                                                                       <  1:444093/150‑1 (MQ=255)                      ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCatta                                                                                                                                    <  1:1022763/148‑1 (MQ=255)                      ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCatta                                                                                                                                    <  2:929957/148‑1 (MQ=255)                       tGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAAATTTTTAAAAT‑AAAACCCCCcttatt                                                                                                                                  >  1:518061/1‑149 (MQ=255)                        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCattatta                                                                                                                                 <  1:168811/149‑1 (MQ=255)                                cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATT‑AACACccccccttttttttt                                                                                                                             >  1:1086635/1‑134 (MQ=255)                                 tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTACATG‑AACACCTCCCTTATTT‑TTTTCCGTgg                                                                                                                      >  2:149965/1‑150 (MQ=255)                                  gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGtt                                                                                                                    <  2:1048933/151‑1 (MQ=255)                                  gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGt                                                                                                                     <  1:1051328/150‑1 (MQ=255)                                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTa                                                                                                              <  1:1326450/150‑1 (MQ=255)                                                tgcgagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTGTTTATC‑TACAACTCCAATTTTA‑TTTTCACAGGTTTGCTTAttttct                                                                                                        >  2:530441/4‑145 (MQ=255)                                                tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAAG‑AAAAACCCCCTTATTA‑TTTTCAGTGGGTTGGTTTTTTCATCt                                                                                                      >  1:1272503/1‑151 (MQ=255)                                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATG‑AAAACCCCCCTTATTT‑TTTTCCATGGGTCGTTTTTTTTATCttt                                                                                                    >  1:946136/1‑149 (MQ=255)                                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATT‑AAAAACCCCCTTTTTT‑TTTTCCGTGGGTCGTTTTTTTCTTCTTTct                                                                                                  >  2:397413/1‑150 (MQ=255)                                                       aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                                                                                               <  2:270901/151‑1 (MQ=255)                                                             tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCATTTTATTAACT‑AAAACCTCCACTTTTATTTTTAATTGGTTCGTTT‑TTTCTTATTTTGAGTGgttct                                                                                         >  2:118768/1‑151 (MQ=255)                                                               gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAAAA‑AACAACTCCCTTATTA‑TTTTCACTGGGTTGTTTATTTCCTCTTTCTGGGGCTTCttt                                                                                       >  1:1257934/1‑151 (MQ=255)                                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                     <  2:508742/151‑1 (MQ=255)                                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                     <  1:1232210/151‑1 (MQ=255)                                                                   tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc                                                                                    <  2:1059176/150‑1 (MQ=255)                                                                   tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGAAATTTTTTACATA‑AAAACCCCCATTATTA‑TTTTTAGTGGTTCGTTTTTTTCATCTTTTCAGGGGTTCtttttt                                                                                    >  1:872930/1‑149 (MQ=255)                                                                    cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATATT‑AAAACCTCCCTTTTTT‑TTTTCCGTGGGTCGTTTTTTTCCTCTTTTGAGTGGtttttttt                                                                                     >  1:85457/1‑148 (MQ=255)                                                                      gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑c                                                                                 <  2:1338217/149‑1 (MQ=255)                                                                          tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTcca                                                                            <  2:51756/150‑1 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATT‑AAAACCCCCCTTTTTT‑TTTTCATTGGGTTGTTTTTTTCATCTTTTTAGTGCTTCtttttttctttcccc                                                                           >  2:629074/1‑141 (MQ=255)                                                                            gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTccacca                                                                         <  2:264272/151‑1 (MQ=255)                                                                              tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACATG‑AACACCACCCCTATTA‑TTTTCAATGGTTTGTTTTTTTCATATTTCGAGGGGTTCTTTTTTT‑CTTccacc                                                                          >  1:396379/1‑148 (MQ=255)                                                                                  ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAAcc                                                                    <  2:520004/150‑1 (MQ=255)                                                                                  ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAAG‑AAAACCCCCCTTATTT‑TTTTCTGTGGGTCGTTTTTTTCATCTTTTGAGGGGttttttttttctt                                                                                >  1:1179253/1‑134 (MQ=255)                                                                                       ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAAACCCCCTTTTTTT‑TTTTCCTGGGTTGGTTTTTTCC‑TCTTTTTGGGGCTTTTTTTTTT‑TTTCCCCCATCAACGGttt                                                               >  2:922355/1‑149 (MQ=255)                                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATTAATA‑AACACCCCCATTTTTT‑TTTTCATTGGGTTGGTTTTTTTATCTTTTCAATGGTTCTTTTTTT‑TTTTCCCCCTAAACTGTTCgggg                                                           >  2:958258/1‑150 (MQ=255)                                                                                            aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAAT‑AACAACCCCCTTATTT‑TTTTCCATGGTTTGTTTTTTTCCTTTTTTTAGGGCTTTTTTTTTACCTTCCCCCCAAACCGGTCCgggtt                                                          >  1:219806/1‑149 (MQ=255)                                                                                               atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAAT‑TAAACCCCCCTTATTT‑TTTTCATGGGGTCGTTTTTTTCATCTTTTTTGTGTTTTTTTTTTT‑CTTTCACCATAAACCGTTCGGGAtttg                                                       >  1:1172787/1‑148 (MQ=255)                                                                                               atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTTATATA‑GACAGCTCCCTTATTT‑TTTTTACAGGTTTGTTTTTTTTCTTTTTCTAGTGGCTCTTTTTTT‑TCTCCCCCATATACCCTTCGGGAtt                                                         >  1:236696/1‑148 (MQ=255)                                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGc                                                      <  1:626543/147‑1 (MQ=255)                                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAGTTTTTTCAAC‑AAAACCCCCATTCTTT‑TTTTCAATGGGTCTTTTTTTTTATCTTTCGCGCGGTTTTTTTTTT‑TTTTTACCCACACCGGTTCGGGGtttgt                                                      >  2:1253250/1‑144 (MQ=255)                                                                                                    gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATG‑AAAACCCCCCTTTTTT‑TTTTCCGTGGGTTGTTTTTTTCCTCTTTTTCGTGGTTTTTTTTTT‑CTTTCCCCCTAACCGGTTCGGGGTTTGCATcc                                                  >  1:21133/1‑150 (MQ=255)                                                                                                    gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATcc                                                  <  2:390857/150‑1 (MQ=255)                                                                                                      ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGtt                                               <  1:267038/151‑1 (MQ=255)                                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTaca                                            <  2:612719/151‑1 (MQ=255)                                                                                                           tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACAAT‑AAAAACCCCCTTATTA‑TTTTCAGTGGGTTGTTTTTTTCATCTTTTTGGTGGTTTTTTTTTA‑CTTCCACCATAAACCGTTCGGGGTTCGCCTCCCTTacac                                           >  1:376765/1‑150 (MQ=255)                                                                                                              tagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                        <  1:62304/150‑1 (MQ=255)                                                                                                                ggAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCt                                     <  1:545939/151‑1 (MQ=255)                                                                                                                        aCCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAAG‑AAAAACCCCCCTATTA‑TTTTCCGTGGGTGGTTTTTTTTCTCTTTTTAGGGGTTCTTTTTTT‑CTTTCCCCCTAAACCGTTCGGGGTTTTCCTCCCTTTTACTTCCCTc                                    >  1:634819/1‑144 (MQ=255)                                                                                                                                   cAATAGTGGCGGTAATTTTATACATG‑AACTCCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCTCCATAACCTGTTCGTGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                    <  2:1316791/149‑1 (MQ=255)                                                                                                                                     aTAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATACCCTGAAATGGCt                <  2:24690/151‑1 (MQ=255)                                                                                                                                        gTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc             <  1:787607/151‑1 (MQ=255)                                                                                                                                          ggCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCa            <  1:1240977/150‑1 (MQ=255)                                                                                                                                               tAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt       >  1:564920/1‑150 (MQ=255)                                                                                                                                                aaTTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc     <  1:888513/151‑1 (MQ=255)                                                                                                                                                  ttttATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg   >  2:424687/1‑151 (MQ=255)                                                                                                                                                    ttATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGa  <  1:898699/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATG‑AACACCTCCATTATTA‑TTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCA‑CTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGA  >  NC_000913/1413586‑1413882

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