Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	44.1%	59.0 / 31.2	59	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (26/7);  new base A (0/26);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.12e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGC  >  NC_000913/5009‑5301                                                                                                                                                       |                                                                                                                                               aaTCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag                                                                                                                                                <  1:859802/151‑1 (MQ=255)    catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGaa                                                                                                                                              <  2:215729/150‑1 (MQ=255)    catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                             >  1:730563/1‑151 (MQ=255)      tcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                            >  1:399945/1‑150 (MQ=255)       caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAATATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                          <  1:989760/151‑1 (MQ=255)         gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaaca                                                                                                                                         >  1:785383/1‑150 (MQ=255)             aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                    <  1:154346/151‑1 (MQ=255)                cTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGt                                                                                                                                 >  1:880151/1‑151 (MQ=255)                  aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtc                                                                                                                               >  2:313253/1‑151 (MQ=255)                         atcaaaaTCGGGCCGGTTTTGTTTTTATGCACCCCGGTTTTTTTATGAAAAATTTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAATTGCGGTAATTAAGAGATGAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCg                                                                                                                         <  2:515403/145‑1 (MQ=255)                              aaTCGGGCGGGTTTTGTTTTTATGCAGCCCGGTTTTTTTATGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTaa                                                                                                                    <  1:37129/150‑1 (MQ=255)                               aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                  >  1:824208/1‑151 (MQ=255)                                           ttGTTTTTATGCAGCCCGGTTTTTTTAAGAAAAAATTATGAAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTACTTTAGAGATTAGGATTGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                       <  2:428824/150‑1 (MQ=255)                                               ttttATGCAGCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                     <  1:496195/148‑1 (MQ=255)                                               ttttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                   >  2:166587/1‑150 (MQ=255)                                                             gtttttttAGAAAAAAATTATGAAAAAAAAGAAAAGGAAAAAAGGAAAAATCCTAAAAAAAGGCGGAAATTTGAAGATAAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                    <  1:240182/149‑1 (MQ=255)                                                              cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                    >  1:594581/1‑150 (MQ=255)                                                               ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                   >  1:722285/1‑150 (MQ=255)                                                               ttttttAAGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                    <  1:507666/149‑1 (MQ=255)                                                                 ttttATGAAAAAATTATGAAAAAAAATGACAGGAAAAAAGGAGAAATTCTCAAAAAATGCGGTATCTTAGAGATTAGGATGGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                 <  2:585572/150‑1 (MQ=255)                                                                     aaaaaaaaaTTAGGAAAAAAAAGGAAAGGAAAAAAGGAGAAATTCTAAAAAAAGGCGGTAATTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                            <  2:307824/145‑1 (MQ=255)                                                                        aagaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                          >  2:370717/1‑150 (MQ=255)                                                                        aagaaATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTTAGAGATTAGGTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                             <  1:697354/147‑1 (MQ=255)                                                                        aagaaATTAGGGAAAAAAATGAGAGGAAAAAAGGAGAAATTCTCAAAAAATGCGGTAACTAAGATTTTAGAATTGCGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                            <  2:828122/148‑1 (MQ=255)                                                                          gaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                       >  1:139543/1‑151 (MQ=255)                                                                          gaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                       >  1:583187/1‑151 (MQ=255)                                                                               taggaaaaaaaaGGACAGGAAAAAAGAAAAAATTCTAAATAAAGGCGGTAATTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                    <  2:484284/142‑1 (MQ=255)                                                                               tATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                   >  2:378210/1‑150 (MQ=255)                                                                                agggAAAAAAAGGACGGGAAAAAAGGAAAAATTCCAAATAAATGCGGTAATTTAGAGATAGGGATGGGGGAAAATAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                  <  2:316849/148‑1 (MQ=255)                                                                                  ggAAAAAAAGGACAGGAAAAAAGGAAAAATTCACAAAAAATGCGTTAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                 <  2:375982/149‑1 (MQ=255)                                                                                     aaaaaaTGACGGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTTAGAGATTAGAATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                            <  1:476359/150‑1 (MQ=255)                                                                                      aaaagaaaCAGGAAAAATGGAAAATTACTCAAAAAACGCGGAAATTAAAAGATTATGATAGGAAAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                           <  2:865888/144‑1 (MQ=255)                                                                                      aaaaaTGAAAGGAAAAAAGGAGAAATTCAAAAAAATGGCGGTAATTGAGAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                           <  2:322066/151‑1 (MQ=255)                                                                                       aaaaGGAGAGGAAAAAGGAAAAAAGTCAAAAAAAAGGGGTTAATTGAGAAATTAGGTGGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                          <  2:386348/151‑1 (MQ=255)                                                                                        aaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                         >  2:772597/1‑151 (MQ=255)                                                                                        aaaTGAAAGGAAAAAAGGAAAATTTCTCAAAAAATGCGTAAACTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                         <  2:722206/151‑1 (MQ=255)                                                                                            gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                     >  1:600451/1‑151 (MQ=255)                                                                                               aGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                  >  1:160305/1‑151 (MQ=255)                                                                                                ggaaaaaaGGAAAAATTCTAAATAAATGCGTTAATTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                  <  2:215074/147‑1 (MQ=255)                                                                                                  gAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                           >  1:367820/1‑123 (MQ=255)                                                                                                  aaaaaaGGAAAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                <  2:197063/149‑1 (MQ=255)                                                                                                   aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                               >  2:964537/1‑150 (MQ=255)                                                                                                      aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                           >  2:673443/1‑151 (MQ=255)                                                                                                      aaGGAAAAATCCTAAAAAAAGGCGTTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                           <  2:416370/151‑1 (MQ=255)                                                                                                            aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                     >  2:751266/1‑151 (MQ=255)                                                                                                             aaTTCCAAATAAATGCGGTAACTAAGAGATCAGTATGGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                    <  2:752328/151‑1 (MQ=255)                                                                                                                  tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                               >  1:451853/1‑151 (MQ=255)                                                                                                                    aaTAAATGGGGAAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                             <  1:756284/151‑1 (MQ=255)                                                                                                                       aaaGGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                           <  1:752336/150‑1 (MQ=255)                                                                                                                                aaCTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                     <  2:590741/147‑1 (MQ=255)                                                                                                                                   ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGAAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                >  2:809492/1‑149 (MQ=255)                                                                                                                                   ttAGAGATAAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                <  1:306566/149‑1 (MQ=255)                                                                                                                                    tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg             >  2:621472/1‑151 (MQ=255)                                                                                                                                     agagATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg            >  1:199028/1‑151 (MQ=255)                                                                                                                                       agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc           >  1:134904/1‑150 (MQ=255)                                                                                                                                       aaaTTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt              <  2:134754/145‑1 (MQ=255)                                                                                                                                             ggATTGGGAAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg    <  1:633091/151‑1 (MQ=255)                                                                                                                                             ggATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac      >  2:908884/1‑149 (MQ=255)                                                                                                                                             ggATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca     <  1:929060/150‑1 (MQ=255)                                                                                                                                                  ggggAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc  <  2:822698/146‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGC  >  NC_000913/5009‑5301

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