Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,056,115	0	A	C	42.6%	81.1 / 18.1	55	G75G (GGT→GGG)	glnL	sensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base A (24/7);  new base C (0/23);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 2.69e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGG  >  NC_000913/4055968‑4056262                                                                                                                                                    |                                                                                                                                                     gTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCt                                                                                                                                                   <  2:590668/151‑1 (MQ=255)   gCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTg                                                                                                                                                 <  1:196109/151‑1 (MQ=255)   gCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTg                                                                                                                                                 <  2:771416/151‑1 (MQ=255)      gTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑ccc                                                                                                                                             >  1:174441/1‑151 (MQ=255)        aGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCg                                                                                                                                           <  2:834140/151‑1 (MQ=255)              tCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTcc                                                                                                                                      >  1:419617/1‑150 (MQ=255)               ccggACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGAAAACCCTTG‑CCCCGCCTCCa                                                                                                                                     <  1:287832/147‑1 (MQ=255)                         gcggcgCTGGTTCTCCATCGGACCCACCCCCAGCAGGATCTGCCCGCCCGTCATGCGCGGGGCCGTCACAAAAAGGATATGCGAGCGCCCGCCGATGACCAGCGTCACTTCGTTAGCGAAAAACCCTGG‑CCCCGCCTCCAGACTTTCTTGc                                                                                                                          <  1:526583/151‑1 (MQ=255)                           ggcgCTGGTTATCCACCGGAGCCTTCCCCAGCAGGATCATGCCGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCCCCCGTCGATGACCAGCGTCACTTCTTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTg                                                                                                                           <  2:702198/148‑1 (MQ=255)                                 ggTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCCTCGATGACCAGCGGCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                                  >  1:878587/1‑151 (MQ=255)                                   ttATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCg                                                                                                                 >  1:709687/1‑150 (MQ=255)                                   ttATCCATAGAGCCCACCTCCACCAGGATATGGCCCTCCGAAGTACGCTGGTCCTTCACAGAAAGCATATGCCGGCGCCCCTCGATCCCAAGCTTTCTTTCGTAAGCGTAAAATCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                                  <  1:344535/149‑1 (MQ=255)                                   ttATCCACCGGAGCCACCTCCAGCAGGATCATGCCGTCCGGCATCCGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCATTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCg                                                                                                                 <  1:594792/150‑1 (MQ=255)                                       ccATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAtt                                                                                                            >  2:8345/1‑151 (MQ=255)                                       ccATCGGACCCATCTCCAGCAGGATCATGCCGTCCGGCATCCGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGAAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAtt                                                                                                            <  2:531045/151‑1 (MQ=255)                                         agcGGGGCCATCCCAAGCGGTTTCACGCCGTCCGCCATCCGCTGGGCCGACACAAAAAGGATATGCGAGCGCCCATCGAGGACCAGCTTCACTTCTTTATCGGTAAACCCTTC‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTa                                                                                                          <  2:711255/149‑1 (MQ=255)                                            ggAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCGGGGCCGTCACAAAAAGGATATGCGAGCGCCCTTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTaa                                                                                                         <  2:712142/149‑1 (MQ=255)                                            ggAGCCATCCCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                       <  2:528277/151‑1 (MQ=255)                                            ggAGCCATCCCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCAAAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGAAAACCCTTC‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGata                                                                                                              <  2:88543/144‑1 (MQ=255)                                              aGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                       >  2:314336/1‑149 (MQ=255)                                                ccttctcCAGCAGGACCCTGCCGCCCGGCATCCGCTGGGCCGTCACAGAAGGGATATGCGAGCGCCCGTGGATGACCAGCGTCCCTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGaa                                                                                                   <  2:599942/148‑1 (MQ=255)                                                  aTCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGt                                                                                                 >  2:882286/1‑151 (MQ=255)                                                     cccgaagggaccttGCCTTCCGGCACACGCGGGGCCTTCACGAAAGGGATAGCCGAGCCCCCGCCGACTACCACCGTCTCTGCATTAGCGGAAAACCCTCG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGc                                                                                              <  1:389356/138‑1 (MQ=255)                                                        agcagAATCATGCCGCCGGGCATACGGTGGGCCGTCACAGAAAGGATATGCGACCGCCCGTCGATGACCAGCGTCACTTCGTTATGGTTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                                                                                            <  2:676482/150‑1 (MQ=255)                                                          cagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAAc                                                                                         >  1:187532/1‑151 (MQ=255)                                                            ggATCATGCCGTCCGGCATCCGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAg                                                                                       <  2:418674/151‑1 (MQ=255)                                                            ggATCATGCCGCCCGGCATACGCTGGGCCGTCACAGAAAGGATAGGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                                        <  2:336731/150‑1 (MQ=255)                                                            ggATCATCCCGTCCGGCATCCGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                                        <  2:523703/150‑1 (MQ=255)                                                             gATCTTGCCGCCGGGCTTACGCGGGCCCGTCACAAAAAGGATATGCGAGCCCCCTTCGATGCCCGGCGTCATTTCGTTAGCGGTAAACCCTGG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                                        <  1:403136/149‑1 (MQ=255)                                                                 aTGCCGTCCGGCATCCGCTGGGCCGTCACAGAAAGGTTATGCGGGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCg                                                                                  <  2:60026/151‑1 (MQ=255)                                                                 aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCg                                                                                  >  1:426783/1‑151 (MQ=255)                                                                     cgcccGGCATACGCTGGCCCGTAACAGAAAGGATATGCGAGCGCCCTTCGATGCCCAGCGTCCCTTCGTTATCGAAAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTaa                                                                              <  1:475512/148‑1 (MQ=255)                                                                            cATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtgt                                                                        >  2:418692/1‑150 (MQ=255)                                                                            cATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtgt                                                                        >  1:331272/1‑150 (MQ=255)                                                                                cGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCa                                                                    >  1:311940/1‑150 (MQ=255)                                                                                 gCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCaa                                                                   >  1:221146/1‑150 (MQ=255)                                                                                   tGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAAc                                                                 >  1:766772/1‑150 (MQ=255)                                                                                       ccGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAAttt                                                            >  1:304750/1‑151 (MQ=255)                                                                                             caGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCggg                                                       >  2:157744/1‑150 (MQ=255)                                                                                              aGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGAAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCggg                                                       <  2:536278/149‑1 (MQ=255)                                                                                                     atatGCGGGCGCCCGCCGATGCCCAGCGTCATTTCGTTATCGGAAAACCCTGG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTggg                                              <  1:794587/151‑1 (MQ=255)                                                                                                     atatGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTgg                                               <  2:840649/150‑1 (MQ=255)                                                                                                         gCGAGCGCCCGTCGATGACCAGCTTCATTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTggg                                              <  1:318778/147‑1 (MQ=255)                                                                                                         gCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCg                                            >  1:118979/1‑149 (MQ=255)                                                                                                            aGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGc                                         >  1:145926/1‑149 (MQ=255)                                                                                                               gcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttg                                    >  1:595735/1‑151 (MQ=255)                                                                                                                 ccGTCGATGCCCAGCGTCCCTTCGTTATCGGTAAACCCTGG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                  <  1:815132/151‑1 (MQ=255)                                                                                                                     cGATGACCAGCGTCCCTTCTTTATCGGTAAACCCTTC‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGttgttg                                 <  1:8345/148‑1 (MQ=255)                                                                                                                         gACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGc                           >  2:312536/1‑150 (MQ=255)                                                                                                                           ccAGCGTCACTTCGTTATGGTTAAACCTTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCa                          <  2:131613/149‑1 (MQ=255)                                                                                                                                        gTTCTCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATgg           >  1:366941/1‑151 (MQ=255)                                                                                                                                        gTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATgg           >  1:381563/1‑151 (MQ=255)                                                                                                                                          tATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAt         >  1:461811/1‑151 (MQ=255)                                                                                                                                              ggTAAACCCTTC‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGc      <  2:419726/150‑1 (MQ=255)                                                                                                                                                tAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGcc     >  2:868761/1‑149 (MQ=255)                                                                                                                                                gaaaaCCCTGCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCg       <  1:314244/147‑1 (MQ=255)                                                                                                                                                 aaaaCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAgg  >  2:802101/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGG  >  NC_000913/4055968‑4056262

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