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

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGC  >  NC_000913/5012‑5304                                                                                                                                                       |                                                                                                                                                  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGaa                                                                                                                                                 <  1:1036266/150‑1 (MQ=255) catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGaa                                                                                                                                                 <  1:770698/150‑1 (MQ=255)    caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaac                                                                                                                                             <  1:955847/151‑1 (MQ=255)           aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                      >  1:198954/1‑151 (MQ=255)               aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTtc                                                                                                                                  >  2:655664/1‑151 (MQ=255)                ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGA‑TAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                >  1:237094/1‑151 (MQ=255)                   aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCAt                                                                                                                              >  2:496756/1‑151 (MQ=255)                           aaTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                      >  2:1012508/1‑151 (MQ=255)                               aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                  >  1:306638/1‑151 (MQ=255)                                ggCGGGGTTTGTTTTTATGCAGCCCGGTTTTTTTATAAAAAATTAAGGGAAAAAAATGACAGGAAAAAA‑GG‑AAAAATTCAAAAAAAATGCGGTAACTT‑AGAAATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAAtct                                                                                                                    <  1:571299/148‑1 (MQ=255)                                               tATGCAGCCCGTTTTTTTTATGAAGAAATTATGGAGAAAAATAACAGGGAAAAA‑GG‑AGAAATTCTAAATAAATGCGTTAACTT‑AGAGATTAGGATTGGGGAGAATAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                     <  1:877086/148‑1 (MQ=255)                                                   cgccccGGTTTTTTTAAAAAAAAATTATGGAAAAAAATGACAGGAAAAAA‑GG‑AAAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                 <  2:989470/145‑1 (MQ=255)                                                      ccgggtttttttAAAAAAAAATTAGGAAAAAAAAGGAAAGGAAAAAA‑GA‑AAAAATTCTAAAAAAAGGCGGTAACTA‑AAAAATTAGGTTGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                              <  1:284301/142‑1 (MQ=255)                                                            ttttttATGAAGAAATTATGAAAAAAAAAGACAGGGAAAAA‑GG‑AAAAATTCTAAATAAATGCGTTAATAT‑AGAAATGAGGATGGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                  <  1:883142/138‑1 (MQ=255)                                                            ttttttAGAAAAAACTTATGAAAAAAAATGA‑GGGGAAAAA‑GGAAGAAATTCTAAATAAAGGGGATATTTA‑AGAAATTAGGCTGGCGCAAAACAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                         <  2:487816/147‑1 (MQ=255)                                                             tttttATGAAAAAATTATGGAAAAAAATGACGGGAAAAAA‑GG‑AAAAATTCTCAAAAAATGCGTTAACTT‑AGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                     <  2:552879/150‑1 (MQ=255)                                                             tttttATAAAAAAATTATGAAAAAAAATGACGGGAAAAAA‑GG‑AAAATTTCAAAAAAAAGGCGGTAACTT‑AGAGATTAGGATGGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                     <  2:1004958/150‑1 (MQ=255)                                                              ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                    >  2:1008096/1‑150 (MQ=255)                                                                   tGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                               >  1:788706/1‑150 (MQ=255)                                                                     aaaaaaTAATGGAAAAAAATGACAGGAAAAAA‑GG‑AAAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                 <  1:655651/143‑1 (MQ=255)                                                                        aatttATGGAAAAAAATGA‑AGGGAAAAAGGG‑AAAATTTCTAAAAAAAGGGGGTACTTA‑AGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                         <  2:1014179/148‑1 (MQ=255)                                                                          aTTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                       >  2:494994/1‑151 (MQ=255)                                                                             aTGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                     >  1:596321/1‑150 (MQ=255)                                                                                 agaAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                  >  1:51049/1‑149 (MQ=255)                                                                                 agaAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                  >  2:1041579/1‑149 (MQ=255)                                                                                 agaAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                  >  2:367815/1‑149 (MQ=255)                                                                                 aaaaaaaTGACGGGGAAAAG‑GG‑AAAAATTCTCAAAAAATGCGTTAGCAG‑AAAGATTAGGTTCGGGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                     <  1:972004/144‑1 (MQ=255)                                                                                     aaaTGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                            >  1:568962/1‑151 (MQ=255)                                                                                     aaaTGACAGGGAAAAA‑GG‑AAAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                            <  1:120002/151‑1 (MQ=255)                                                                                      aaTGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                            >  1:203775/1‑150 (MQ=255)                                                                                       aggACAGGAAAAAA‑GA‑AAAAATTCTAAAAAAAGGCGGTAATTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                          <  2:992946/149‑1 (MQ=255)                                                                                          aCAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                       >  1:157969/1‑151 (MQ=255)                                                                                             ggaaaaaa‑GG‑AAAAATTCTCAAAAAAGGCGGTAACTT‑AGAAATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                      <  1:819837/146‑1 (MQ=255)                                                                                              ggAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                   >  2:952276/1‑151 (MQ=255)                                                                                               aaaaaa‑GG‑AAAAATTCACAAAAAATGCGGTAATTT‑AGAGATTGGGTTTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                  <  1:496747/150‑1 (MQ=255)                                                                                                aaaaa‑GG‑AAAAATTCTAAAAAAATGCGGTAATTT‑AGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                  <  1:148324/150‑1 (MQ=255)                                                                                                  aaa‑GG‑AAAAATTCTCAAAAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                               <  2:1053128/151‑1 (MQ=255)                                                                                                         agaAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                          >  1:618819/1‑151 (MQ=255)                                                                                                         aaaaaTTCTAAAAAAAGGCGGTAACTT‑AGAAATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                            <  1:177738/147‑1 (MQ=255)                                                                                                           aactccaaaaaaaaTGCGTTAAGTG‑AGATATGAGGAGTGGGAAAAAAAACTCCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                          <  2:874759/138‑1 (MQ=255)                                                                                                            aaTTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                       >  2:758132/1‑151 (MQ=255)                                                                                                            aaTTCTCAAAAAGTGCGTTAACTA‑GGAGTTGGGGGTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                       <  1:879596/151‑1 (MQ=255)                                                                                                                ctcAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                   >  1:377512/1‑151 (MQ=255)                                                                                                                caaaaaaaaTCCGGTAACTA‑AGGATTTGGGAGTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                   <  2:767147/145‑1 (MQ=255)                                                                                                                caaaaTAAATGCGTTAACTT‑AGAGATGAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                    <  2:1045393/147‑1 (MQ=255)                                                                                                                 tcAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                  >  1:903833/1‑151 (MQ=255)                                                                                                                   aaTAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                >  2:355413/1‑151 (MQ=255)                                                                                                                   aaTAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAAAATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                <  2:37141/151‑1 (MQ=255)                                                                                                                    aataattgggTAACCTAAGAGATTGGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                <  2:788793/143‑1 (MQ=255)                                                                                                                          gCGTTAATTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                         <  2:397832/151‑1 (MQ=255)                                                                                                                            ggTAATTT‑AGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                        <  1:410627/150‑1 (MQ=255)                                                                                                                               aattt‑AGAAATTAGGTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTc                     <  1:1006270/147‑1 (MQ=255)                                                                                                                                     agagATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg               >  2:317580/1‑151 (MQ=255)                                                                                                                                      aaagtaatgATTGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                 <  1:401596/140‑1 (MQ=255)                                                                                                                                           tAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac         >  1:729722/1‑151 (MQ=255)                                                                                                                                            aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac         >  1:328179/1‑150 (MQ=255)                                                                                                                                             ggATTGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgca             <  2:864713/145‑1 (MQ=255)                                                                                                                                               aTTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc     >  2:764513/1‑151 (MQ=255)                                                                                                                                                  ggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc  <  2:237097/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGC  >  NC_000913/5012‑5304

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