Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,816	0	G	A	51.8%	11.0 / 79.9	54	intergenic (+75/+144)	yneM/dgcZ	inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (23/3);  new base A (0/28);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 2.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-05

TATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622669‑1622966                                                                                                                                                          |                                                                                                                                                       tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGaa                                                                                                                                                      <  2:537534/150‑1 (MQ=255)   tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGaaaa                                                                                                                                                    <  1:130255/150‑1 (MQ=255)    ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGaaaaaa                                                                                                                                                  <  1:66302/151‑1 (MQ=255)      ccGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGaaaaaaa                                                                                                                                                 >  2:847764/1‑150 (MQ=255)        gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAAcc                                                                                                                                              >  2:431329/1‑151 (MQ=255)               gAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGt                                                                                                                                        >  1:580362/1‑150 (MQ=255)                     ttttttttttTGTTTTTCGGGCCGTGTATTACAGCCAAAAGGGGGATGAATAATGAAGGGAAATAACCCCTCACCTACCCGCCCCT‑TTG‑TTACAGTTGTGTACAGGGGGCT‑TGT‑TTTTTAG‑GAGGG‑AAAAAAAAAACCGCCAGTAAACcg                                                                                                                                  <  2:353583/150‑1 (MQ=255)                     tttttttttCTGTTTTTCGGGCCGTTTTTTCCACCCAAAAAGGGGAGAACAAATAAAGGGAAAAAACCCCCCCCCAACCCGCCCCT‑TTT‑TTACGTTTGTGTA‑AAGGGGCCTTGT‑TTTTTAG‑GCCGGAAAAAAAAAAACCGCCAGTAAAcc                                                                                                                                   <  2:547933/150‑1 (MQ=255)                             tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGaa                                                                                                                          >  1:766143/1‑150 (MQ=255)                               tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGc                                                                                                                       >  1:701271/1‑151 (MQ=255)                                  ttttctggCCGCTTTTTCCACCCAAAAAGGGGATGATTAATGAACGGAAAAACCCCCTCCCCTACCCGCCCCT‑TTT‑TTACATTTGTGAACAGGGGGCT‑TGT‑TTTTTAG‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAAt                                                                                                                         <  2:375680/146‑1 (MQ=255)                                   ttttggcccGCTTTTTCCACCCAAAAAGGGGATGACTAAGGAAGGGA‑AAAACCCCCCCCCCAACCCGCCCCTTTT‑TTACGGTTGTGAACAGGGGCCT‑TGT‑TTTTTAT‑GAGGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                                    <  1:23413/143‑1 (MQ=255)                                    ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                   >  2:136365/1‑150 (MQ=255)                                      ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                                >  2:589985/1‑151 (MQ=255)                                        ggCCGCTTATTTCACCCACAAAGGGGATGAAAAATGAAGGGAAATAACCCCCCCCCAACCCGCCCCT‑TTT‑TTACATTTGTGAACAGGGGGCT‑TGT‑TTTTTAG‑GCCGG‑CAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                                <  1:844592/149‑1 (MQ=255)                                        ggCCGCTTATTCCACCCAAAAAGGGGATGACTAATGAACGGAAAAAACCCCCCCCCAACCCGCCCCT‑TTT‑TTACAGTTGTGTACAGGGGGCT‑TGT‑TTTTTAT‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                                 <  2:155699/148‑1 (MQ=255)                                            gcttATTCCACCCAAAAAGGGGATGAAAAATGAAGGGAAATACCCCCTCCCCAACCCGCCCCT‑TTT‑TTACAGTTGTGAACAGGGGGCT‑TGT‑TTTTTAT‑GAGGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg                                                                                                           <  1:50298/147‑1 (MQ=255)                                              gTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                        >  2:907363/1‑151 (MQ=255)                                               ttttccAGCAAAAAAGGGGAGGACAAAGGAACGGAAATAACCCCTCCCCAACCCGCCCCT‑TTG‑TTACAGTTGTGTAAAGGGGGCT‑TGT‑TTTTTAT‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGa                                                                                                          <  2:538134/143‑1 (MQ=255)                                                 ttccacccAAAAAGGGGATGACAAATGAGCGGAGATAACCCCCCCCCAACCCGCCCCT‑TTT‑TTACATTTGTGTAAAGGGGGCT‑TGT‑TTTTTAG‑GAGGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtg                                                                                                      <  2:74413/144‑1 (MQ=255)                                                 ttacAGCCACAAAGGGTATCACTAAAGAACGGAAATAATCCCTCTCCAACCCCCCCCC‑TTT‑TTAGAGTTGTATAAAGGGGACT‑TTT‑TTTTTAG‑GGCGA‑AAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtg                                                                                                      <  2:372067/147‑1 (MQ=255)                                                    cacccACAAATGGGATAACAAATGACGGGAGAAAACCCCCCCCCT‑ACCGGCCCCTTTG‑TTACAGTTGTGTACAGGGGGCC‑TGA‑TTTTTAT‑GCCGG‑CAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgtg                                                                                                    <  2:180658/146‑1 (MQ=255)                                                     aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                                 >  1:676391/1‑151 (MQ=255)                                                     aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAAAGCTTGCATGGGTAGATTTGTGctt                                                                                                 >  2:253603/1‑151 (MQ=255)                                                       ccACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                               >  2:187862/1‑151 (MQ=255)                                                         aaaaaTGGGATGACTAATGAGGGGAAATAACCCCTCCCCT‑ACCGGCCCCTTTT‑TTACGGTTGTGAAAAGGGGGCT‑TGT‑TTTTTAT‑GGGGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtt                                                                                                  <  1:904857/144‑1 (MQ=255)                                                           aaaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                           >  1:292059/1‑151 (MQ=255)                                                           aaaTGGAATGACTAATGAACGGAAATAACCCCCCCCCT‑ACCCGCCCC‑TTGTTTACAGTTGTGAAAAGGGGGCC‑TTA‑TTTTTAT‑GCCGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                                 <  2:577370/145‑1 (MQ=255)                                                                            aaggaaaaTAACCCCCCCCCAACCCGCCCCT‑TTT‑TTAAAGTGGTGTAAAGGGGGCC‑TGT‑TTTTTAT‑GACGG‑CAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtt                                                                           <  2:607364/143‑1 (MQ=255)                                                                                   aaaacccccccccTACCCGCCCCT‑TTG‑TTACAGTTGTGTACAGGGGGCT‑TGT‑TTTTTAG‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCa                                                                   <  2:44010/140‑1 (MQ=255)                                                                                   aTAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                                    >  2:300230/1‑150 (MQ=255)                                                                                         cccccccTACCCGCCCCT‑TTT‑TTACTGTGGTGAAAAGGGGGCT‑TTT‑TTTTTAG‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              <  1:960755/145‑1 (MQ=255)                                                                                          cTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                             >  2:337891/1‑150 (MQ=255)                                                                                              ccTACCCGCCCCT‑TTT‑TTACTGTTGGGAACAGGGGGCT‑TTT‑TTTTTAG‑GGGGG‑CAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                                          <  1:376068/149‑1 (MQ=255)                                                                                                 acccGCCCCC‑TTT‑TTACGTTTGTTAAAAGGGGGCT‑TGT‑TTTTTAT‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                                       <  1:898248/147‑1 (MQ=255)                                                                                                 aaCCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                                     >  2:296991/1‑151 (MQ=255)                                                                                                     ggCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  >  2:929908/1‑150 (MQ=255)                                                                                                             tG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                           >  1:24486/1‑150 (MQ=255)                                                                                                                ttACAGTTGTGTACAGGGGGCT‑TGT‑TTTTTAG‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                                         <  1:300705/150‑1 (MQ=255)                                                                                                                   aaGTTGTGTACAAGGGGCT‑TTT‑TTTTTAT‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                     <  1:667078/150‑1 (MQ=255)                                                                                                                             aaaGGGGCC‑TTTTTTTTTAG‑GCCGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                               <  1:55752/147‑1 (MQ=255)                                                                                                                               aGGGGCC‑CTA‑TTTTTTTAGGCGG‑CAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                           <  2:189461/150‑1 (MQ=255)                                                                                                                                   gCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                     >  2:910083/1‑151 (MQ=255)                                                                                                                                    cc‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg                    >  1:810643/1‑151 (MQ=255)                                                                                                                                     cttGA‑TTTTTAG‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                     <  2:940090/148‑1 (MQ=255)                                                                                                                                      cttgttttttAT‑GACGG‑CAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                     <  2:889008/145‑1 (MQ=255)                                                                                                                                       tGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  >  2:401279/1‑151 (MQ=255)                                                                                                                                          ttttttAT‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  <  1:495327/148‑1 (MQ=255)                                                                                                                                          ttttttAG‑GAGGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa                 <  1:952314/149‑1 (MQ=255)                                                                                                                                           tttttAT‑GAGGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt                   <  2:292086/147‑1 (MQ=255)                                                                                                                                           tttttAT‑GACGG‑AAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa                 <  1:596412/149‑1 (MQ=255)                                                                                                                                              ttAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg            >  2:158888/1‑151 (MQ=255)                                                                                                                                               tAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg            >  2:480350/1‑150 (MQ=255)                                                                                                                                                aT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg            >  1:335472/1‑149 (MQ=255)                                                                                                                                                          gAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt  >  2:73953/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       TATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGA‑TTTTTAT‑GACGG‑CGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622669‑1622966

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