Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	47.5%	77.1 / 27.8	61	L79V (TTA→GTA)	maeB	malic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (25/7);  new base C (0/29);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 3.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577994‑2578291                                                                                                                                                          |                                                                                                                                                        tggttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTaa                                                                                                                                                         <  2:90159/150‑1 (MQ=255) tggttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTaa                                                                                                                                                         <  1:338709/150‑1 (MQ=255)            cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CAccgc                                                                                                                                          <  2:419897/151‑1 (MQ=255)            cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CAccgc                                                                                                                                          >  2:564787/1‑151 (MQ=255)              gcgACAACTTAAATAAATTTGTCCGGGTCGAGTGCGTCAACTTCAATGTCAAATACAGCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCACCC‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CAccgccg                                                                                                                                        <  1:776783/151‑1 (MQ=255)                gACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGt                                                                                                                                       >  2:907251/1‑150 (MQ=255)                  cAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTAcc                                                                                                                                    >  1:76638/1‑151 (MQ=255)                   aaCTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGACTTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCg                                                                                                                                   >  1:686387/1‑151 (MQ=255)                         aaTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTag                                                                                                                               >  2:522710/1‑149 (MQ=255)                         aaTAAATTTGCCCGGGTCGATTTCTCCAACTTCATTGTAAAATACATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTGCCTCCCATCC‑CCGTTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAC‑CCCCA‑C‑CACCGCCGTACCGtt                                                                                                                                 <  1:630485/147‑1 (MQ=255)                               tttGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTCCCATCC‑CCGTTTTG‑CCTG‑CCA‑GCGCGCCATTGTTGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCAc                                                                                                                        <  2:22181/150‑1 (MQ=255)                               tttGCCGGGGGCGATTTCGTCATCTTCAATGTAAAATACACCAATCCCGGCAAATTTCTAAAACAGACCCCCCTTGCCTCCCATCC‑CCGGTTTG‑CCTG‑CCA‑GCGCGCAAATGTTGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCAc                                                                                                                        <  2:1490/150‑1 (MQ=255)                                   tCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                    >  1:912151/1‑150 (MQ=255)                                      gggTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGccacc                                                                                                                 >  1:332277/1‑150 (MQ=255)                                       ggTCGAGTTCGTCACCTACAATGTAAAATAAATCAACCCCGGCAATTTTCTAAAAAAGACCCCCCTTCCCTCCCA‑CC‑CCGGTTTGCCCTC‑CCG‑GCGCGCCAATGTGCCCAAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAg                                                                                                               <  2:183460/151‑1 (MQ=255)                                        gTCGATTTCGTCAACTTCAATGTCAAATACATCAACCCCGGCGATTTTCTTAAACAGACCGCCCTTCCCTCCCATCC‑CCGGTTTG‑CCTC‑CCG‑GCGCGCCAATGTTGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAgg                                                                                                              <  1:385109/151‑1 (MQ=255)                                            aGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTAc                                                                                                          >  2:250920/1‑151 (MQ=255)                                              ttcctaaaCTTAAATGTAAAATACATCAACCCCGGCGATTTTCTAAAACAGACCGCCCTTGCCTCCCCCCA‑CCGGTTTG‑CCGG‑CCG‑GCGCGCCAATTTTGCCAAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                            <  2:565039/141‑1 (MQ=255)                                              ttCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                                        >  1:726440/1‑151 (MQ=255)                                                cgccaccTTCAATGTAAAATACATCAACCCCGGCGATTTTATAAAACAAACCCCCTTTCCCTTCCCTCC‑CCGTTTTC‑CCTC‑CCA‑GCGCACAAATTTCCCCACC‑CCCCA‑G‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                                        <  1:892498/143‑1 (MQ=255)                                                    aactaaaatgaaaaaTACACCACCCCCGGCAATTTTCTAAAACAGACCCCCCTTGCCTCCCATCC‑CCGTTTTG‑CCGC‑CCG‑GCGCCCAAATTTTCCCTAC‑CCCCA‑C‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                                       <  1:334535/134‑1 (MQ=255)                                                    aaCTTCATTTTCAATTCCACCAACCCCGGCGATTTTCTTAACCAGACCCCCCTTCCTTCCCATCC‑CGGTTTTG‑CCTG‑CCG‑CCGCCCAATTGTTGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                   <  2:313385/150‑1 (MQ=255)                                                      cTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                                >  2:163480/1‑151 (MQ=255)                                                       ttaaaTGTAAAATACATAACCCCCGGCAAATTTCTAAAACAGAACCCCCTTGCCTCCCACCC‑CCGTTTTG‑CCTG‑CCA‑GCGCGCCAATTTTGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                   <  2:436310/144‑1 (MQ=255)                                                          aaTGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                            >  2:280231/1‑151 (MQ=255)                                                          aaTGTAAAATACATAAACCCGGGCAATTTTCTAAAACAGACCCCCCTTCCCTCCCACCC‑CCGTTTTG‑CCTC‑CCG‑GCGCGCAATTTTTCCCTAC‑CCCAA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                            <  2:839610/151‑1 (MQ=255)                                                                aaaTACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                        >  2:228968/1‑149 (MQ=255)                                                                 aaTACATCATCCCCGGCGAATTTCTAAAAAAGAAGGCCTTGGCCTTCCATCC‑CCGCTTGG‑CCCG‑CCC‑GCGCGCGATTGTTCCCAAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                     <  2:956521/151‑1 (MQ=255)                                                                         acccccGGCAAATTTCTAAAACAGAACCCCCTTCCCTTCCATCC‑CCGTTTTC‑CCTC‑CCA‑GCCCGCAATTGTTGCCTAC‑CCCCA‑G‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                                <  2:64463/145‑1 (MQ=255)                                                                         aaccccGGCGAATTTCTAAAACAGACCCCCCTTCCCTCCCACCC‑CCGTTTTG‑CCTG‑CCA‑GCGCGCAATTGTTGCCAAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTaa                                                                              <  2:463124/147‑1 (MQ=255)                                                                           ccccgcaaaaTTTCTAAAAAAAACCCCCCTTCCCTCCCACCC‑CCGTTTTG‑CCTG‑CCG‑GCGCGCAATTGTTCCATAC‑CCCAA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                             <  2:699228/141‑1 (MQ=255)                                                                           acccGGGCAATTTCCTAAACAACACCCCCCTGCCCCCCCCCCCCCGGTTTC‑CCGC‑CCA‑CCGCGCAATTGTGCCCTAC‑CCCCA‑G‑CGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt                                                                                   <  1:716290/143‑1 (MQ=255)                                                                            cccGGCGATTTTCTAAAACAGAACGCCCTTCCCCCCCATCC‑CCGTTTTC‑CCTG‑CCA‑GCGCGCCAATTTGCCCTCC‑CCCCA‑G‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTaa                                                                              <  1:126481/147‑1 (MQ=255)                                                                            cccGGCAATTTTCTAAAACAG‑ACCCCCCTGCCCTCCAACC‑CCGGTTTT‑CCCGCCCA‑GCGCCCCAATGTCGCCTAC‑CCCCG‑G‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                                  <  1:139513/143‑1 (MQ=255)                                                                            cccGCCAAATTTCTCAAACACACCCCCCTTGCCTTCCATCC‑CCGGTTAG‑CCGC‑CCC‑GC‑CGCCAAATTTCCCCTA‑CCCCACG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                                  <  2:468869/143‑1 (MQ=255)                                                                               ggCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                                       >  2:665616/1‑151 (MQ=255)                                                                                 cGATTTTCTAAAACAGACCCCCTTTCCCTCCCACCC‑CCGTTTTG‑CCTC‑CCA‑GCGCGCCAATGTGCCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                      <  2:37031/150‑1 (MQ=255)                                                                                   aaTTTCTAAAACAGACCGCCTTTGCCTCCCACCC‑CCGTTTTG‑CCTG‑CCG‑GCGCGCCAATGTTGCCTAC‑CCCCA‑C‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                    <  2:751125/150‑1 (MQ=255)                                                                                     tttttaaaacaaacCCCCCTTCCCTCCCACCC‑CCGGTTTC‑CCGG‑CCG‑GCGCCCCAATTTTCCCAAC‑CCCCC‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                  <  1:428639/144‑1 (MQ=255)                                                                                     tttCTAAAAAAGACCCCCCTT‑CCCTCCCCCC‑CCGGTTTGCCCTC‑CCA‑GCCCCCAATTGTTCCCAAC‑CCCCG‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                    <  1:845649/148‑1 (MQ=255)                                                                                      ttCTAAAACAGACCGCCCTT‑CCCTCCCCCC‑CCCGTTTG‑CCCGCCCA‑GCGCCCAATTGTTGCCAAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                                 <  1:665912/150‑1 (MQ=255)                                                                                           aaaCAGACCGCCCTTCCCTTCCATCT‑CCGGTTTC‑CCTG‑CCG‑CCGCTCAAGTGTTCCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                               <  2:38028/147‑1 (MQ=255)                                                                                              cAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           >  1:752919/1‑148 (MQ=255)                                                                                                   cccccTTGCTTCCCATCC‑CCGTTTTC‑CCCC‑CCA‑GCCCGCAAATTTTCCTTAC‑CCCCA‑C‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAg                                                    <  2:481/148‑1 (MQ=255)                                                                                                      ccTTCCCTCCCATCC‑CCGGTTTG‑CCTG‑CCG‑GCGCGCAAATGTTGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                                <  2:76636/151‑1 (MQ=255)                                                                                                        ttGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                              >  2:702044/1‑151 (MQ=255)                                                                                                        ttGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGGGCg                                               >  1:728679/1‑150 (MQ=255)                                                                                                           ccTCCCATCC‑CCGTTTTC‑CCTG‑CCA‑GCGCGCCAATTTCGCCTAC‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                            <  2:332284/150‑1 (MQ=255)                                                                                                            cTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                          >  1:178123/1‑151 (MQ=255)                                                                                                             ttCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                         >  2:461165/1‑151 (MQ=255)                                                                                                             ttCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                         >  1:57819/1‑151 (MQ=255)                                                                                                                ccctcccggttttC‑CCTC‑CCA‑GCGC‑CCAATTTTTCCCAACCCCCG‑G‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGcc                                       <  1:564776/141‑1 (MQ=255)                                                                                                                cATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                                      >  1:218120/1‑151 (MQ=255)                                                                                                                     ccttttttG‑CCCT‑CCA‑GCGCTCCGTTTTTTCAAAA‑CCCCG‑C‑CACCGCCGTGCCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                  <  2:660798/146‑1 (MQ=255)                                                                                                                       cGGTTTT‑CCCG‑GCAGGCGCGCAATTGTT‑CCCAA‑CCCCA‑GCCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTg                                     <  1:535038/147‑1 (MQ=255)                                                                                                                          tttG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                               >  1:15900/1‑149 (MQ=255)                                                                                                                          ttgg‑ccgg‑ccg‑GCGCCCCAATTTTGCTTAA‑CCCCG‑C‑CTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                               <  2:784744/146‑1 (MQ=255)                                                                                                                               tgccgccccgcccgcCATTGTTCCCTAC‑CCCCA‑C‑CTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAgg                            <  2:686403/138‑1 (MQ=255)                                                                                                                                ctg‑ccA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                          >  2:96167/1‑149 (MQ=255)                                                                                                                                        gcgcgcCATTGTTCCCTAC‑CCCAA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAg                  <  1:125924/151‑1 (MQ=255)                                                                                                                                               aaTGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc           >  2:323787/1‑151 (MQ=255)                                                                                                                                               aaTGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc           >  2:856414/1‑151 (MQ=255)                                                                                                                                                 tGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTCACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtg         >  1:219215/1‑151 (MQ=255)                                                                                                                                                         aa‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc  >  1:594472/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                        TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCA‑CCGGTTTG‑CCTG‑CCA‑GCGCGCCAATGTTGCCTAA‑CCCCA‑G‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577994‑2578291

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