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

AGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2578001‑2578291                                                                                                                                                  |                                                                                                                                                      aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑c                                                                                                                                               >  1:783672/1‑150 (MQ=255)  gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CAc                                                                                                                                             <  1:617560/151‑1 (MQ=255)   cgcgGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CAcc                                                                                                                                            >  2:998772/1‑151 (MQ=255)   cgcgGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CAcc                                                                                                                                            >  2:859008/1‑151 (MQ=255)      ggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CAccgcc                                                                                                                                         <  2:427457/151‑1 (MQ=255)       gcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CAccgccg                                                                                                                                        <  1:56847/151‑1 (MQ=255)            aaCTTAAAAAAATTTGTCGGGGTCGAGTTCGTCAACTTCAATGT‑AAAATACACCAACCCCGGCGATTTTCTTAAACAGACCCCCCTTCCCTTCCATCCCCGTTTT‑GCCTGCCGGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCg                                                                                                                                   <  1:168616/151‑1 (MQ=255)                  aataatTTTGTCGGGGTCGTGTTCGCCAACTTCAATGT‑AAAATACACCAACCCCGGCGAATTTCTTAAACAGACCCCCCTTGCCTCCCATCCCCGTTTT‑GCCTGCCAGCGCGCAATTG‑TTGCCTACCCCCAG‑CACCGCCGTACCGtt                                                                                                                                 <  2:187580/147‑1 (MQ=255)                   aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTag                                                                                                                               >  1:817742/1‑148 (MQ=255)                     aattttGTCCGGGTCGAGTTCGTCAACTTAAATGT‑AAAATACACCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTCCCTCCCATCCCCGGTTT‑GCCTGCCAGCGCGCAAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATc                                                                                                                          <  1:620772/148‑1 (MQ=255)                      aaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATc                                                                                                                          >  1:820163/1‑150 (MQ=255)                      aaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATc                                                                                                                          >  1:555386/1‑150 (MQ=255)                         ttGTCCGGGTCGGGTTTCTCAACTTCAATTTAAAAATACATCAACCCCGGCAAATTTTTAAAACAGACCGCCCTTCCCTTCCACCCCCGTTTT‑GCCTGCCAGCGCGCAATTT‑TTGCCTACCCCCAG‑CACCGCCGTACCGTTAGAGATCAc                                                                                                                        <  1:452548/150‑1 (MQ=255)                         ttGTCCGGGTCGATTTCGTCAACTCCAATTA‑AAAATACATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTCCCTTCCACCCCCGTTTT‑GCCTGCCAGCGCCCCAATG‑TTCCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCAcc                                                                                                                       <  2:507191/150‑1 (MQ=255)                           gTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                    >  1:447690/1‑151 (MQ=255)                                ggTCGTGTTCTTCAACTTAAATGT‑AAAATACATAAACCCCGGAAATTTTTTAAAACAGACCCCCCTTCCCTCCCATCCCCGTTTT‑GCCTGCCAGCCCCCCAATT‑TTCCCAACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAg                                                                                                               <  1:178754/151‑1 (MQ=255)                                   cgattccGTCAATTTAAATTT‑AAAAAACATAAACCCGGGCGAATTTCTAAAACAGACCGCCCTTTCTTTCCACCCCGGTTTT‑GCGTCCCACCGCGCCATTT‑TTCCCTACCCCAAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                             <  1:898968/144‑1 (MQ=255)                                         cgccAACTAAAATGT‑AAAATACACCACCCCCGGCGATTTCCAAAAACAAACCGCCTTTCCCTCCCA‑CCCCGGTTTGGCCTGCCGGCGCCCCATTG‑TCCCCAACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                                        <  2:1029633/146‑1 (MQ=255)                                          gTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                                     >  2:632547/1‑151 (MQ=255)                                            cAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                   >  2:1018800/1‑151 (MQ=255)                                             aaCTTCAATGT‑AAAATACACCAACCCCGGCAAATTTCTAAAACAGACCGCCCTTCCCTCCCACCCCCGTTTT‑GCCTCCCGGCGCGCCATTG‑TTCCCTACCCCCAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                   <  1:71519/150‑1 (MQ=255)                                              acctacaataaaaaaCACATCAACCCCGGCAATTTTCATAAACACACCGCCCTTACCTCCCACCCCCGTTTT‑GCCGGCAACCGCGCC‑ATG‑TCCCCGAACCCCCGCCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                                        <  2:201424/133‑1 (MQ=255)                                                ttCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                                >  1:843871/1‑150 (MQ=255)                                                ttCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                                >  2:884906/1‑150 (MQ=255)                                                ttCAATGT‑AAAATACATCACCCCCGGCGAATTTCTAAAAAAGACCGCCCTTGCCTCCCACCCCCGGTTT‑GCCTGCCAGCGCGCCATTG‑TTCCCTACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                                       <  1:585823/143‑1 (MQ=255)                                                 tCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                              >  1:1032012/1‑151 (MQ=255)                                                    aatgaaaaaTACACCAACCCCGGAAATTTTCAAAAAAAGACCGCCTTTGCTTCCCATCCCCGTTTG‑CCCGGCCGGCGCGCAATTG‑TGCCCTACCCCCAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGta                                                                                               <  1:652552/142‑1 (MQ=255)                                                     tGT‑AAAATACACCACCCCCGGCAATTTTCTAAAACAAACCCCCCTTCCCTCCCCCCCCCGTTTT‑GCCTCCCGGCGCCCCATTT‑TGCCCTCCCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                                  <  2:305777/143‑1 (MQ=255)                                                       t‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                        >  1:313228/1‑151 (MQ=255)                                                          aaaaacaaaatccccGGGGATTTTCTAAAACAAACCCCCTTTGCCTTCCACCCCGGTTTT‑CCCTCCCCGCGCGCAATTG‑TTCCCTACCCCAGG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                        <  1:861334/137‑1 (MQ=255)                                                           aaTACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                     >  2:525512/1‑151 (MQ=255)                                                             tACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt                                                                                   >  1:735092/1‑151 (MQ=255)                                                              aCATCAACCCCGGCAAATTTCTAAAACAGACCCCCCTTCCCTCCCATCCCCGGTTT‑GCCTGCCAGCGCCCAAATG‑TTCCCTACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGc                                                                                    <  1:966116/149‑1 (MQ=255)                                                                   acccccgggaattttCTTAAACAAACCCCCCTTCCCTCCCTCCCCCGTTTT‑GCCTGCCCGCCCCCAATTG‑TGGCCTCCCCCACG‑CTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                                <  2:555199/136‑1 (MQ=255)                                                                    accccGGCAATTTTCTAAAACAGACCGCCCTTCCCTCCCACCCCCGTTTT‑CCCGGCCGCCGCG‑CAATGTTTCCCTACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTaa                                                                              <  2:620419/147‑1 (MQ=255)                                                                      cccGGCGAATTTCTAAAACAGACCGCCCTTGCCTCCCACCCCCGTTTT‑GCCGGCCGGCGCGCAAATT‑TGGCCAACCCCCAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                           <  1:523314/150‑1 (MQ=255)                                                                           cGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt                                                                     >  2:819309/1‑151 (MQ=255)                                                                             attttCTTAAAAAAACCCCCCTTCCTTCCCACCCCCGTTTT‑CCCCGCCGGGCCGAAATTG‑TTCCCAACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt                                                                     <  1:323397/147‑1 (MQ=255)                                                                             attttCTTAAAAAAAACCCCGTTGCCTCCCACCACGGTTTC‑GCCTGCCACCGCCCCAATT‑TCCCCAACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt                                                                     <  2:769123/147‑1 (MQ=255)                                                                              aTTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCAGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                  >  1:1005294/1‑151 (MQ=255)                                                                                  cTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                              >  2:544427/1‑151 (MQ=255)                                                                                   taaaaCAGACCCCCTTTGCTTCCCACCCCCGTTTT‑GCCTCCCAGCGCCCAATTG‑TTCCCAACCCCCAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                               <  1:1005208/147‑1 (MQ=255)                                                                                   taaaaCACACCCCCCTTCCCTCCCCTCCCCGTTTG‑TCCTCCCACCGCCCAATTG‑TTTCATACCCCCGC‑AACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                                 <  1:820907/145‑1 (MQ=255)                                                                                      aaCAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           >  2:14861/1‑150 (MQ=255)                                                                                           accccccTTCCCTCCCACCCCCGTTTT‑GCCTGCCAGCGCCCAAATG‑TTCCCAACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACaa                                                     <  1:840613/147‑1 (MQ=255)                                                                                           accGCCCTTCCTTCCCATCCCCGGTTG‑CCCTGCCAGCGCCCAAATG‑TTGCCAACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                      <  1:618877/148‑1 (MQ=255)                                                                                               cccTTCCCTCCCATCCCGGTTTT‑GCCTGCCAGCGCGCAAATT‑TGCCCTACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTg                                                 <  1:341908/151‑1 (MQ=255)                                                                                                  ttcccTCCCACCCCCGGTTT‑CCCTCCCAGCGCCCCAATG‑TTGCCTACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                                <  1:998785/146‑1 (MQ=255)                                                                                                   tGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc                                             >  2:82909/1‑151 (MQ=255)                                                                                                    gCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc                                             >  2:581660/1‑150 (MQ=255)                                                                                                    ccttcccctccccGTTTT‑GCCTCCCAGGGCCCCAATT‑TTGCCTCCCCCCCC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTCCCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                               <  1:110691/137‑1 (MQ=255)                                                                                                         ccacccccGTTTT‑GCCTGCCAGCGCGCCAATT‑TTGCCAACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                         <  1:210916/143‑1 (MQ=255)                                                                                                             ccccGTTTT‑CCCTGCCGGCGCGCAATTG‑TTGCCTACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGt                                   <  1:384612/149‑1 (MQ=255)                                                                                                             ccccGTTTT‑CCCTGCCACCGCCCCAATG‑TTGCCAACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTg                                     <  2:440412/147‑1 (MQ=255)                                                                                                                cGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                  >  2:585887/1‑149 (MQ=255)                                                                                                                 ggTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                               >  1:343765/1‑151 (MQ=255)                                                                                                                             cAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                    >  2:810689/1‑151 (MQ=255)                                                                                                                                     caagttttcccTACCCCCAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGAt                <  2:732388/139‑1 (MQ=255)                                                                                                                                      aaTG‑TTCCCAACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcg              <  1:289236/148‑1 (MQ=255)                                                                                                                                                tAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc  >  2:923827/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                      AGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGT‑CAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATG‑TTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2578001‑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: ‑