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

ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413586‑1413885                                                                                                                                                      |                                                                                                                                                          aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTta                                                                                                                                                           <  2:1083892/150‑1 (MQ=255) aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTta                                                                                                                                                           <  1:738979/150‑1 (MQ=255)       tttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATg                                                                                                                                                     <  2:540298/150‑1 (MQ=255)        ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                                                                                                    <  2:986432/150‑1 (MQ=255)         tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTATGaa                                                                                                                                                   >  1:409516/1‑150 (MQ=255)             ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTAATGaaaaacc                                                                                                                                              >  1:625606/1‑146 (MQ=255)             ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATaaa                                                                                                                                                   >  2:1084306/1‑146 (MQ=255)                tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCTcc                                                                                                                                            >  2:512248/1‑150 (MQ=255)                tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAAcccc                                                                                                                                             >  2:888986/1‑147 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTTAAACCTcc                                                                                                                                            >  1:1180615/1‑149 (MQ=255)                   tctctGGAAACGCCTGCTTCTTTTTTTCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCATTTATTAAGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACTCCTCCa                                                                                                                                           <  1:203470/148‑1 (MQ=255)                      ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCatta                                                                                                                                        <  1:947131/148‑1 (MQ=255)                              gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAAAAACTCCCTTAtttttttcctt                                                                                                                             >  2:582197/1‑146 (MQ=255)                                 tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATATTAAAACCCCCCTTAtttttt                                                                                                                                  >  1:246494/1‑143 (MQ=255)                                  gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGt                                                                                                                          <  1:377640/150‑1 (MQ=255)                                        tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAAAACCTCCCTTATTATTTTCAGTGgtttgttt                                                                                                                    >  2:1184183/1‑150 (MQ=255)                                        tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGttt                                                                                                                    <  2:321692/150‑1 (MQ=255)                                            tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTc                                                                                                                <  1:147872/150‑1 (MQ=255)                                             aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTcc                                                                                                               <  2:601236/150‑1 (MQ=255)                                                  agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                                         <  1:785400/151‑1 (MQ=255)                                                   gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                                                                                                        <  2:491859/151‑1 (MQ=255)                                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCtt                                                                                                          <  1:547612/148‑1 (MQ=255)                                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                                                                                                       <  2:927912/151‑1 (MQ=255)                                                       aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCTATTTTTTAAAATAACACCACCCTTTTTATTTTCAGCGGGTGGtttttttcctttttttcg                                                                                                     >  1:307335/1‑138 (MQ=255)                                                              cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCCCCCTTTTTTTTTTTACTGGTTTGTTTTTTTCCTCTTTTTAGTGcttctt                                                                                             >  1:331481/1‑151 (MQ=255)                                                               gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACACCTCCCTTTTTTTTTTCTGTGGGTTGTTTTTTTCCTTTTTTTTGTGCtttttt                                                                                            >  1:316261/1‑151 (MQ=255)                                                                cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt                                                                                           <  1:94947/151‑1 (MQ=255)                                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAAGAAAACCTCCCTTATTTTTTTTCGTGGTTCGTTTTTTTCCTTTTTTCAGGGCTTCttttt                                                                                          >  2:795122/1‑151 (MQ=255)                                                                  ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTTAAATGAAAACCCCCCTTTTTATTTTTTTTGGTTCGTTTTTTTCTTCTTTCTTGTGCTTCtttttt                                                                                         >  1:128934/1‑150 (MQ=255)                                                                        tttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCCCCCTTTTTATTTTCCGTGGTTTGTTTTTTTCCTCTTTTTAGGGCttttttttttttt                                                                                     >  2:668193/1‑144 (MQ=255)                                                                         ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGATTTTTTACATGAAAACCGCCGGTATTATTTTTCGTGGGTCGTTTATTCTATTTTTACAGTTCtttttttt                                                                                          >  2:892317/1‑143 (MQ=255)                                                                         ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACAACCCCCTTTTTTTTTTCAGGGGTTTGTTTTTTTCCTCTTTCTGGTGGttttttttttctt                                                                                     >  2:1089481/1‑143 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACCCCCCTTTTTTTTTTCCCTGGGTCGTTTTTTTTCTCTTTTTTGTTCTTCTTTTTTAC‑tttccc                                                                                >  1:1138437/1‑146 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTccac                                                                                <  1:1022239/150‑1 (MQ=255)                                                                            gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGGAAACCTCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCCTCTTTTGAGTGGCTTTTTTTTTC‑TTccccc                                                                               >  1:401228/1‑150 (MQ=255)                                                                               aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑aa                                                                          <  1:1035652/151‑1 (MQ=255)                                                                                gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAACACCTCCCTTTTTTTTTTTAGTGGTTCGTTTTTTTTCTCTTTTTTGTGCTTTTTTTTTTT‑TTTCACCAT‑aa                                                                          >  1:679670/1‑150 (MQ=255)                                                                                     gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑Gtt                                                                   <  2:150605/151‑1 (MQ=255)                                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATCAAAAACCCCATTATTTTTTTTAATGGGTTGTTTTTTTCATCTTTTTTGTGCTTTTTTTTTTT‑TTTCACCATAAACCG‑GTTCggg                                                               >  2:742204/1‑151 (MQ=255)                                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGc                                                         <  1:755691/147‑1 (MQ=255)                                                                                                    gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAAAAAACACCCCCCTCATTTTTTTCTATGGTTCGTTTTTTTCCTCTTTTTTGGGCTTCTTTTTTTC‑TTCCCCCCT‑AAAAGCCTTCGGGA‑TCGCATcc                                                     >  2:1165940/1‑150 (MQ=255)                                                                                                       cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACATTAACACCCCCCCTTTTTTTTTCCGTGGGTCGTTTTTTTCCTCTTTTCAGGGCTTTTTTTTTTC‑TTTCACCCT‑AACCC‑GGTCGGGGTTCTCATCCCtt                                                  >  2:279483/1‑150 (MQ=255)                                                                                                              tagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTc                                           <  2:188/150‑1 (MQ=255)                                                                                                              tagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATATTAACACCCCCCTTTTTTTTTTCCCTGGGTTGTTTTTTTCATCTTTTGAGGGGTTTTTTTTTTCTTTCCCCCCA‑AACCG‑GTTTGGG‑TTCGCCTccccttcacct                                            >  1:128606/1‑140 (MQ=255)                                                                                                                 gAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTc                                       <  2:1049391/151‑1 (MQ=255)                                                                                                                            tccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATc                              <  2:153308/149‑1 (MQ=255)                                                                                                                             ccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCTCCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCt                            <  2:975040/150‑1 (MQ=255)                                                                                                                                     aTAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCt                   <  1:559778/151‑1 (MQ=255)                                                                                                                                       aGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCttt                 <  1:940592/151‑1 (MQ=255)                                                                                                                                        gTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑ATCCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                <  1:328744/151‑1 (MQ=255)                                                                                                                                          ggCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAc              <  2:187642/151‑1 (MQ=255)                                                                                                                                            cGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCg            <  1:66184/151‑1 (MQ=255)                                                                                                                                             ggTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt           >  2:504272/1‑151 (MQ=255)                                                                                                                                                aaTTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTa         <  2:82025/150‑1 (MQ=255)                                                                                                                                                  ttttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg      >  2:460106/1‑151 (MQ=255)                                                                                                                                                   tttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg      >  2:572793/1‑150 (MQ=255)                                                                                                                                                     tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   >  1:830264/1‑151 (MQ=255)                                                                                                                                                      ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  2:445283/151‑1 (MQ=255)                                                                                                                                                      ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  >  2:365229/1‑151 (MQ=255)                                                                                                                                                      ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   <  1:605163/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                          ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCG‑GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413586‑1413885

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