Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,080,856	0	T	G	43.4%	76.7 / 12.1	53	intergenic (+43/‑500)	putP/efeU	proline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base T (9/21);  new base G (23/0);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 5.73e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATG  >  NC_000913/1080709‑1081003                                                                                                                                                      |                                                                                                                                                     tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCCCTATCATTCGGGTCCGCCGTCACGGTTGCAGGAAAGCTAAGTGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGtt                                                                                                                                                     <  1:982381/149‑1 (MQ=255)     gTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGATTTAGCA‑GCAGGGTTggccc                                                                                                                                                >  2:367959/1‑147 (MQ=255)      tttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCA‑GC‑GGGGTGCG‑CTcc                                                                                                                                             >  2:902783/1‑151 (MQ=255)        tGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCGCA‑GC‑GGGTTGGG‑CTcccc                                                                                                                                           >  2:1104824/1‑150 (MQ=255)         gCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTAGCA‑GC‑GGGGTGCG‑Cccccct                                                                                                                                          >  2:353858/1‑149 (MQ=255)               tAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTGTTTGTCCTCTCCA‑GC‑GGGGTGCG‑CCCCCtt                                                                                                                                          >  1:43988/1‑145 (MQ=255)                 aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCA‑GC‑GGGGTGCG‑CTCCCCTAAtgtg                                                                                                                                    >  2:74720/1‑149 (MQ=255)                 aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CCCCCTTAATGTGcc                                                                                                                                  >  2:695147/1‑151 (MQ=255)                 aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGcc                                                                                                                                  <  2:499386/151‑1 (MQ=255)                     gcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTc                                                                                                                                <  1:1327321/149‑1 (MQ=255)                      ccGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTTTTTTGCTTCTGCC‑GC‑GGGGTGGGCCCCCCTTAAT‑TGGCTCGcc                                                                                                                             >  1:989537/1‑151 (MQ=255)                       cGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTTCTTCCGCA‑GC‑GGGGTGGG‑CTCCCTTTAAATGCCCCGccc                                                                                                                            >  1:404643/1‑150 (MQ=255)                        gTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGGTGCG‑CCCCCCTTATTTGCCTcccc                                                                                                                             >  2:978747/1‑149 (MQ=255)                         tCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTTTTTTGCTTTCGCA‑GC‑GGGGTTCG‑CCCCCCTTATGTGcc                                                                                                                                  >  1:682476/1‑143 (MQ=255)                          cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCata                                                                                                                          <  1:33134/150‑1 (MQ=255)                            gcggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATaa                                                                                                                       <  2:779258/151‑1 (MQ=255)                               gcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAAtt                                                                                                                    <  2:475270/151‑1 (MQ=255)                                 gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                                                  <  1:1402504/151‑1 (MQ=255)                                  aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCCGCA‑AC‑GGGGTTCG‑CCCCCTTAATGGGCCTCGCCCAATaaaat                                                                                                                    >  2:517254/1‑146 (MQ=255)                                  aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGGTCTACA‑GA‑GGGGTGGG‑CGCCCCTAAAGTGGCCTGCCCTATaaa                                                                                                                      >  1:1366367/1‑146 (MQ=255)                                   tGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCA‑GC‑GGGGTGCG‑CCCCCTTAATGTGCCTCGCCCTATAAATTTaa                                                                                                                 >  2:12559/1‑150 (MQ=255)                                     cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATgg                                                                                                              <  2:671670/151‑1 (MQ=255)                                     cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATgg                                                                                                              <  1:1383551/151‑1 (MQ=255)                                     cAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTGGTGTCAGCA‑GG‑GGGGAGCG‑CTCCCCTTAAATGTCCCTCCACATAAATTTaa                                                                                                                 >  1:756837/1‑148 (MQ=255)                                                       gCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGCTCAGCA‑GC‑GGGGTTCG‑CCCCCTTTATTTTGCCCGCCCTATTAAATGGATGGGGCAGGGGGCGCGCCggg                                                                                            >  1:850239/1‑151 (MQ=255)                                                          gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCA‑GC‑GGGGTGCG‑CTCCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAggggg                                                                                          >  2:147712/1‑150 (MQ=255)                                                          gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCA‑GC‑GGGGTGCG‑CCCCCTTTATGTGGCTCGCCATATAAATTGAAAGGTGCAGGGGGGGCGCAggggg                                                                                          >  1:762179/1‑150 (MQ=255)                                                          gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                                         <  2:1302452/151‑1 (MQ=255)                                                           aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCA‑GC‑GGGGTGCG‑CTCCCTTTATGTGCCCCCCCCTATAAATTTAATGGGGCAGGGGGCGCGCAGGGggcg                                                                                        >  2:932369/1‑151 (MQ=255)                                                           aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                                         <  1:451461/150‑1 (MQ=255)                                                               gcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCATCA‑GC‑GGGGTGCG‑GCCCCTTAATGTGTCCCGCCATATTAAATGAAATGTGCAGGGAGCGCGCAGCgggggg                                                                                       >  2:502132/1‑148 (MQ=255)                                                                 aCTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGGTTTACCC‑GC‑GGGGTGCG‑CCCCCTTAATTAGGCTCCCCCCATTAATTTAATGGTGGAggg                                                                                                       >  1:1325471/1‑130 (MQ=255)                                                                      cATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                              <  2:287713/150‑1 (MQ=255)                                                                          cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAACA‑GC‑GGGGTGCG‑CTCCCTTAATGTGCCTCGCCCTATAAAATGGATGGTGCAGGGAGCGCGCAGGGGGCGGCCAAAcgccccc                                                                           >  2:533477/1‑149 (MQ=255)                                                                          cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCA‑CG‑GGGGTGGG‑CGCCCCTGATGGGGCTCCCCCTATTAATTGAATGATGCCGGGAGCGCGGAGGGGGCGGCCAACcgcggccg                                                                          >  1:278491/1‑150 (MQ=255)                                                                          cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGAGAGGGCGCAGGGGGCGGCCAATcgccgccg                                                                          <  2:180379/150‑1 (MQ=255)                                                                                cgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CCCCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGGGGCCCATCGCCGCCGccccc                                                                     >  2:579799/1‑149 (MQ=255)                                                                                  ccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCctgct                                                                 <  1:991564/151‑1 (MQ=255)                                                                                   cgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTTAGCA‑GC‑GGGGTGGG‑CTCCCTTTAACTGCCCCGCCCCATAAAATGAAAGGTGGCGGGGGCGGGCAGGGGGCGGCCCAACGCCGCCGccccc                                                                     >  2:604061/1‑146 (MQ=255)                                                                                       aCGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTccc                                                            <  1:365369/151‑1 (MQ=255)                                                                                              cAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGGTGCG‑CTCCCCTTATGTGCCTCGCCCTATAAATTGAATGGTGCAAGGGGGGCGCCGGGGGGGGGCAATCGCCGCCGCCCCCTGCCGGCCCCGcc                                                        >  2:1076577/1‑148 (MQ=255)                                                                                                 gAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCggg                                                  <  2:717629/151‑1 (MQ=255)                                                                                                     gCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCACCC‑GGGTTCCG‑‑TCCCTTAATGTGCCTCGCCCTATAAAATGAATTGTGCAGGGGGAGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCCGCTGTCCCCGCCTTCgggg                                                 >  2:618952/1‑148 (MQ=255)                                                                                                         aGGGACTTAGCCTGCGGCGGTTTTGTTTTGGTTTAGCA‑GC‑GGGGTGCG‑CTCCCCTAATTTGGCTCGCCATATAAATTGAATGGGGCCGGGGGGGCGCCAGGGGGGGCCAAACGCCCCCGCCCCCCGCCGTTCCGGCCCTCgggg                                                 >  1:1221169/1‑144 (MQ=255)                                                                                                            gACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAg                                       <  1:596537/151‑1 (MQ=255)                                                                                                                    cTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGAtttt                                <  2:228118/150‑1 (MQ=255)                                                                                                                           ggTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccac                        >  1:24758/1‑151 (MQ=255)                                                                                                                              tttgtttgGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacca                      >  1:797258/1‑150 (MQ=255)                                                                                                                                   ttgGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAccc                >  2:766162/1‑151 (MQ=255)                                                                                                                                    tgGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCg               >  1:161931/1‑151 (MQ=255)                                                                                                                                             ca‑gc‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGttat      <  2:623410/151‑1 (MQ=255)                                                                                                                                                gc‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGttatta    <  2:198121/151‑1 (MQ=255)                                                                                                                                                 c‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTtattat   <  2:1119810/151‑1 (MQ=255)                                                                                                                                                   gggTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATg  <  1:367946/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCA‑GC‑GGGTTGCG‑CTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACACCCGAGCTGTTATTATG  >  NC_000913/1080709‑1081003

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