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

GCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTG  >  NC_000913/2578004‑2578289                                                                                                                                                  |                                                                                                                                                    gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCAccg                                                                                                                                         <  1:709562/151‑1 (MQ=255)     cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGt                                                                                                                                     >  1:765213/1‑151 (MQ=255)      gACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCTTCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCAccgccg                                                                                                                                      >  1:447465/1‑149 (MQ=255)      gACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGt                                                                                                                                     >  2:142444/1‑150 (MQ=255)       aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTAc                                                                                                                                   >  2:645267/1‑151 (MQ=255)         aaCTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCg                                                                                                                                 >  1:264745/1‑151 (MQ=255)            ttaaataatTTTGTCGGGGTGGATTTCGTCAACTTCAATGTCAAATACACCAACCCCGGGGAA‑TTTCTAAAA‑CAGAACC‑CCC‑TTGCCTCCCACCCCCGTTTTG‑CCGG‑CCAGCGCGCC‑ATGTTTCCCAACCCCCAGCCCCGCCGTACCGt                                                                                                                                <  1:438118/146‑1 (MQ=255)            ttCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGtt                                                                                                                               >  1:304977/1‑150 (MQ=255)             tCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGtt                                                                                                                               >  2:64023/1‑149 (MQ=255)                aaaaaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCGAA‑TTTCTAAAA‑CAGAACG‑CCC‑TTGCCTTCCATCCCCGGTTTG‑CCTG‑CCAGCGCGCCAATGTGGCCTAC‑CCCCAGCACCGCCGTACCGTTAgaga                                                                                                                          <  2:50468/149‑1 (MQ=255)                          ccGGTTCGATTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCGAA‑TTTCTTAAA‑AAAACCC‑CCC‑TTGCCTCCCACCCCCGTTTTG‑CCTC‑CCAGCGCGCAAATTTTCCCAAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCg                                                                                                                    <  2:925839/147‑1 (MQ=255)                          ccGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGccac                                                                                                                >  2:685850/1‑151 (MQ=255)                          ccGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGccac                                                                                                                >  1:248683/1‑151 (MQ=255)                            gggTCGTTTTCTTCAACTTAAATTTAAAATACACCACCCCCGG‑GAA‑TTTTTTAAACCAGAACC‑CCC‑TTCCCTTCCACCCCCGTTTCC‑CCGC‑CCGGCCCGCCATTGTCCCCAAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGccacc                                                                                                               <  1:172553/150‑1 (MQ=255)                                cGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                          >  2:607246/1‑151 (MQ=255)                                       gTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                                   >  2:771686/1‑151 (MQ=255)                                             atCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGATTG‑CCTG‑CCAGCGCGCCAATGTTTCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                              >  1:536672/2‑150 (MQ=255)                                              tCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                            >  1:83283/1‑151 (MQ=255)                                                aatgaaaaaTAAACCAACCCCGGC‑AA‑TTTTCTAAA‑AAAAACC‑CCC‑TTTCCCTCCCACCCCGGTTTGCCCTC‑CCGGCGCCCCATTGTTCCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGta                                                                                             <  1:787376/142‑1 (MQ=255)                                                aaTGTAAAAAACA‑CAACCCCCGGGAATTTTCTAAAA‑CAGACCC‑CCC‑TTGCCTCCCATCCCGGTTTTG‑CCTG‑CCAGCGCGCCAATTTTCCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGta                                                                                             <  2:1028042/148‑1 (MQ=255)                                                      aaaTACATCAACCCCGGCGAA‑TTTCTAAAA‑CAGACCC‑CCC‑TTCCCTCCCACCCCGGTTTTC‑CCTC‑CCAGCGCCCCAATGTCGCTTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAtt                                                                                         <  2:571422/146‑1 (MQ=255)                                                      aaaTACA‑CAAACCCGGGAAATTTTCTAAAA‑CAGACCC‑CCT‑TTCCCTCCCACCCCCGTTTTC‑CCTG‑CCAGCGCCCAATTGTTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                      <  2:618736/149‑1 (MQ=255)                                                       aaTACTCCAACCCCGGTTAA‑TTTTTAAAA‑AAGACCC‑CCT‑TTGCCTCCCACCCCGGTTTTG‑CCGG‑CCGGCCCCCCAATTTTGCTTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                   <  2:945583/151‑1 (MQ=255)                                                       aaTACATCAACCCCGGCGAT‑TTTCTAAAA‑CAGACCC‑CCC‑TTGCCTCCCTTCCCCGTTTTG‑CCGG‑CCAGCGCGCAATTGTTGCCAAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                   <  2:115572/151‑1 (MQ=255)                                                       aaTACACCAACCCCGGCTAA‑TTCCTAAAA‑CAGCACC‑CCC‑TTCCCTCCCCCCCCCGTTTTG‑CCTG‑CCGGCCCACCTTTGTGTCCTAC‑CCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                   <  1:119366/151‑1 (MQ=255)                                                            aTCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                              >  1:630030/1‑151 (MQ=255)                                                             accaccccgggcaattttCTAAAA‑CAAAACC‑CCT‑TTGCCTTCCACCCCGGTTTTG‑CCTG‑CCAGCGCCCAATTTTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                               <  2:1030979/135‑1 (MQ=255)                                                                 tCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                         >  1:780506/1‑151 (MQ=255)                                                                   ccGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                        >  1:722222/1‑150 (MQ=255)                                                                     ggCGAA‑TTTCTTAAA‑CAGACCG‑CCC‑TTGCCTCCCATCCCCGTTTTG‑CCTG‑CCAGCGCCCCATTTTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                                      <  2:644581/150‑1 (MQ=255)                                                                      gCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                    >  1:999926/1‑151 (MQ=255)                                                                       cGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt                                                                   >  1:741147/1‑151 (MQ=255)                                                                             tttttAAA‑AAAAACC‑CCCTTTGCCTTCCATC‑CCGGTTTGCCCTG‑CAGGCGCCCCAGTTTCCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                <  2:78740/147‑1 (MQ=255)                                                                             tttctAAA‑CAAAACCCCCC‑TTCCCTCCCACCCCCGTTTTC‑CCTC‑CCGGCGCGCAATTTTTCCCTAC‑CCCCGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                <  1:938146/146‑1 (MQ=255)                                                                             ttCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                               >  2:1001125/1‑150 (MQ=255)                                                                                 tAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                          >  2:154949/1‑151 (MQ=255)                                                                                  aaa‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                          >  2:1033303/1‑150 (MQ=255)                                                                                  aaa‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                          >  2:812756/1‑150 (MQ=255)                                                                                      acaaaCC‑CCC‑CTTCCCTCCCTCCCCCGTTTGCCCGG‑CCAGCGCCCCAATTTTCCCTAC‑CCCCGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAg                                                       <  1:156532/148‑1 (MQ=255)                                                                                        aaacccccc‑TTCCCTCCCCCCCCCGTTTTG‑CCTG‑CCGCCGCCCAATTGTTGCCAAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                     <  1:6048/145‑1 (MQ=255)                                                                                          acccccc‑ttcccttCCATCCCCGGTTTG‑CCTC‑CCAGCGCCCCATTGTTGCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                    <  2:76438/146‑1 (MQ=255)                                                                                           accccc‑TTTCCCTCCCCCCCCGGGTTG‑CCCGCCCGGCCCCCCATTTTTCCCTAC‑CCCAGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACaa                                                   <  1:639527/148‑1 (MQ=255)                                                                                            ccccc‑TTGCCTCCCCCCCCGGTTTTC‑CCCG‑CCAGCGCCCAAGTGTTCCCAAC‑CCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAg                                                  <  2:588985/148‑1 (MQ=255)                                                                                            ccccc‑TGGCCTCCCATCCCCGGTTTG‑CCTG‑CCGGCGCCCAAATGTTCCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAg                                                  <  2:693924/148‑1 (MQ=255)                                                                                                  ttGCTTTCCATCACGGGTTGG‑CCGC‑CCAGGGCCCAAATGTGCCCAAC‑CCCCAGAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAg                                                  <  2:881285/145‑1 (MQ=255)                                                                                                     ccTCCCACCCCCGTTTTG‑CCTG‑CCAGCGCGCAATTGTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                            <  2:349866/148‑1 (MQ=255)                                                                                                      cTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                        >  2:161871/1‑151 (MQ=255)                                                                                                          cacccccGGTTTG‑CCTG‑CCAGCGCGCCATTGTTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                                    <  2:470101/146‑1 (MQ=255)                                                                                                              aCCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                >  1:133831/1‑151 (MQ=255)                                                                                                                   tttG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                             >  2:493831/1‑149 (MQ=255)                                                                                                                   tttG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                             >  2:175572/1‑149 (MQ=255)                                                                                                                        cctg‑ccAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                        >  1:887038/1‑150 (MQ=255)                                                                                                                        cctg‑ccAGCGCGCAAATGTTCCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAgg                          <  2:868430/148‑1 (MQ=255)                                                                                                                               aGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                  >  1:586216/1‑150 (MQ=255)                                                                                                                                ggacccaatttttcccTAC‑CCCCAGAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCa                 <  2:940328/136‑1 (MQ=255)                                                                                                                                 cgcgcAAATGTTCCCCAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAg                      <  2:572347/144‑1 (MQ=255)                                                                                                                                    cccAATGTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcg            <  2:305098/150‑1 (MQ=255)                                                                                                                                     caattGTTCCCTAC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATc             <  2:777968/145‑1 (MQ=255)                                                                                                                                      cAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg          >  2:76514/1‑151 (MQ=255)                                                                                                                                       aaTGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc         >  2:298507/1‑151 (MQ=255)                                                                                                                                        aTGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGTTTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg          >  2:546488/1‑149 (MQ=255)                                                                                                                                         tGTTCCCAAC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATc             <  1:390697/145‑1 (MQ=255)                                                                                                                                           ttcccaaa‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtgt      <  1:525645/144‑1 (MQ=255)                                                                                                                                              ccTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTg  >  1:732289/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAACG‑CCC‑TTGCCTTCCATCACCGGTTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTG  >  NC_000913/2578004‑2578289

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