Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,574,051	0	G	T	41.7%	103.0 / 21.4	60	V384V (GTC→GTA)	pqqL	putative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (10/25);  new base T (25/0);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 3.83e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCG  >  NC_000913/1573902‑1574200                                                                                                                                                      |                                                                                                                                                       aGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGa                                                                                                                                                       <  2:567941/151‑1 (MQ=255)     gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGGTTTTAcatc                                                                                                                                                   >  1:597549/1‑151 (MQ=255)           aaaGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGACATCATCg                                                                                                                                               >  2:443988/1‑149 (MQ=255)             aGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGtt                                                                                                                                           <  2:936791/151‑1 (MQ=255)                tGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGttct                                                                                                                                         <  2:590896/150‑1 (MQ=255)                tGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTACATCATCCAGgtctt                                                                                                                                        >  1:292470/1‑151 (MQ=255)                    atGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTcagc                                                                                                                                    <  1:673687/151‑1 (MQ=255)                       tctcTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACCTCCTTGGGGTCTTtagc                                                                                                                                    >  1:505931/1‑148 (MQ=255)                         tctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGGGAATTTACCTCATCGCGTTCTTTAGCAGaa                                                                                                                                >  2:900114/1‑150 (MQ=255)                            tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGaaaa                                                                                                                              <  1:972112/149‑1 (MQ=255)                             ccGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTCCGCAGAAAAAcc                                                                                                                           >  1:596483/1‑151 (MQ=255)                             ccGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCCAGTTTTTTAGCAGaaaaa                                                                                                                             >  2:703388/1‑149 (MQ=255)                              cGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTTTTTCGCAGaaaaa                                                                                                                             >  2:753513/1‑148 (MQ=255)                                  gACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATAGAGTTCTTCAGCAGAAAAAACACGc                                                                                                                       >  1:3385/1‑150 (MQ=255)                                  gACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACAACATCGAGTTTTTTTGCAGAAAAAACATGCt                                                                                                                      >  2:113600/1‑151 (MQ=255)                                    cAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTGCACAATCGCGTTCTTCAGCAGAAAAAACATggt                                                                                                                      >  2:568765/1‑147 (MQ=255)                                        aaCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTTAGCAGAAAAAACAAGCTg                                                                                                                     >  2:767669/1‑146 (MQ=255)                                        aaCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCaa                                                                                                                <  2:55062/151‑1 (MQ=255)                                         aCGGCGTATTATTTAATGAGCTGGATGCCAAGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACACCATCGAGTTTTTTTCCAGAAAAAACAAGCTGATCa                                                                                                                 >  2:855142/1‑149 (MQ=255)                                          cGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTTTTTAGCAGAAAAAACATGCTGTTCaa                                                                                                                >  1:447696/1‑149 (MQ=255)                                            gCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGGTCTTTAGCAGAAAAACCATGCTGAt                                                                                                                   >  2:734717/1‑144 (MQ=255)                                              gTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTCAGCAGAAAAAACATGCTGATTAATGG‑tt                                                                                                          >  1:754104/1‑150 (MQ=255)                                              gTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAAATTTGAAATCATCGAGGTTTTCTGCAGAAAAAACAAGCTGATCaa                                                                                                                >  2:214888/1‑145 (MQ=255)                                                     ttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAAct                                                                                                   <  2:353419/150‑1 (MQ=255)                                                     ttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTCTTTTGCAGAAAAAACATGCTGATCAATGG‑TTGCCAAc                                                                                                    >  2:508866/1‑149 (MQ=255)                                                      tAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTCTTCAGCAGAAAAAACAAGCTGATCAAAGG‑TTGCCaa                                                                                                     >  1:384355/1‑147 (MQ=255)                                                       aaTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGGTTTTACATCCACCAGGTTTTTCGCAGAAAAAACACGCTGATCAATGG‑GTGccca                                                                                                     >  2:971949/1‑144 (MQ=255)                                                          gAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCATCGAGTTTTTTAGCCGAAAAAACATGCTGATCAATGG‑TTGCCAACTCTGcc                                                                                              >  2:481216/1‑150 (MQ=255)                                                               ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTa                                                                                          <  2:18732/149‑1 (MQ=255)                                                                      aGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTTCATCATCGAGGTTTTCAGCAGAAAAAACAAGCTGTTCAATGG‑TTTCCTACCCTTCCCTTAAtttatt                                                                                   >  2:1028267/1‑144 (MQ=255)                                                                         cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCgcag                                                                              <  2:603297/151‑1 (MQ=255)                                                                            cTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAGTTGACAACATCGAGGTCTTCAGCAGAAAAAACATGCTGATCAATTG‑TTGGCCAATCTGCCCTTTAAGCATTCccagc                                                                             >  2:440821/1‑149 (MQ=255)                                                                             tGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGGGATTTTGAATCAATACGGTCTTTAAGAGAAAAAAAAAACTGATCAATGGTTTGCCAA‑TCTCCCCTTAAttccattccaac                                                                             >  2:303009/1‑137 (MQ=255)                                                                                tCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGGAATTTTACATCATCGAGTTTTTTAACAGAAAAACCAAGCTGGTCAAAGG‑TTGCCAAACCTGCCATTAATGCATTCCCAGCAccc                                                                         >  2:606486/1‑149 (MQ=255)                                                                                 cAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTg                                                                       <  2:607279/150‑1 (MQ=255)                                                                                  aGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCa                                                                     <  1:468291/151‑1 (MQ=255)                                                                                          tAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCCAGGTCTTCCGCAGAAAAAACAAGCTGATTAATTG‑TTGCCCACTCTGCCATTAATGCATTCGCCGCATCCTGCAtatt                                                                 >  2:153276/1‑147 (MQ=255)                                                                                                aCGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGt                                                       <  2:83839/151‑1 (MQ=255)                                                                                                 cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTg                                                      <  1:125045/151‑1 (MQ=255)                                                                                                     cAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAtt                                                  <  2:690806/151‑1 (MQ=255)                                                                                                      aGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTTTTTTTCAGAAAAAACATGCTGATCAATTG‑TTGCCAACCCTTCCCTTAATTCATTTGCCGCAACCTGCATATTATCGTCTCCGGCAttt                                                 >  2:1017771/1‑151 (MQ=255)                                                                                                       gCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTa                                                <  1:665322/151‑1 (MQ=255)                                                                                                          tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTaca                                              <  1:581251/150‑1 (MQ=255)                                                                                                          tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTaca                                              <  1:440831/150‑1 (MQ=255)                                                                                                           gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACg                                            <  2:434032/151‑1 (MQ=255)                                                                                                            cTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGCTAGATTTGGCATTATCCGGTTCTTCCGCAGAAAAAACACACTGATCAATCG‑GTGGCAACTCTGCCATTAATACATTTTCAACATCCTCCATAATTTCGtctccttaattt                                                 >  1:341274/1‑137 (MQ=255)                                                                                                                   aCCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACa                                    <  1:5028/151‑1 (MQ=255)                                                                                                                       cATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGa                                <  1:214893/151‑1 (MQ=255)                                                                                                                        aTTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATCATCGAGGTTTTTTGCCGGAAAAACATGCTGGTTAATTG‑TTTGCAACTCTTCCCTTTATGCATTCGCAGCATCCTGCATATTATTTTCTCGGGCCTTTTCACGGgaaaaaa                                    >  2:95554/1‑144 (MQ=255)                                                                                                                          ttttCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCCCGTTTTTTAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCCTTAATGCATTTGCCGCATCCTGCATATTATTCTCTCGTGCATTTACACGGGAAAAAAACGaat                              >  2:834538/1‑148 (MQ=255)                                                                                                                           tttCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAGAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTg                             <  2:111475/150‑1 (MQ=255)                                                                                                                                   aGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCggg                    <  2:940548/151‑1 (MQ=255)                                                                                                                                          gCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtt             >  2:1006604/1‑151 (MQ=255)                                                                                                                                           cGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAttt            >  1:870005/1‑151 (MQ=255)                                                                                                                                             ggTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTa          <  2:188392/151‑1 (MQ=255)                                                                                                                                             ggTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTa          <  1:919224/151‑1 (MQ=255)                                                                                                                                               tAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAAc        >  2:516847/1‑151 (MQ=255)                                                                                                                                                    ttGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTg    <  2:384346/150‑1 (MQ=255)                                                                                                                                                     tGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg  >  1:22786/1‑151 (MQ=255)                                                                                                                                                     tGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg  >  1:789502/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       AGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGG‑TTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCG  >  NC_000913/1573902‑1574200

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