Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,437	0	T	C	44.6%	87.5 / 20.3	65	I38T (ATC→ACC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (26/10);  new base C (0/29);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 3.29e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACG  >  NC_000913/2421288‑2421587                                                                                                                                                        |                                                                                                                                                       attattTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATc                                                                                                                                                       >  2:51168/1‑151 (MQ=255)  ttattTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATcc                                                                                                                                                      >  1:125921/1‑151 (MQ=255)   tattTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATcc                                                                                                                                                      <  1:791218/150‑1 (MQ=255)    attTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCAc                                                                                                                                                    >  2:678447/1‑151 (MQ=255)      ttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACt                                                                                                                                                   <  1:801161/150‑1 (MQ=255)      ttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACt                                                                                                                                                   <  2:709731/150‑1 (MQ=255)       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)                                                                                                                                                   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)                                                                                                                                                        |                                                                                                                                                       ATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCG‑GAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACG  >  NC_000913/2421288‑2421587

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