Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	44.1%	99.8 / 21.2	61	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (31/2);  new base C (0/26);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 9.54e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421296‑2421592                                                                                                                                                         |                                                                                                                                                      aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTAccc                                                                                                                                                    >  1:262155/1‑151 (MQ=255)   ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTAcccc                                                                                                                                                   >  2:1011821/1‑150 (MQ=255)     tAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCccgcc                                                                                                                                                >  2:985322/1‑151 (MQ=255)      aaaaaaaaCAACGGAGACCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCccgccg                                                                                                                                               <  2:886450/147‑1 (MQ=255)       aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCccgccg                                                                                                                                               >  1:188834/1‑150 (MQ=255)          aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATa                                                                                                                                            >  1:41681/1‑150 (MQ=255)             acaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAgc                                                                                                                                        >  2:210443/1‑151 (MQ=255)                   gagaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCAc                                                                                                                                  >  2:29564/1‑151 (MQ=255)                             atggccacaCCCTCCCTTTATTATTCGTAAAAAAACCCTTCGTTGGCTTACTTTTATCTAAATAAAGAAAAAAAAAATAAACCCCCGAC‑GGA‑ATCTTTTTATCAAT‑GTA‑AGG‑CCCCCCTCCCCCGCCGATAAAGCGCGCAc                                                                                                                                  <  2:561448/136‑1 (MQ=255)                                  cacaACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                    >  1:138799/1‑150 (MQ=255)                                    cacccggCCGCTATTTTTCGTAAAAAGACCCTTCGTTTGCGTACTTTTAGCGAAATTAGGAAAAAAAAAATAAACACCCGTC‑GGG‑ATATTGTTATAAAT‑GTG‑ACG‑TTCCACTTCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                                   <  1:984940/143‑1 (MQ=255)                                    cacccTGCCGTTATTTTTCGTAAAAAAACCCTTCGTTTGCGTACGTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGA‑ATTTTTTTATAAAT‑GGG‑ACA‑ATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                    <  1:896053/145‑1 (MQ=255)                                    caACCTGCTTCTTTTTTTCTTAAAAAAACCCTTCTTTGGCGTTGTTTTATCGAATTGAAGAAAAAAGAATTTAACACCCG‑C‑AGGAATATTTTTACCATT‑TTG‑ACG‑CCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  <  1:667828/150‑1 (MQ=255)                                    acaCCTGCCGTTATTTTTTGTAAAAAAACCCGTGGTTTGCTTAGTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGAC‑AGG‑TTATTGTTAAAAGG‑GTG‑ACT‑ACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                                   <  1:800475/147‑1 (MQ=255)                                     accctcccttttttattCATAAAAAAACCCTTCTTTTGGGTACTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTA‑AGG‑ATTTTGTTATAAAT‑GAG‑ACA‑ACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                                                                                       <  2:103133/132‑1 (MQ=255)                                     aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  2:1138400/1‑150 (MQ=255)                                       ccggcctttattttTCGAAAAAAAACCCTTTGTTTGCTTCCTTTTATCGAAATAAAGGAAAAAAAAATAAACACCCGTA‑AGA‑ATTTTGTTATAAGT‑GTG‑ACG‑ACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  1:633337/142‑1 (MQ=255)                                       ccTGCTGTTATTATTGTTAAAAATCCCCTTGGTTTGAGTCATTTTATCGAAATAAAGAAAAAAGAAATTAA‑AACCCGA‑GGG‑TTTTTTTTAATAATGGTG‑ACG‑ACCCCTTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                                                                              <  2:708849/151‑1 (MQ=255)                                           ccGTTATTATTATTAAAAAAACCCTTCTTTTGCGTAGGTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGGC‑AGG‑ATATTGTTAAAAAG‑GTG‑ACA‑CCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            <  1:18736/149‑1 (MQ=255)                                                 tttttCTTAAAAAAACCCTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTA‑AGA‑ATATTTTTATCAAT‑GTG‑CCA‑ACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            <  1:649113/140‑1 (MQ=255)                                                 ttattCTAAAAAAGATCTTTCGTTTGCTTTTGTTTATCGGAATAAAGAAAGAAAAAATTACCAGCCGAA‑AGG‑ATTTTGTAATAAAT‑GGG‑ACA‑ACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                     <  2:11578/150‑1 (MQ=255)                                                      cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                                >  2:980953/1‑150 (MQ=255)                                                         ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            >  2:1015337/1‑151 (MQ=255)                                                         aaaaaaacccTTCTTTTGCTTCCTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGA‑ATTTTGTTATAATT‑GTG‑ACG‑CTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            <  1:341014/143‑1 (MQ=255)                                                           aaaaacctttCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGCC‑AGG‑ATTTTTTAAAAAAT‑GGG‑ACG‑CCCCCCTGCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          <  1:443680/143‑1 (MQ=255)                                                            aaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          >  1:829552/1‑150 (MQ=255)                                                               accctttttttGCTTTCGTTTATAGAAATTAAGAAAAAAAAAAAAAACCCCCGTA‑AGG‑ATTTTTTTAAAAAT‑GAG‑ACG‑CCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      <  2:31089/143‑1 (MQ=255)                                                               aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                       >  2:148152/1‑150 (MQ=255)                                                                 cctttttttGCGTCGTTTTATCGAAATAAAGGAAAAAAAAATAAACACCCTTA‑AGG‑ATTTTTTTATCAAT‑GTG‑CGG‑ACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa                                                                                     <  1:489520/144‑1 (MQ=255)                                                                   ttCTTTTGCGTACTTTTATGGAAATTAAGAAAAAAAAAATAAACACCCGCC‑AGA‑ATATTTTAATAAAT‑GTG‑CCG‑ACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                   <  1:296377/150‑1 (MQ=255)                                                                     cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 >  1:947666/1‑150 (MQ=255)                                                                      gTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               >  1:261339/1‑151 (MQ=255)                                                                      gTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               >  2:748706/1‑151 (MQ=255)                                                                      gTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               >  2:869573/1‑151 (MQ=255)                                                                       tttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              >  2:552250/1‑151 (MQ=255)                                                                       tttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              >  1:546200/1‑151 (MQ=255)                                                                         tGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                            >  1:896797/1‑151 (MQ=255)                                                                             tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  2:912268/1‑149 (MQ=255)                                                                                tttttcgggAATTAGGAAAGAA‑AAATTAACCACCGGCGGGG‑TTATGGTAATCAAT‑‑TG‑ACGCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                        <  2:928554/141‑1 (MQ=255)                                                                                 tttATCGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGG‑ATTTTTTTATAAAT‑GTG‑ACG‑ACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                     <  1:237287/150‑1 (MQ=255)                                                                                 tttATAGAAATAAAGAAAAAAAAAATTAACAACCTTC‑AGA‑TTTTTTTTAAAAAT‑GTG‑ACG‑TCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                     <  2:137632/150‑1 (MQ=255)                                                                                    aTCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  1:444604/1‑149 (MQ=255)                                                                                       ggAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  1:951225/1‑148 (MQ=255)                                                                                         aaTAAAGAAAAAAAAAATAAACACCCGTC‑AGA‑ATTTTTTTATAAAT‑GTG‑ACA‑CCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                            <  1:665639/151‑1 (MQ=255)                                                                                         aaTAAAGAAAAAAAAAATAAACACCCGTA‑AGG‑ATTTTTTTATAAAT‑GTG‑ACG‑ACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 <  1:347877/146‑1 (MQ=255)                                                                                             aagaaaaaaaaatatAACACCCGTA‑AGT‑ATTTTTTTATAAAT‑GAG‑ACA‑CCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                          <  2:798380/135‑1 (MQ=255)                                                                                             aaGGAAAAAAAAATAAACACCCGTC‑AGG‑ATTTTTTTATAAAT‑GTG‑ACT‑CCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                          <  1:124744/149‑1 (MQ=255)                                                                                             aaGAAAAAAAAAATAACCACCCGTA‑AGA‑ATTTTGTTATAAGT‑GAG‑CCG‑ACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         <  2:237014/150‑1 (MQ=255)                                                                                                     aaaaaTAAACACCCGTC‑AGG‑ATTTTTTTATCAAT‑GGG‑ACG‑CCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                <  1:1015474/149‑1 (MQ=255)                                                                                                      gaaATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               >  1:65093/1‑151 (MQ=255)                                                                                                      aaaaTTAACACCCGTC‑AGA‑ATTTTTTTACAAAT‑GTG‑ACA‑CCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                <  2:994529/149‑1 (MQ=255)                                                                                                            accacccGAC‑AGA‑ATTTTGTTAT‑AAT‑GGGAACG‑ACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         <  2:822537/146‑1 (MQ=255)                                                                                                              caacCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                        >  1:450088/1‑150 (MQ=255)                                                                                                                  cGTC‑AGG‑ATATTGTTATCAAT‑GGG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCGGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                                   >  2:75101/1‑151 (MQ=255)                                                                                                                     c‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                                 >  1:168875/1‑150 (MQ=255)                                                                                                                                 ttATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                      >  1:1017816/1‑151 (MQ=255)                                                                                                                                  taacAAT‑GTG‑ACG‑CCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                      <  2:357129/147‑1 (MQ=255)                                                                                                                                    tCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGGTATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                   >  1:486968/1‑151 (MQ=255)                                                                                                                                      aaT‑GTG‑CCG‑ACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                   <  2:775053/149‑1 (MQ=255)                                                                                                                                                g‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  2:714548/1‑150 (MQ=255)                                                                                                                                                        tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  2:22143/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                      ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGG‑ATATTGTTATCAAT‑GTG‑ACG‑ATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421296‑2421592

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