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

TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592                                                                                                                                                        |                                                                                                                                                      tCATCCTAAGTAAAACAACGGAGACCCTGCAATGGCACAAGCTGCCGCTATTATTCGGATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCCCTa                                                                                                                                                       <  2:472874/150‑1 (MQ=255)     ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCccg                                                                                                                                                  <  1:98435/151‑1 (MQ=255)         aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGa                                                                                                                                              <  1:1131593/151‑1 (MQ=255)                       gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  1:730399/1‑150 (MQ=255)                        acccgaaaaggccacaCCCTGCCGTTTTTATTCATAAAAAAACCCTTCTTTTGCGTACTTTTATGG‑AAATAAAGAAAAAAAAAATAAACACCCGTCAGGA‑TTTTGTTATAAAGGTGCCGCCCCACTCCCCCGCCGATAAAGCGCGCACTAg                                                                                                                               <  1:1136742/139‑1 (MQ=255)                         aCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                              >  2:877002/1‑151 (MQ=255)                           cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             >  2:1180756/1‑150 (MQ=255)                              aaaggccacccccTGCCTTTTTTTTTCATAAAAAAACCCTTCTTTTGCTACCTTTTACAG‑AAATTAAGAAAAAAAAAATAAAAACCCGAAAGGA‑TTTTTTTAAAAATGTGACGTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            <  2:34803/137‑1 (MQ=255)                                 ggccacaCCCTGCCGTTATTATTCATAAAAAAACCCTTCTTTGGCGTATGTTTATCG‑AAATAAAGGAAAAAAAAATTAACACCCGTAAGGA‑TTTTTTTATCAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                      <  1:643679/148‑1 (MQ=255)                                       aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  1:1090785/1‑150 (MQ=255)                                           tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            >  2:79405/1‑151 (MQ=255)                                           tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            >  2:79378/1‑151 (MQ=255)                                             ccttttttattCTTAAAAAAACCCTTCTTTTGCGTACGTTTATCG‑AAATAAAGAAAAAAAAAATAAACACCCGCCAGAA‑TTTTGTTAACAATGTGACGCCCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             <  1:609264/142‑1 (MQ=255)                                             ccGTTATTATTCTAAAAAAGACCCTTCGTTTGGTTACGTTTATCG‑AAATAAAGGAAGAAAAAATAAACACCCGGCAGAA‑TATTGTTATAAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                          <  2:223403/151‑1 (MQ=255)                                                 tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGCTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      >  1:280787/1‑151 (MQ=255)                                                   ttattCTAAAAAAGACCCTTTTTTTGCGTACTTTTATGG‑AAATAAAGGAAAAAAAAATTACCACCCGT‑AGGATTTTTGTTATAAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      <  1:685597/149‑1 (MQ=255)                                                        cGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                                >  2:826165/1‑150 (MQ=255)                                                        cGTAAAAAGACCCTTCGTTTGCGTCCGTTTATGG‑AAATTAAGAAAAAAAAAATAAACAACCGTCAGTA‑TTTTGTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:1147537/151‑1 (MQ=255)                                                           ataAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            >  1:669556/1‑151 (MQ=255)                                                           aaaaaaacccTTCTTTTGCTTTTGTTTATGG‑AAATAAAGAAAAAAAAAATAAACACCCGTCAGGA‑TTTTTTTATAAATGTGCCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            <  1:653731/143‑1 (MQ=255)                                                             aaaaacccTTCTTTTGCGTACTTTTATGG‑AAATAAAGAAAAAAAAAATAAACACCCGTCAGGA‑TTTTTTTATAATTGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          <  1:393050/145‑1 (MQ=255)                                                             aaaaacccTTCGTTTGCGTATTTTTATGG‑AAATTAAGAAAAAAAAAATAAACACCCGTCAGGA‑TTTTTTTATAATTGGGAACACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          <  2:562877/145‑1 (MQ=255)                                                             aaaGACCCTTCTTTGGCGTACTTTTATCG‑AAATAAAGGAAAAAAAAATAAACACCCGTAAGGA‑TTTTTTTATAAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCATCGt                                                                                           <  2:803865/150‑1 (MQ=255)                                                             aaaGACCCTGCTTTGGCGTACTTTTATGG‑AAATAAAGAAAGAAAAAATAAACACCCGAAAGGA‑TTTTTTTATCATTGGGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:972214/146‑1 (MQ=255)                                                               acacccTTTGTTTGCGTAC‑TTTTTCGAAAATAAAGGAAAAAAAAATAAAC‑ACCGGTAGGATTTTTTTTATAATTGTGACTTCCCCCTCCGCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        <  1:682407/147‑1 (MQ=255)                                                               aaacctttCTTTGGCGTTCTTTTATGG‑AAATAAAGAAAAAAAAAATAAACACCCGTCAGGA‑TTTTGTTATAAATGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        <  1:892985/145‑1 (MQ=255)                                                                gAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                       >  2:668356/1‑151 (MQ=255)                                                                       cGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 >  1:584391/1‑150 (MQ=255)                                                                          tgtgtttcGTTTGTAG‑AAATGAAGGAAGAAGAATTCAACGACCCTGGGGT‑TATGGTCTGCAAGTTGATAATCGTCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa                                                                                  <  2:1060108/139‑1 (MQ=255)                                                                          tggCGTTCTTTAATCG‑AAATTAAGGAAAAAAATATTACCACCCGTAAGAA‑TATTGTTAAAAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              <  2:252813/148‑1 (MQ=255)                                                                           tGCGTCCTTTTATCG‑AAATAAAGAAAGAAAAAATAAACACCCGTAAGTT‑TTTTTTTATAAATGTGACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                             <  2:815520/150‑1 (MQ=255)                                                                             cGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  1:879120/1‑151 (MQ=255)                                                                              gTACTTTTATGG‑AAATAAAGAAAAAAAAAATAAACACCCGTCAGGA‑TTTTGTTATCAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  1:1130047/150‑1 (MQ=255)                                                                                  ttttATCG‑AAAAAAAGAAAAAAAAAAAAAACACCCGTCAGAT‑TTTTGTAATAATTGGGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  1:432131/145‑1 (MQ=255)                                                                                   tttATGG‑AAATAAAGAAAAAAAAAATAAACACCCGAAAGAT‑TTTTTTTAAAAATGTGACCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                     <  2:391166/150‑1 (MQ=255)                                                                                         g‑gAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  1:279622/1‑148 (MQ=255)                                                                                           gAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              >  1:803772/1‑150 (MQ=255)                                                                                            aaTTAAGAAAAAAAAAATAAACACCGGCAAGGT‑TTTTTTTAAAAATGTGACGCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              <  2:1067640/149‑1 (MQ=255)                                                                                            aaTAAAGAAAAAAAAAATAACCACCCGTCAGTT‑TTTTTTTAAAAAGGTGCCGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             <  2:402693/150‑1 (MQ=255)                                                                                                aaGGAAGAAAAAATTAACACCCGTCAGGA‑TTTTGTTATAAGTGTACCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                                  <  1:465138/141‑1 (MQ=255)                                                                                                aaGAAAAAAAAAATAACCACCCGTAAGGT‑TTTTGTTATAATGGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                          <  1:936045/149‑1 (MQ=255)                                                                                                aaGAAAAAAAAAATAAACACCCGTAAGAA‑TTTTTTAATAAATGTGCCAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAATATTATTCAGCATGTAGAGAATAAcc                                                        <  1:551835/151‑1 (MQ=255)                                                                                                       aaaaaaTAACCACCCGTAAGGA‑TTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                                 <  2:518697/148‑1 (MQ=255)                                                                                                         aaaataaccccccGTCAGTA‑TTTTTTTATAAATGTACCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:242389/140‑1 (MQ=255)                                                                                                               aacaccCTTAAGGA‑TTTTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                            <  2:172373/148‑1 (MQ=255)                                                                                                               aacaacCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         >  1:857308/1‑151 (MQ=255)                                                                                                               aaaacccttaaGGT‑TATTGTTATCATGGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         <  1:466498/141‑1 (MQ=255)                                                                                                                  aacCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                       >  2:1122506/1‑150 (MQ=255)                                                                                                                    ccGTAAGAA‑TATTTTTATCAATGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                    <  2:675985/151‑1 (MQ=255)                                                                                                                      gTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  >  2:60059/1‑151 (MQ=255)                                                                                                                       acAGGA‑TATTGTTAAAAATGTGACGATCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  <  1:741978/149‑1 (MQ=255)                                                                                                                        cAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                                 >  2:448553/1‑150 (MQ=255)                                                                                                                         aGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                               >  1:1058879/1‑151 (MQ=255)                                                                                                                              tatTTTTAACAATGTGACTCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCa                            <  2:416914/150‑1 (MQ=255)                                                                                                                                ttGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                          >  2:676843/1‑150 (MQ=255)                                                                                                                                   ttATAAAGGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCa                            <  1:466885/145‑1 (MQ=255)                                                                                                                                           gtgcctctcctccaCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGata              <  1:745555/138‑1 (MQ=255)                                                                                                                                             gACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc            >  1:45268/1‑151 (MQ=255)                                                                                                                                               cGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  2:531268/1‑151 (MQ=255)                                                                                                                                                gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  2:314905/1‑150 (MQ=255)                                                                                                                                                  tCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACg       >  1:179727/1‑151 (MQ=255)                                                                                                                                                   ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt      >  2:607598/1‑151 (MQ=255)                                                                                                                                                    cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:897418/1‑151 (MQ=255)                                                                                                                                                     aCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:484112/1‑150 (MQ=255)                                                                                                                                                       tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:288425/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑