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

ACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATG  >  NC_000913/1413595‑1413884                                                                                                                                             |                                                                                                                                                        aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAac                                                                                                                                                >  2:71311/1‑150 (MQ=255) aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTaa                                                                                                                                                 >  2:264763/1‑149 (MQ=255)    ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAACAcc                                                                                                                                             >  2:494447/1‑150 (MQ=255)      gtgatgtgtGGAAAAGTATGGTTCTTTTTTACTTGAGAGCATTTTTTCGCATTGTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGaaaaaat                                                                                                                                            >  2:87513/9‑144 (MQ=255)      cTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATTAAAAACTcc                                                                                                                                          >  2:132987/1‑151 (MQ=255)        gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTcc                                                                                                                                          <  2:402723/149‑1 (MQ=255)           ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACcccc                                                                                                                                          >  2:646774/1‑146 (MQ=255)           ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACATTAAAACCCCCatt                                                                                                                                       >  1:893446/1‑149 (MQ=255)             ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCatta                                                                                                                                      <  2:617885/148‑1 (MQ=255)                  aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcc                                                                                                                              <  1:632772/151‑1 (MQ=255)                  aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcc                                                                                                                              <  1:514668/151‑1 (MQ=255)                     gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAAGAACACCTCCATttttttttt                                                                                                                               >  2:707111/1‑146 (MQ=255)                         gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAAAACCTCCCTTTTTTTTTTCCGTg                                                                                                                          >  2:582535/1‑148 (MQ=255)                          cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATTAAAACCTCCATTATTATTTTCCGTgg                                                                                                                         >  1:468905/1‑148 (MQ=255)                            tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCg                                                                                                                     <  1:801254/150‑1 (MQ=255)                                   tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAAGAACACCTCCCTTTTTATTTTCAGTGGGTTGtttttt                                                                                                               >  1:572427/1‑149 (MQ=255)                                    aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTcc                                                                                                             <  2:614162/150‑1 (MQ=255)                                    aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                                                                                                            <  1:169785/151‑1 (MQ=255)                                           agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                                                                                                     <  1:865812/151‑1 (MQ=255)                                           agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                                                                                                     <  1:30223/151‑1 (MQ=255)                                                tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTg                                                                                                 <  2:18023/150‑1 (MQ=255)                                                    tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttc                                                                                             <  2:687890/150‑1 (MQ=255)                                                     cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCTCCCTTTTTTTTTTTCGTGGGTTGTTTTTTTCATCTTTTTTGTGcttctt                                                                                           >  2:611950/1‑151 (MQ=255)                                                     cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTTcttctt                                                                                           <  2:898667/151‑1 (MQ=255)                                                     cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttct                                                                                            <  2:595988/150‑1 (MQ=255)                                                        aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCAATCTTTCGAGTGCTTCttttt                                                                                        <  1:420823/151‑1 (MQ=255)                                                           cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCTCCCTTATTATTTTCCGTGGTTCGTTTTTTTCCTCTTTTTAGTGCTTCttttt                                                                                        >  1:199831/1‑148 (MQ=255)                                                           cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTATAAAAAAAAACCTCCATTTTTATTTCCAGCGGGTTGTTTTTTCCATTTTTTGCGGGGTTAtttttt                                                                                       >  1:573503/1‑148 (MQ=255)                                                               tttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTGGGTGGttttttttttctttc                                                                                 >  1:904978/1‑144 (MQ=255)                                                                      aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAGTAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑a                                                                         <  1:352403/150‑1 (MQ=255)                                                                      aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑aa                                                                        <  2:458415/151‑1 (MQ=255)                                                                       gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AAc                                                                       <  1:430600/151‑1 (MQ=255)                                                                               tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCg                                                                <  1:522004/150‑1 (MQ=255)                                                                                          tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCCCCATTATTATTTTCCGTGGGTCGTTTTTTTCCTTTTTTTTGTGCTTTTTTTTTTCTTTCACC‑AT‑AAACCGGTCGGGGttttcctcc                                                    >  1:477413/1‑144 (MQ=255)                                                                                            tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCCCCCTTTTTATTTTCAGTGGGTTGTTTTTTTCATCTTTCGAGTGCTTTTTTTTTTCTTCCACC‑AT‑AACCGGGTCGGGGTTTGCCTcc                                                    >  1:462345/1‑149 (MQ=255)                                                                                             ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATAAACACCCCCCTTTTTTTTTTCCGTGGGTCGTTTTTTTCATTTTTTTTGGGCTTCTTTTTTTCTTTCCCC‑AT‑AACCCGGTCGGGATTCGCCTCCGtt                                                 >  2:78135/1‑151 (MQ=255)                                                                                              cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCTCCCTTTTTATTTTCCGTGGGTTGTTTATTTCATTTTTTCAGGGCTTTTTTTTTTCTTTCACC‑ATAAACCGGTT‑GGGGTTTGCATCCCtt                                                 >  1:593778/1‑150 (MQ=255)                                                                                              cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTGCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTa                                                <  2:209949/151‑1 (MQ=255)                                                                                                  tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCCCCCCTATTTTTTTCAGGGGGTCGTTTTTTTCATCTTTCTAGTGCTTTTTTTTTTCTTCCCCC‑AT‑AACCCGTTTGGGGTTTTCATCCGtt                                                 >  2:39965/1‑146 (MQ=255)                                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTc                                          <  1:599592/151‑1 (MQ=255)                                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTc                                          <  2:522094/151‑1 (MQ=255)                                                                                                           tGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATTAAAACCTCCCTTATTTTTTTCAGTGGTTCGTTTATTTCCTCTTTCGAGGGCTTCTTTTTTACTTCCACC‑AT‑AAACCGTTCGGGATTCGCCTCCGtt                                                 >  2:661356/1‑137 (MQ=255)                                                                                                            gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAACACCCCCATTATTATTTTCAGGGGGTTGTTTTTTTCATCTTTTGCGTGGTTTTTTTTTACTTTCACCAAA‑AACCGGTTGGGGA‑TCGCATTCCTTacaccttccccgg                                    >  1:529236/1‑140 (MQ=255)                                                                                                             ccACCTTCCTTTTCAATAGTGGCGGTGATTTTTTTAAATAACACCCCCCTTTTTTTTTTTCGGGGGTCGTTTTTTTTTTCTTTTTGGTGGTTCTTTTTTTCTTTCACC‑CT‑AACCAGGTTGGGGGTCGCATCCCTTTTACTTCCCTc                                      >  2:58269/1‑146 (MQ=255)                                                                                                                ccttccttTTCAATAGTGGCGGGAAATTTTTACAATAACACCTCCCTTTTTTTTTTCAGTGGGTTGTTTTTTCCCTCTTTTGGGGGGTTTTTTTTTACTTTCCCC‑AT‑AACCCGGTCGGGGTTCGCATCCGTTACACTTTGCCCCGTaa                                 >  1:757682/1‑148 (MQ=255)                                                                                                                 cttccttTTCAATAGTGGCGGTAATTTTTTTAAAGAAAACCTCCCTTATTTTTTTTAGTGGGTTGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCACC‑CT‑AAACCGGTCGGGGTTTGCATCCGtt                                                 >  2:885205/1‑131 (MQ=255)                                                                                                                  ttccttTTCAATAGTGGCGGTAATTTTTTACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTTCATCTTTCGAGTGCTTTTTTTTTACTTTCACC‑AT‑AAACGGGTCGGGAATCGCATCCGTTACACTTCGCTCGGTAATATcc                            >  1:544830/1‑151 (MQ=255)                                                                                                                    ccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCt                           <  2:843912/150‑1 (MQ=255)                                                                                                                     cttTTCAATAGTGGCGGTAAATTTTTTTATTAACACCTCCCTTTTTATTTTCTGTGGGTCGTTTTTTTCCTCTTTTCAGTGCTTTTTTTTTTCTTTCCCC‑CT‑AAACGGTTCGGGGTTTGCCTCCGTTTCACTTCCCTCCGTAATATTCt                           >  2:126260/1‑149 (MQ=255)                                                                                                                       tttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGa                         <  1:286924/149‑1 (MQ=255)                                                                                                                         tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                      <  2:780721/150‑1 (MQ=255)                                                                                                                              aGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                 <  1:646857/150‑1 (MQ=255)                                                                                                                                tGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc               >  1:608447/1‑150 (MQ=255)                                                                                                                                    ggTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt          >  1:829954/1‑151 (MQ=255)                                                                                                                                     gTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt         <  2:271027/151‑1 (MQ=255)                                                                                                                                     gTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt         <  2:215772/151‑1 (MQ=255)                                                                                                                                            tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg  >  2:592371/1‑151 (MQ=255)                                                                                                                                            tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg  >  2:877581/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                                        ACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACC‑AT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATG  >  NC_000913/1413595‑1413884

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