Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	44.3%	74.5 / 43.1	61	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (34/0);  new base A (0/27);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 6.33e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGC  >  NC_000913/1755722‑1756017                                                                                                                                                        |                                                                                                                                                   aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGt                                                                                                                                                    >  2:702275/1‑150 (MQ=255)       ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAAccg                                                                                                                                             >  1:52386/1‑151 (MQ=255)        gACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAAccgc                                                                                                                                            >  1:509121/1‑151 (MQ=255)         aCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAAccgcc                                                                                                                                           >  1:163201/1‑151 (MQ=255)                  cGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTg                                                                                                                                  >  2:289274/1‑151 (MQ=255)                   gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGc                                                                                                                                 >  1:892609/1‑151 (MQ=255)                      tCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTg                                                                                                                              >  1:439151/1‑151 (MQ=255)                        tGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGAt                                                                                                                            >  1:554655/1‑151 (MQ=255)                          aaGAGTTTTTAGCAAAAATGCTGGCCGCTGGCATGAACTTTATGCTTCTAAACTCCCCTCATGGTGACTATGAAGAACAGGGCCAGCCCATCCAAATTCGGCGAAAGGTGAGGAG‑A‑AAAACTGGAAAAACCGCCGCTATCCTGCTTGa                                                                                                                             <  1:782964/148‑1 (MQ=255)                          aaGAGATTTTAGAAAAAATGCGGGACGCTGGCATGAACTTTATGCTTCTGATCTTCTCTCATGGTGACTAGGCAGAACACGGTCAGCGCATTAAGAATCTGCGCAACGTGATGAG‑A‑AAAACTGGAAAAACCGCCGCTATCCTGCTTGAt                                                                                                                            <  2:356330/149‑1 (MQ=255)                          aaGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATAc                                                                                                                          >  1:852604/1‑151 (MQ=255)                                  ttagaaaaaaTGCTGGACGCTGGCATGAACGTTATGCGTTTAATCTTCTCTCTTGGTGATTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑A‑AAACCTGGAAAAGCCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                    <  2:49849/143‑1 (MQ=255)                                           aTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                                         >  1:150217/1‑151 (MQ=255)                                            tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCg                                                                                                        >  2:209418/1‑151 (MQ=255)                                                ggACGCTGGCATGAACGTTATGCGTCTGACCTTCTCTCGTGGTGACTAGGCAGACCCGGGTCAGCGCATTAAGAATCTGCGCAACGTAATGAG‑A‑AAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                     <  2:250239/150‑1 (MQ=255)                                                ggACGCGGGCAGAACCGTTATGCGTCTGACCTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCACCGTGATGAG‑A‑AAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                     <  2:119696/150‑1 (MQ=255)                                                ggACGCGGGAAGGAACTTTATGCGTCAAACCTCCCCTCATGGTAACTATGCAGAACAGGGCCAGCCCATTCAGAATCTGCGCAACGTGAGGAG‑A‑AAAATTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCa                                                                                                      <  2:410757/149‑1 (MQ=255)                                                      tGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGaa                                                                                               >  1:38258/1‑150 (MQ=255)                                                       ggCATGAACGTTATGCGTCTAATTTTCTCTCGTGGTGATTAGGAAGAACGCGGCCAGCGAATTCAGATTCTGCGCAACGTGAGGAG‑A‑AAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATg                                                                                                 <  2:196883/147‑1 (MQ=255)                                                        gCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACt                                                                                            >  2:125105/1‑151 (MQ=255)                                                          atgacctttagccttCTGAACTTCTCTAGGGGTGACAATCAAAAACCCGACCAGCCCATACAAATTCTCCGAAACGTGCTGAG‑A‑AAAACTGAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGa                                                                                                <  1:63153/131‑1 (MQ=255)                                                          ataaaCGTTATGCTTCTAAACTTCTCTCATGGTAACTAAGCAAAACACGGCCACCCCATTCAGAATCTGCACAACGTAATGAA‑A‑AAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGa                                                                                                <  1:37434/142‑1 (MQ=255)                                                          aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                          >  2:575302/1‑151 (MQ=255)                                                            gAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                         >  1:547054/1‑150 (MQ=255)                                                             accGTTATGCGTCTGCCCTTCCCACAGGTTGAATATGAAAAACACGGACAGCGCATCCAAATTCGGCGCAACGGGAGAAG‑A‑AAAATGGGAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGaa                                                                                        <  1:462125/148‑1 (MQ=255)                                                                gTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                    >  1:792368/1‑151 (MQ=255)                                                                      cGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                              >  2:198155/1‑151 (MQ=255)                                                                        tttGAATTCCTCTCAGGGTGACTAGCCAAAACACGGCCAGCGCTTTCAAAGTGTGCACAAGGTGA‑GAGAA‑AAAACTGAAAAAGCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAac                                                                               <  1:691267/146‑1 (MQ=255)                                                                        tctaaccTTCTCTCAGGGTGACTATGCAAAACCCGGTCAGCGAATTCAAAATCTGCCAATCGTGATGAG‑A‑AAACCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCgg                                                                                   <  2:985794/138‑1 (MQ=255)                                                                         cTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgacg                                                                           >  2:19198/1‑151 (MQ=255)                                                                           gAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                         >  1:934978/1‑151 (MQ=255)                                                                            aaCTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                        >  1:475134/1‑151 (MQ=255)                                                                                   ctcATGGTGACAAGACAGACCAGGGCCAGCGCATAAAAATCCGGCCCAAGGTGATAAA‑A‑AAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                        <  1:764244/144‑1 (MQ=255)                                                                                     cATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                                 >  1:22197/1‑149 (MQ=255)                                                                                       tattaactatcaaaaaCACGTTCACCCCATCAAAAATCTGCGCAACGTGAGAAG‑A‑AAAATTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                    <  2:82464/130‑1 (MQ=255)                                                                                        ggTGACTAGAAAAAACACGGTCAGCGCATCAAGAATCGGCGAAACGTGATGAG‑A‑AAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAg                                                               <  1:209351/148‑1 (MQ=255)                                                                                            aCTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                         >  2:671129/1‑150 (MQ=255)                                                                                              tATGCAAAACACGGTCAGCGCATTCAGAACCGGCAAAACGTGATGAG‑A‑AAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                          <  1:848530/147‑1 (MQ=255)                                                                                               aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                     >  2:594308/1‑151 (MQ=255)                                                                                               aTGCAGAACACGGCCAGCGCATACAAAATCGGCGAAAGGTGAGGAA‑A‑AAACCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                       <  1:540433/149‑1 (MQ=255)                                                                                                   aaacccggctCAGCGCATTCAAAATCTGCGCAGGGTGAGGAG‑A‑AAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                   <  1:319396/140‑1 (MQ=255)                                                                                                    gAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                                 >  1:733063/1‑150 (MQ=255)                                                                                                        acGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACttt                                            >  2:78502/1‑151 (MQ=255)                                                                                                          ggccAGCGCATTCAGAACCCGCGCAA‑GGGAAGAA‑CAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                           <  1:263692/147‑1 (MQ=255)                                                                                                           gTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcac                                         >  1:741027/1‑151 (MQ=255)                                                                                                                cgaaTTAAGAATCTGCGAAACGTGATGAA‑A‑AAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACt                                     <  1:681540/147‑1 (MQ=255)                                                                                                                 gcATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                                    >  1:437116/1‑150 (MQ=255)                                                                                                                       aGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                             >  1:526894/1‑151 (MQ=255)                                                                                                                               cgcAACGTAATGAA‑A‑AAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATc                      <  1:702139/150‑1 (MQ=255)                                                                                                                                 cAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                   >  2:843014/1‑151 (MQ=255)                                                                                                                                  accGTGATAAG‑A‑AAAATTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                             <  1:919097/138‑1 (MQ=255)                                                                                                                                  aaCGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                   >  2:387081/1‑150 (MQ=255)                                                                                                                                    cGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc                >  2:129049/1‑151 (MQ=255)                                                                                                                                      ggatgaG‑A‑AAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCaa                 <  1:401747/147‑1 (MQ=255)                                                                                                                                             acaaaaaCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGt                          <  2:843926/131‑1 (MQ=255)                                                                                                                                              c‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg       >  1:866237/1‑151 (MQ=255)                                                                                                                                                aaaaCTGTAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa          <  2:38695/147‑1 (MQ=255)                                                                                                                                                aaaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg      >  1:973384/1‑151 (MQ=255)                                                                                                                                                aaaaCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg            <  2:682463/145‑1 (MQ=255)                                                                                                                                                aaaaCTGAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa           <  2:554768/146‑1 (MQ=255)                                                                                                                                                  aaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt    >  2:374847/1‑151 (MQ=255)                                                                                                                                                    cTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGc  >  2:814363/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                   ACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAG‑C‑AAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGC  >  NC_000913/1755722‑1756017

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