Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,200	0	G	A	48.2%	47.2 / 35.1	56	S217F (TCC→TTC)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (23/6);  new base A (0/27);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.56e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTT  >  NC_000913/3556051‑3556341                                                                                                                                                        |                                                                                                                                              ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAggg                                                                                                                                               <  1:877157/150‑1 (MQ=255) ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAggg                                                                                                                                               <  1:253656/150‑1 (MQ=255)       gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACt                                                                                                                                        <  1:677640/151‑1 (MQ=255)         gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTcc                                                                                                                                      <  1:744470/151‑1 (MQ=255)             cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGc                                                                                                                                  <  1:159013/151‑1 (MQ=255)              aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGcc                                                                                                                                 >  2:284919/1‑151 (MQ=255)                                 cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                              >  2:973214/1‑151 (MQ=255)                                 cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                              >  2:937789/1‑151 (MQ=255)                                     tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                          >  1:547348/1‑151 (MQ=255)                                       tccccaccgaaAAAAAGCCGTTGGGTCATATTTTCCCGTCCGACTCCAAGCGAGCCGTTATTACCCGCCCCC‑GGGCCGCGCGGCGGTTTATGAATTTTCG‑TTTCATGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                          <  1:761265/142‑1 (MQ=255)                                        cACCGACAAAAAAAAAGCTTTGGGTAATATTTTCCCTTCCGACTTCAAGCGAGACGTTATTCCCCGCCCCC‑GGGCCGCGCGGCAGTTTATAAATTTTCT‑TTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                        <  2:358814/150‑1 (MQ=255)                                              ccaaaaaaaaCCTTTGGGGGATATTTTCCCTTCCGCCTCAAGGCGAAGCGGTTTTCCCCGGCCCG‑TGGTGCCGGGGCAGGTTA‑GAATTTTTTGGTGCATGCGGAGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                   <  1:684083/145‑1 (MQ=255)                                                 cAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                              >  1:1019160/1‑151 (MQ=255)                                                   aaaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                             >  1:905169/1‑150 (MQ=255)                                                        gCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                        >  2:898057/1‑150 (MQ=255)                                                         cGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                        >  2:99538/1‑149 (MQ=255)                                                          gTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                     >  1:651304/1‑151 (MQ=255)                                                           ttgggTAATATTTTCACTTTCGCCTCCAAGCGAGACGTTATTCCCCGGCCCC‑TGGTCGCGGGGCAGTTTATGAATATTCT‑GTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                      <  2:635694/146‑1 (MQ=255)                                                           ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                    >  2:303602/1‑151 (MQ=255)                                                             cGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                   >  2:897420/1‑150 (MQ=255)                                                              gttGATATTTTCTCTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCC‑GGGTCGCGCGGCAGTTTATGAATATTCT‑GTTCATGCGGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                   <  2:773286/147‑1 (MQ=255)                                                              ggTGATTTTTTCCCTTCCGACTTCAAGCGAGACGGTATTCCCCGGCCCC‑TGGTCGCGCGGCAGTTTATGAATTTTCT‑GTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                     <  2:13122/147‑1 (MQ=255)                                                               gtatttttttCACTTTCGACTTCAGGCAAGCCGGTTTCCCCCGGCCCC‑GGGTCGCGCGGAGGGTTATAAATATTCT‑GTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                  <  2:869424/143‑1 (MQ=255)                                                                   ttttttCTCTTTCGCCTTCAAGCGAGACGTTATTCCCCGGCCCC‑TGGTCGCGCGGCAGTTTATGAATTTTCT‑GTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                               <  2:888057/146‑1 (MQ=255)                                                                   ttttttCACTTGCGACTACAAGCGAGACGGTATTCCCCGCCCCC‑GGGCCGCGGGGCGGGTTATAAATATTCT‑TTTCAGGGAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                              <  1:238735/147‑1 (MQ=255)                                                                     ttttCCCTTTCGACTACAAGCGAGACGTTATCCCCCGGCCCC‑GGGCCGCGCGGCGGGTTATAAATATTTT‑TTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                            <  2:603108/149‑1 (MQ=255)                                                                     ttttCACTTTCGACTCCAAGCGAGACGGTATTCCCCGGCCCC‑TGGTCGCGCGGCAGTTTATGAATATTCT‑GTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                            <  2:518827/149‑1 (MQ=255)                                                                      tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                         >  2:512200/1‑151 (MQ=255)                                                                       ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                        >  2:350965/1‑151 (MQ=255)                                                                            tttCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                   >  2:562181/1‑151 (MQ=255)                                                                             ttCGTCTACAAGCGAGTCGTTTTTCCCCGCCCCC‑TGGCCGCGCGGGAGGTTATGAATTTTCT‑TTTAATGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                   <  1:831033/150‑1 (MQ=255)                                                                               cGACTTCAAGCGAGACGTTA‑TCCCCGGCCCCTTGGTCGCGCGGCAGTTTATAATTTTTTT‑TTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                   <  1:797829/148‑1 (MQ=255)                                                                               cGACTCCAAGCGAGACGGTATTCCCCGCCCCC‑GGGCCGCGCGGCAGGTTATGAATTTTCT‑GTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                  <  2:510241/149‑1 (MQ=255)                                                                                     cAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                          >  1:101399/1‑151 (MQ=255)                                                                                      aaGCGAGACGGTATTCCCCGCCCCG‑TGGGCGCGCGGCAGTTTATGAATATTCT‑GTTCATGAAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                         <  1:167796/151‑1 (MQ=255)                                                                                          gagaCGTTATTCCCCGCCCCC‑TGTCCGCGCGGCAGTTTATGAATTTTTT‑TTTCATGAAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                       <  1:445317/149‑1 (MQ=255)                                                                                          gagaCGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                     >  2:464192/1‑151 (MQ=255)                                                                                          gagaCGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                     >  1:653045/1‑151 (MQ=255)                                                                                                 tATTCCCCGCCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                <  1:336425/149‑1 (MQ=255)                                                                                                    tACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                           >  1:538304/1‑151 (MQ=255)                                                                                                     tcccggcacccggcgcgcgcgGGCGGTTATTACTTTTCT‑GTTCAGGAAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                            <  2:72982/135‑1 (MQ=255)                                                                                                     ccccGCCCCC‑TGGTCGCGGGGCATTTTATAAATTTTCT‑TTTCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                          <  1:1024485/150‑1 (MQ=255)                                                                                                      cccGGCCCC‑TGGTCGCGGGGCAGGTTTTGAATTTTCT‑GTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                           <  1:299739/149‑1 (MQ=255)                                                                                                      cccGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                         >  2:1007204/1‑151 (MQ=255)                                                                                                      cccGGCCCC‑TGGGCGCGCGGCAGTTTATGAATATTCT‑GTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                           <  1:978486/149‑1 (MQ=255)                                                                                                           cccc‑GGGTCGCGCGGCAGTTTATGAATATTCT‑GTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCa                                       <  1:388357/148‑1 (MQ=255)                                                                                                              c‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                 >  1:552778/1‑151 (MQ=255)                                                                                                                    cgggcgcaagtttATGAATTTTCT‑TTTCGTGCAGAAAAACCCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCaa                                      <  1:31664/130‑1 (MQ=255)                                                                                                                     gcgcgGCAGTTTATAAATTTCCT‑GTCCATGGGGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                            <  1:62794/150‑1 (MQ=255)                                                                                                                     gcgcgGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                           >  1:130083/1‑151 (MQ=255)                                                                                                                           cAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                     >  1:812972/1‑151 (MQ=255)                                                                                                                               ttATGAATATTCT‑GTTCATGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                  <  2:424323/150‑1 (MQ=255)                                                                                                                                   aaaTATTTT‑GTGCATAAAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                   <  1:973193/144‑1 (MQ=255)                                                                                                                                      tgttCT‑TTCCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg           <  1:169537/148‑1 (MQ=255)                                                                                                                                        ttCT‑GTTCATGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc         <  2:131143/150‑1 (MQ=255)                                                                                                                                              ttCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag      <  2:1015783/148‑1 (MQ=255)                                                                                                                                               tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt  >  1:474360/1‑151 (MQ=255)                                                                                                                                                  ggCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg        <  1:354114/141‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                              CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTT  >  NC_000913/3556051‑3556341

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: ‑