Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	52.6%	37.9 / 85.8	57	L79V (TTA→GTA)	maeB	malic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (22/5);  new base C (0/30);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 2.31e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-07

CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTT  >  NC_000913/2577999‑2578288                                                                                                                                                     |                                                                                                                                                   cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑Ag                                                                                                                                             <  1:786463/151‑1 (MQ=255)      gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCAccg                                                                                                                                        >  2:35891/1‑151 (MQ=255)      gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCAccg                                                                                                                                        >  1:552674/1‑151 (MQ=255)        ggcgTCAACTTCAATAAATTTGTCCGGGTCGAGGTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGCATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCAccgcc                                                                                                                                      >  1:849501/1‑151 (MQ=255)            aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTAc                                                                                                                                  <  1:771999/151‑1 (MQ=255)             cAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTAcc                                                                                                                                 >  1:676490/1‑151 (MQ=255)                  tCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGtt                                                                                                                              >  2:879608/1‑149 (MQ=255)                     aaaaaTTTGTCCGGGTCGGGTTCGTCAACTTAAATGTAAAATACACCAACCCCGGCGAATTTCTT‑AAACAGACCGCCCTGGCCTCCCACCCCCGTTTTG‑CCTGCAAGCGCGCAAA‑TGTGG‑CCTACCCCC‑AGCACCGCCGTACCGtt                                                                                                                              <  1:903346/144‑1 (MQ=255)                             gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                 >  1:517613/1‑151 (MQ=255)                             gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                 >  1:679727/1‑151 (MQ=255)                                tgggTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGccacc                                                                                                              >  1:324352/2‑151 (MQ=255)                                 gtgTCGAGTTCGTCAACTTCAATGTAAAATAAATCACCCCCGGCAATTTTTTA‑AAACAGACCCCCTTTGCCTCCCAACCCCGGTTGC‑CCGGCCAGCGCGCCAA‑TTTTC‑CCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGccacca                                                                                                             <  2:258244/149‑1 (MQ=255)                                 ggttCGATTCCGCCAACTCCAATTTAAAATACACAAACCCCGGCGATTTTTTA‑AAACAAACCCCCCTCGCCTCCCACCCCGGTTTTG‑CCTGCCGGCGCCCCAA‑TGTTC‑CCCCCCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGccacca                                                                                                             <  1:531453/148‑1 (MQ=255)                                         ttCGTCAACTTCAATGTAAAATACATCAACCCCGGCAATTTTCTA‑AAACAGAACCCCTTTGCCTTCCATCCCCGGTTTG‑CCTGCCAGCGCGCAAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                        <  2:1059796/148‑1 (MQ=255)                                               aaCTTCATGGCCAAAATCATCATTCCCGCAGTTTTTCAT‑AACCAGACCTCCTTGGCCTCCAAACCCGGATTTG‑CCTACCGCCTCGCGAA‑TGTCC‑CCTACCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                               <  2:693371/151‑1 (MQ=255)                                                     attGTCAAATACATCACCCCCGGCGATTTTCTA‑AAACAGACCGCCCTTCCCTCCCATCCCCGGTTTC‑CCTCCCAGCGCGCCAA‑TGTTG‑CCTACCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtata                                                                                          <  1:360183/148‑1 (MQ=255)                                                     aaTGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                         >  1:377800/1‑151 (MQ=255)                                                           aaaTACATCAATACCGGCGAATTTCTT‑AAACAGCACGCCCTTACCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGCTACCTCGGGCGGTATATTTGt                                                                                     >  1:947745/1‑149 (MQ=255)                                                            aatccacaaaccccGGCAACTTCCTA‑AAACAAACCGCCCTCCCCTCCCCCCCCCGTTTTC‑CCGCCCTGCCCCCCAT‑TTTTG‑CCTACCCCC‑ACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                                   <  1:907553/139‑1 (MQ=255)                                                             atcaggcAACCCCGACGAATGTTTA‑AAACAAAGCCCCCTTGCTTTCCACCGCCGTTTAG‑CTTCCCGCCGCCCAAT‑TGTTC‑CCTCCCCCC‑ATCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                    <  2:325412/142‑1 (MQ=255)                                                                 accAACCCCGGCAAATTTCTA‑AAACAGACCGCCTTGCCCTCCCCTCCCCGTTTGG‑CCGGCCGGCCCGCCAA‑GGTGG‑CTTACCCCC‑GGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                              <  1:709444/148‑1 (MQ=255)                                                                 aacAACCCCGGCGAATTTCTA‑AAACAGACCCCCCTTCCCTTCCACCCCCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                              <  2:436364/148‑1 (MQ=255)                                                                      ccccGGCGAA‑TTCTT‑AAAAAAAACCCCCTTTCCTTCCCACCCCCGTTTGCCCTGCCGGCGCCCCAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTa                                                                            <  1:1039527/146‑1 (MQ=255)                                                                       cccGGCGAATTTCTA‑AAACAGACCGCCCTTGCCTCCCATCCCCGTTTTG‑CCTGCCAGCGCGCCAT‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                       <  2:560737/151‑1 (MQ=255)                                                                       cccGGCAAATTTCTT‑AAACAGACCGCCCTTGCCTCCCATCACCGGTTTG‑CCTGCCAGCGCGCAAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                       <  2:180817/151‑1 (MQ=255)                                                                        ccgggaaaTTTTTTAAAACAGACCGCCCTTGCCTCCCACCCCCGTTTTG‑CCTGCCAGCGCCCAAA‑TGTTC‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                         <  2:301811/143‑1 (MQ=255)                                                                        ccGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                       >  1:476293/1‑150 (MQ=255)                                                                          ggCGATTTTCTA‑AACCAGGTCCCCTTTGCCTCCCATCCCATGTTTC‑CCGGCCGCCGCGCAAT‑TGTTG‑CCCAACCCC‑AGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                                     <  2:11287/150‑1 (MQ=255)                                                                          ggCGAATTTCTA‑AAACAGACCCCCCTTGCCTCCCACCACCGTTTTG‑CCGGCCAGCGCGCCAA‑TGTTC‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                                    <  2:546188/151‑1 (MQ=255)                                                                              aaTTTTTA‑AAACAGACCGCCCTTCCCTCCCATCCCCGGTTTG‑CCTGCCGGCGCGCCAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                   <  2:997168/148‑1 (MQ=255)                                                                                     taaaaCAGCACCCCCTT‑CCCTCCCTCCCCGGTTTT‑CCCTCCGGGCCCCCAA‑ATTTGCCCTACCCCC‑GCCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                           <  1:613100/148‑1 (MQ=255)                                                                                       aaacaaacCCCCCTTCCCTCCCACCCCGGGTTTG‑CCGCCCAGCGCCCAAA‑TGTTG‑CCTACCCCC‑GGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTCTTAACGGGTCTTTTTCGATTTCa                                                        <  1:907238/151‑1 (MQ=255)                                                                                       aaaCAGACCGCCCTTGCCTTCCACCCCCGGTTTG‑CCTCCCAGCGCGCCAA‑TGTTG‑CCTACCCCC‑ACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                         <  2:542469/150‑1 (MQ=255)                                                                                        aaCAGACCGCCTTTCCCTCCCACCCCCGGTTGG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCaa                                                       <  2:734197/151‑1 (MQ=255)                                                                                        aaCAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                        >  1:10776/1‑150 (MQ=255)                                                                                         aCAGACCGCCCTGCCCTCCCACCCCCGTTTTC‑CCTCCCCGCCCCCCAT‑TTTTC‑CCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAg                                                      <  1:4843/151‑1 (MQ=255)                                                                                         aCAGACCCCCCTTCCTTCCCACCCCCGTTTTC‑CCTCCCGGCGCCGCAT‑TGTTC‑CCTACCCCC‑GGAGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                              <  2:696930/143‑1 (MQ=255)                                                                                           agacccctcttcccttccACCCCCGTTTTG‑CCGGCC‑GCCCCCCAATTTTTG‑CGTACCCCA‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                    <  1:582481/143‑1 (MQ=255)                                                                                            gAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                    >  2:792762/1‑150 (MQ=255)                                                                                             accccccTTGCCTTCCATCCCCGTTTTG‑CCTGCCGGCGCGCAAA‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                   <  2:333778/146‑1 (MQ=255)                                                                                             accGCCCTTCCCTCCCTCCCCGGGTTGG‑CCTGCCAGCGCCCCAT‑TGTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGa                                                     <  2:390407/146‑1 (MQ=255)                                                                                               cccccTTCCCTCCCACCACCGTTTTG‑CCTGCCAGCGCGCCAT‑TGTTG‑CCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                <  2:688351/149‑1 (MQ=255)                                                                                                 cccTTCCCTCCCATCCCCGTTTTC‑CCTGCCAGCGCCCAAA‑TTTTC‑CCAACCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                <  2:290402/149‑1 (MQ=255)                                                                                                    ttGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                           >  2:799006/1‑151 (MQ=255)                                                                                                     tGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc                                          >  2:324483/1‑151 (MQ=255)                                                                                                       ccTCCCACCCCCGTTTTC‑CCTCCAGGCGCGCCAA‑TGTGG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                         <  1:140673/150‑1 (MQ=255)                                                                                                       ccTCCCACCCCCGGTTTG‑CCTGCCAGCGCGCAAA‑TGTTC‑CCAAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCaa                                        <  1:153628/151‑1 (MQ=255)                                                                                                               ccccGTTTTG‑CCTCCCGCCGC‑CCAA‑TTGTG‑CCTAACCCCGCGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGt                                <  1:656740/149‑1 (MQ=255)                                                                                                                          cctgccAGCGCCCAAA‑TGTTG‑CCTACCCCA‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAgg                         <  1:383912/148‑1 (MQ=255)                                                                                                                          cctgccAGCGAGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                       >  1:264153/1‑150 (MQ=255)                                                                                                                            tgccAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgc                    >  1:1028340/1‑151 (MQ=255)                                                                                                                              ccAGCGCCCAAA‑TTTTG‑CCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc                  <  2:369865/151‑1 (MQ=255)                                                                                                                                  cgcgcCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                 >  1:25075/1‑148 (MQ=255)                                                                                                                                     gcCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                 >  2:835621/1‑145 (MQ=255)                                                                                                                                         a‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt       >  1:831287/1‑151 (MQ=255)                                                                                                                                            gTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtgt     >  2:734876/1‑151 (MQ=255)                                                                                                                                               g‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGtt  >  2:642901/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTTCCATCACCGGTTTG‑CCTGCCAGCGCGCCAA‑TGTTG‑CCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTT  >  NC_000913/2577999‑2578288

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