Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	41.9%	125.8 / 25.3	74	L79V (TTA→GTA)	maeB	malic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (40/3);  new base C (0/31);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 8.99e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577999‑2578291                                                                                                                                                      |                                                                                                                                                        cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑cccc                                                                                                                                                    >  1:1060311/1‑149 (MQ=255) cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑cccc                                                                                                                                                    >  2:747695/1‑149 (MQ=255) cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑cccc                                                                                                                                                    >  2:449814/1‑149 (MQ=255) cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑cccc                                                                                                                                                    >  2:1105093/1‑149 (MQ=255)          cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGt                                                                                                                                       >  2:259046/1‑151 (MQ=255)            aCAACTTCAATAAATTTGCCCGGGTCGAGTTCGCCCAATTCAATGTCACATACATCAATCTCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTCTGCCT‑GCCAGCG‑CGTCAATGTTG‑CCTA‑A‑CCCCA‑GC‑CACGCTGTa                                                                                                                                      >  2:992673/1‑150 (MQ=255)              aaCTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCg                                                                                                                                   >  2:390468/1‑151 (MQ=255)               aCTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGt                                                                                                                                  >  1:252162/1‑151 (MQ=255)                   cAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGt                                                                                                                                  >  2:643297/1‑147 (MQ=255)                    aaaaattttGTCGGGGTCGATTTCGAAAACTTAAATGACAAATACAACAACCCCGCCGAATTTTTTAAACAGACCCCCCTTC‑CCTCCCATCCCCGTTTTGCCG‑GCAGGCG‑CGCAAGTGTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTag                                                                                                                               <  2:1269779/143‑1 (MQ=255)                          tttGTCCGGGTCGAGTTCGTCATCTTCAATCAAAAATCCATCATCCCGGGCGAATTTCTTAAACGGAACCCCCTTG‑CCATCCACCCCCGGTTTGCCT‑CCCCGGG‑CGCTAATGTTG‑CCTC‑C‑CCCCA‑GC‑TCCGCCGTACCGTTagag                                                                                                                             <  2:624261/145‑1 (MQ=255)                          tttGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCAcc                                                                                                                       >  2:1206683/1‑151 (MQ=255)                          tttGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCAcc                                                                                                                       >  2:737448/1‑151 (MQ=255)                              tCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                    >  1:290788/1‑150 (MQ=255)                               cggggTCGGTTTCGTAAACTTAAGTATAAAATCCACCAACCCCGGAAAATTT‑TTAAAAAAAACCCCCTTGCCCTCCCATCCCGGTTTTCCCT‑GCGGGCG‑CCCCAGTTTTG‑CCTA‑C‑CCCCA‑GC‑CCCGCCGTACCGTTAGAGATCACCGcca                                                                                                                   <  1:1219778/148‑1 (MQ=255)                               ccGGGTCGATTCCGTCAACTAAAATGTCAAATCCACCAACCCCGGCGAATTTCTAAAACAGACCCCCCTTC‑CCTCCCACCCCCGTTTTGCCT‑GCCAGCG‑CCCAAATGTTG‑CCAA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGc                                                                                                                     <  2:618996/148‑1 (MQ=255)                                cGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGccacc                                                                                                                 >  1:310213/1‑151 (MQ=255)                                  gtgcgggTTCGTCAACTTAAATGTCAAATAAATCAACCCCGGCGATTTTCTTAAAAAGACCCCCCTTC‑CCTTCCACCCCGGTTTTGCCG‑GCCAGCG‑CGCAATTTTTG‑CCAA‑C‑CCCCA‑CC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAg                                                                                                               <  1:983251/145‑1 (MQ=255)                                  ggTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGccacc                                                                                                                 >  2:1122955/1‑149 (MQ=255)                                    tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                             >  1:202020/1‑151 (MQ=255)                                    tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                             >  2:865790/1‑151 (MQ=255)                                      gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                            >  2:772103/1‑150 (MQ=255)                                      gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                            >  2:89752/1‑150 (MQ=255)                                           cgccaccTTCAATCAAAAACACATCAATCCCGGCGATTTTCTAAAACAGACCGCCCTTC‑CTTTCCCCCACCGGTTTGCTT‑GCAAGCG‑GGCCAAGGTGG‑CTTA‑C‑CCCCA‑GA‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                           <  1:169879/140‑1 (MQ=255)                                              cAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                                     >  1:683480/1‑149 (MQ=255)                                              cAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                   >  2:745076/1‑151 (MQ=255)                                               aaCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                                  >  1:1291244/1‑151 (MQ=255)                                               aaCTTAATTGTAAAATACATCAACCCCGGGAAATTTCTAAAAAAGACCGCCCTTC‑CCTCCCACCCCCGGTTTGCCT‑GCCGGCG‑CGCCAATGTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                                     <  1:522995/148‑1 (MQ=255)                                                aCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg                                                                                                 >  1:1213411/1‑151 (MQ=255)                                                  ttaaaTGTAAAATACATCAACCCCGGCGATTTTCTAAAACAGACCGCCTTTG‑CCTCCCACCCCCGGTTTGCCT‑GCCGGCG‑CGCCAATGTTG‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                                    <  1:280449/143‑1 (MQ=255)                                                     atgtaaaaaCACACCAACCGCGGTGAGTTTCTAAAAAAGAGCCCCCTTG‑CTTCCCCCCCCCGTTTTCCCC‑GCGCGCG‑ACCAGTTGTTG‑CC‑A‑A‑CCCCC‑GCAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtata                                                                                             <  2:1354477/144‑1 (MQ=255)                                                     aaTGTAAAATACATCAACCCCGGCAATTTTCTAAAACAGACCCCCCTTG‑CCTCCCATCCCCGGTTTCCCT‑CCCAGCG‑CGCAATTTTTG‑CCAA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtata                                                                                             <  1:1244162/150‑1 (MQ=255)                                                           aaaTACATCACCCCCGGCAATTTTCTAAAACAGAACGCCTTTG‑CCTTCCCTCCCCGGTTGGCCT‑GCCAGCG‑CGCCAATGTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                          <  2:1068837/147‑1 (MQ=255)                                                           aaaTACATAAACCCCGGCAATTTTCTAAAAAAGACCCCCCTTC‑CCTCCCACCCCGGTTTTGCCG‑GCCAGCG‑CGCAAATGTTG‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg                                                                                         <  1:1295362/148‑1 (MQ=255)                                                                 atcttccccgcggAATTTCTTAACCACACCG‑CCTTT‑CCCTCCCCCACCAGTTTTCCC‑CCCAGGG‑CCCCAATCGTC‑ACTA‑C‑CCCCC‑CCGCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                       <  1:907807/132‑1 (MQ=255)                                                                 aTCAACCCCGGCGATTTTCTAAAACAGACCCCCTTTC‑CCTCCCACCCCCGTTTTCCCG‑CCCGGCG‑CCCCATTTTTC‑CCAA‑C‑CCCCG‑GC‑CCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                 <  2:1299606/150‑1 (MQ=255)                                                                  tCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                                >  1:784335/1‑150 (MQ=255)                                                                   cAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTa                                                                               >  2:489379/1‑150 (MQ=255)                                                                    aaccccGGCAATTTTTTAAAACAGACCCCCTTTG‑CCTTCCTTCCCCGTTTTCCCT‑GCCGGCG‑CGCCATTGTTC‑CCTA‑C‑CCCCG‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTaa                                                                              <  1:831101/147‑1 (MQ=255)                                                                     aTCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                            >  1:905336/1‑151 (MQ=255)                                                                      ccccGGCAATTTTCTAAAACAGACCGCCCTTG‑CCTCCCACCCCCGGTTTGCCT‑GCCAGCG‑CGCAAATGTTG‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                           <  1:200819/150‑1 (MQ=255)                                                                       cccGGCAATTTTCTAAAACAGACCCCCCTTG‑CCTTCCCCCCCGGGTTGGCCT‑CCCCGCG‑CCCAATTGTGC‑CCTA‑C‑CCCAA‑CC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                            <  1:589710/149‑1 (MQ=255)                                                                          ggCGATTTTCTAAAACACACCGCCCTTC‑CCTCCCATCCCCGTTTTGCTG‑CCCGGCC‑CCCAAATTTTC‑CCAA‑C‑CCCCA‑CC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                                        <  1:913830/150‑1 (MQ=255)                                                                          ggCAAATTCCTTAAAAAGAGCCCCCTTG‑CCTCCCACCCCCGGTTGGCCC‑CCCAGCG‑CGCCATTGTTG‑CCAA‑C‑CCCCA‑GA‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                                        <  1:85054/150‑1 (MQ=255)                                                                            cGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt                                                                     >  1:240796/1‑151 (MQ=255)                                                                              aaTTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                   >  1:781454/1‑151 (MQ=255)                                                                              aaTTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                   >  1:371300/1‑151 (MQ=255)                                                                                tttCTAAAACAAACCGCCCTCG‑CCTCCCA‑CCCCGGTTTTCCC‑GCCCGCC‑GCCCAATTTTGCCCTA‑C‑CCCCA‑GC‑TCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                   <  2:376106/149‑1 (MQ=255)                                                                                  cttaaaacaaacCGCCCTTC‑CTTCCCCCCCCCGGTTCCCCT‑GCCACCG‑CCCCAGTGTTC‑CCAA‑C‑CCCCA‑GC‑CCCGCCGTACCGTTAGAGATCACCGCCTCCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                  <  2:33929/144‑1 (MQ=255)                                                                                    ttAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                             >  1:1205473/1‑151 (MQ=255)                                                                                      aacaagacggCCTTTC‑CCTTCCA‑CCCCGGTTTGCCTGGCCGGCG‑C‑CCAATTTTC‑CCCACC‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                             <  2:484300/140‑1 (MQ=255)                                                                                      aaaCAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTC‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                            >  1:750413/1‑150 (MQ=255)                                                                                       aaCAGACCGCCTTTC‑CCTCCCACCCCCGTTTGGCCT‑GCCAGCC‑CCCAAGTGTTG‑CCAA‑C‑CCCCA‑GA‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCaa                                                          <  2:450414/151‑1 (MQ=255)                                                                                       aaCAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           >  1:650437/1‑150 (MQ=255)                                                                                         cAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           >  2:961581/1‑148 (MQ=255)                                                                                          aGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGAG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                       >  1:93066/1‑151 (MQ=255)                                                                                            accGCCCTTC‑CCTCCCCCCCCCGTTTTGCCG‑CCCAGCG‑CGCAATTTTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACaa                                                     <  1:277077/149‑1 (MQ=255)                                                                                                 cctttcccTTCCATCCCCGTTTTGCCT‑GCCAGCG‑CGCAATTGTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                  <  1:1288746/144‑1 (MQ=255)                                                                                                         ttCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                         >  2:1345944/1‑151 (MQ=255)                                                                                                         tcccATCCCCGTTTTGCCT‑CCCAGCG‑CCCAATTGTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                               <  1:1051621/143‑1 (MQ=255)                                                                                                          tCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGc                                        >  2:291966/1‑151 (MQ=255)                                                                                                           ccacccccGGTTTGCCG‑CCCGGCG‑CGCCATTGTTG‑CCCA‑C‑CCCCA‑GC‑CCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                         <  2:253629/143‑1 (MQ=255)                                                                                                               ccctggTTTGCTT‑GCCGGCG‑CCCAATTTTTC‑CCTA‑C‑CCCCA‑CC‑CCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTg                                     <  2:527068/145‑1 (MQ=255)                                                                                                                aCCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                  >  1:444595/1‑151 (MQ=255)                                                                                                                 cggttttGCCT‑‑CCCGCGCCGCCAATGTTG‑CC‑A‑A‑CCCCAGGC‑CCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                  <  1:1194988/146‑1 (MQ=255)                                                                                                                      ttttcgccccggcg‑cgcAAGTGTTC‑CCCA‑C‑CCCCG‑AC‑GCCGCCGTACCGTTAGAGATCGCCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                                      <  2:1014900/131‑1 (MQ=255)                                                                                                                      ttccct‑gccAGCG‑CGCCATTGTTC‑CCTA‑C‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAg                             <  1:737460/147‑1 (MQ=255)                                                                                                                              ccAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc                     >  2:561276/1‑151 (MQ=255)                                                                                                                              ccAGCC‑CGCAAATTTCC‑CCTA‑C‑CCCAA‑GC‑CCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTa                              <  2:922188/142‑1 (MQ=255)                                                                                                                               cAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                    >  2:784335/1‑151 (MQ=255)                                                                                                                                   gcgccccagtttG‑CCCA‑ACCCCCA‑AC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCa                   <  2:253664/140‑1 (MQ=255)                                                                                                                                   cccccaattGTTC‑CCTA‑C‑CCCCA‑GC‑CCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAg                        <  1:459367/136‑1 (MQ=255)                                                                                                                                        cAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg            >  2:92040/1‑151 (MQ=255)                                                                                                                                             ttG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtgtg       >  2:220402/1‑151 (MQ=255)                                                                                                                                                    a‑a‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGCGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc  >  1:461544/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                        CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCTTCCATCACCGGTTTGCCT‑GCCAGCG‑CGCCAATGTTG‑CCTA‑A‑CCCCA‑GC‑ACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577999‑2578291

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