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.3%	106.7 / 11.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 (30/4);  new base C (0/27);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592                                                                                                                                                       |                                                                                                                                                      tCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAc                                                                                                                                                      >  1:385803/1‑151 (MQ=255)   aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                    >  1:696687/1‑151 (MQ=255)     ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                                   >  2:29104/1‑150 (MQ=255)     ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccg                                                                                                                                                  <  2:816427/151‑1 (MQ=255)         aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                               >  1:835473/1‑150 (MQ=255)                      aaacccTGAAAGGGCACACCCTGCCGTTATTTTTCAAAAAAAGACCCTTCT‑TTTGCGTCCTTTTATCGGAATAAAGAAAAAAAAAATAAACACCGGCAAGGATTTTTTTAAAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCAc                                                                                                                                  <  2:401846/146‑1 (MQ=255)                        acccTGAAAGGCCACACCCTGCCGTTTTTATTCGTAAAAAAACCCTTCT‑TTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTAAGGATTTTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTa                                                                                                                                <  1:1062371/148‑1 (MQ=255)                        aaCCTGCAATGACACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  1:443275/1‑149 (MQ=255)                               agtgacacaCCCTGCCTCTATTATGCGTATAAA‑AACCTTTGTTTTGCTTGCTTTTAGCGAAATAAAGAAAAAAAAAATAAACACCCGTGAGGTTATTTTTAGAAGTGTGCCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                        <  1:1086538/146‑1 (MQ=255)                               aaggccacaCCCTCCCGTTATTATTCTTAAAAAGACCCTTCG‑TTTGCGTACGTTTAGCGAAATAAAGAAAAAAAAAATAAACACCCGACAGGATATTGTAATAAATGGGGCAATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                        <  1:408661/146‑1 (MQ=255)                               aaTGCCACACCCTGCCGCTATTATTCGTAAAAAGAACCTTCT‑TTTGCGTACTTTTATCGAAATAAAGGAAGAAAAAATAAACACCCGTAAGGATATTGTTATCAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             <  1:848045/146‑1 (MQ=255)                               aaTGACACACCCTCCCGTTATTATTCATATAAAGAACCTTCT‑TTTGCGTGCTTTTATCGAAATTAAGGAAGAAAAAATAAACACCCGTCAGAATTTTTTTATAATTTTGCCCATCCCCCACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            <  1:976938/147‑1 (MQ=255)                                        ccctccctttttttttCGTAAAAAAACCCTTCT‑TTTGCGTTGTTTTAGCGAAATAAAGAAAAAAAAAATAAACACCCGTAAGGTTATTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  2:1092450/137‑1 (MQ=255)                                             ccGTTATTATTCTTAAAAAGACCCTTCT‑TTTGCGTCCTTTTATGGAAATAAAGAAAAAAAAAATAACCACCCGTCAGGATATTTTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcgc                                                                                                           <  1:203207/150‑1 (MQ=255)                                              cGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                         >  1:457102/1‑151 (MQ=255)                                                 tattatTCTTAAAAAAACCCTTCT‑TTTGCGTACTTTTATCGAAATTAAGGAAAAAAAAATAAACACCCGTCAGGATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                                                                                       <  1:1231409/150‑1 (MQ=255)                                                 tattatTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      >  2:89173/1‑151 (MQ=255)                                                     attCGTATTAAGAACCTTCG‑TGTGCGTACGTTTATCGGAATTAAGAAAGAAGCAATTAACAACCGTAACGATATTGTTATCAATGGGAAGATCCAATACCCCCCCGATAAAGCACGAACTAGCGTAGATATCAACGATGCGCtaa                                                                                                        >  2:1063483/1‑143 (MQ=255)                                                      ttaataaaaaGACCCTTTT‑TTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATTAACACCCGGCAGGATTTTGTAATAAATGTGACTACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                   <  2:360542/142‑1 (MQ=255)                                                      ttCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                                 >  2:1029604/1‑151 (MQ=255)                                                           aaaaaGACCCTTCT‑TTTGCTAACTTTTATCGAAATAAAGGAAAAAAAATTAAACACCCGACAGGATATTTTTATCAATGTGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                             <  1:1139064/148‑1 (MQ=255)                                                             aacaacctttCT‑TTTGCGTCTTTTTATCGGAATAAGGGAAAAAAAAATAAACCCCCGAAAGTATATTTTTATCAAAGGCGCCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           <  1:182826/142‑1 (MQ=255)                                                             aaaGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  1:102342/1‑150 (MQ=255)                                                             aaaGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  1:1118905/1‑150 (MQ=255)                                                              aaGACCCTTCG‑TTTGCGTATTTTTATGGAAATAAAGAAAAAAAAAATTAACACCCGTAAGGATATTGTTAAAAATGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                         <  1:873327/151‑1 (MQ=255)                                                               aGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        >  1:919050/1‑151 (MQ=255)                                                                       cG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 >  2:41655/1‑150 (MQ=255)                                                                          tttGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAAACAACCGCCAGGATTTTGTTATGAATGGGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                    <  1:437533/145‑1 (MQ=255)                                                                            tGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                            >  2:521189/1‑151 (MQ=255)                                                                             gCGTTCGTTTAGAGAAATAAAGAAAGAAAAAATAAACACCCGTCAGTATATTTTTAAAAATGTAACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                            <  2:203126/150‑1 (MQ=255)                                                                                tACTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATTTTTTTATAAATGGGACGTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  1:90662/149‑1 (MQ=255)                                                                                       aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  1:151930/1‑149 (MQ=255)                                                                                         cGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  2:1089769/1‑149 (MQ=255)                                                                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  2:329752/1‑148 (MQ=255)                                                                                          ggAATGAAGAAAAAAAAAATAAACACCCGACAGAATTTTTTTATCAGTGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              <  1:940929/151‑1 (MQ=255)                                                                                            aaTAAAGAAAAAAAAAATAAACACCGGAAAGAATTTTTTAATAAATGTGCCCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             <  2:927471/150‑1 (MQ=255)                                                                                            aaTAAAGAAAAAAAAAATAAACACCCGTAAGTTTTTTTTTATAAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             <  2:640179/150‑1 (MQ=255)                                                                                                    aaaaaaaaaTAAACACCCGAAAGGATTTTTTTAAAAATGTAACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                          <  2:233470/139‑1 (MQ=255)                                                                                                     aaaaaaaaTAAACACCCGACAGGATTTTGTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                   <  1:102818/146‑1 (MQ=255)                                                                                                        agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                >  1:972339/1‑151 (MQ=255)                                                                                                         aaaaTAACCACCCGTAAGGTTTTTTTAATAAATGTCACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:481402/150‑1 (MQ=255)                                                                                                           aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                             >  1:417270/1‑151 (MQ=255)                                                                                                           aaTAAACACCCGTAAGGTTATTTTAATCAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                             <  2:370222/151‑1 (MQ=255)                                                                                                              tAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                          >  1:761088/1‑151 (MQ=255)                                                                                                               aacaacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         >  1:230763/1‑151 (MQ=255)                                                                                                                  acccGTAAGGATATTGTTATAAATGTGACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                       <  2:514444/148‑1 (MQ=255)                                                                                                                  aacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                       >  2:640719/1‑150 (MQ=255)                                                                                                                    ccGTAAGGTTTTTTTTATCAATGTGCCGCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     <  1:119543/150‑1 (MQ=255)                                                                                                                        cAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                                 >  2:34075/1‑150 (MQ=255)                                                                                                                          ggATATTTTTATAAATGTACCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                                 <  2:1192205/148‑1 (MQ=255)                                                                                                                               ttGTTATAAATGTAACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                          <  1:1136701/150‑1 (MQ=255)                                                                                                                                tGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                        >  2:1249930/1‑151 (MQ=255)                                                                                                                                       attgtgctgATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                   <  2:1256774/141‑1 (MQ=255)                                                                                                                                         tgtgACCACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGa                <  1:564384/150‑1 (MQ=255)                                                                                                                                         tgagacgacCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGa                <  1:1114100/147‑1 (MQ=255)                                                                                                                                             aCGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg           >  1:966037/1‑151 (MQ=255)                                                                                                                                               gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  1:151284/1‑150 (MQ=255)                                                                                                                                                   cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGACTTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  1:562332/1‑151 (MQ=255)                                                                                                                                                   cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:955144/1‑151 (MQ=255)                                                                                                                                                     cTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGa    >  1:220144/1‑150 (MQ=255)                                                                                                                                                      tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:989601/1‑151 (MQ=255)                                                                                                                                                       aCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt      >  1:79778/1‑146 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCG‑TTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592

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