Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,661	0	T	C	41.3%	92.4 / 13.6	63	I44V (ATC→GTC)	dppB	dipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (30/7);  new base C (0/26);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTAT  >  NC_000913/3705513‑3705806                                                                                                                                                        |                                                                                                                                                  aCGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATcc                                                                                                                                                 >  1:123293/1‑151 (MQ=255)   gAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATccc                                                                                                                                                >  2:156962/1‑150 (MQ=255)     aCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACg                                                                                                                                             <  1:79700/151‑1 (MQ=255)        ctTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTc                                                                                                                                          >  2:457221/1‑151 (MQ=255)               aaCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTcgcccgcc                                                                                                                                   <  1:540527/151‑1 (MQ=255)                aCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTcgcccgcc                                                                                                                                   >  2:877997/1‑150 (MQ=255)                          ggCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCAc                                                                                                                        >  2:243615/1‑151 (MQ=255)                             ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                     >  2:235488/1‑151 (MQ=255)                             ttttCATTGAAATGCCTAGACCGCCATGCATACCGCCCA‑AAATGTAAGGGAGATACTGCTGCCACATGGTTTTATCTAAGCCGTGTTCGG‑CCAGCAGCTGCGCGGGC‑CGCTCGGGGGAGATCCCCCGTCCGCCCGCCATGATCATCACCg                                                                                                                      <  2:397935/150‑1 (MQ=255)                              tttCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCg                                                                                                                      >  1:371383/1‑149 (MQ=255)                                 cATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATc                                                                                                                  >  2:195541/1‑150 (MQ=255)                                 atttGAAGTGCTTAGATCGCCAT‑CAAACCCCCCCAAAAGGAAATGGAGATCCTCCTGCCCCATCGTTTTACCTAGGCCAATTCCAC‑CCAGCGTCTGGGCGGGC‑CGTTCGGGGGAGCCCCCCCTTTGGCCCGCCATGATCATCACCGGATc                                                                                                                  <  1:624130/148‑1 (MQ=255)                                      aaatccttagcccGCCATCCATAACCCCCA‑AAATGTAATGGAGATCCTGCCCCCACATCGGTTTATCTAAGCCGATTCCAC‑CCAGCAGCTGCGCGTGG‑CGCTCGGGGGAGATCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCg                                                                                                             <  1:90280/138‑1 (MQ=255)                                      aaaTGCGTAGCCCGCCGTGCATACCCCCCA‑AAATGTAAGGGAGATACTGCTGCCACTTCGTTTTATCT‑AGCCGGATTCCGCCCACCAGCGGCGCGTGA‑CGTTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATc                                                                                                                  <  1:5196/145‑1 (MQ=255)                                      aaaTGCCTAGATCGCGATGCATAACCCCCA‑AAATTTAATGGAGATACTGCTGCCCCATCGTTTTATCAAGCCCGATTTCAG‑CCAGCGGCGGCGGTTGC‑CGCTCTGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCg                                                                                                             <  2:390422/150‑1 (MQ=255)                                      aaaTGCCTAATCCCCCAGGCATAACCCCCA‑AAAGGTAAGGGAAATCCTGCTCCCACATCGTTTTATCTACGCCGAGTTCAG‑CCAGCAGCTGCGCGTGC‑CGCTTTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCg                                                                                                             <  2:404067/150‑1 (MQ=255)                                       aaTGCCTAGATCGCCAGGCAAAACCCCCA‑AAATGTAAGGGAGATCCTGCTGCCACTTGGGTTTATCAAAGCCGAGTTCAG‑CCAGCAGGTGCGCGTGA‑CGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                           <  2:596036/151‑1 (MQ=255)                                        acccgcagaacacgaTGCATAACTCCAA‑AAATGTACGGCACATCCTCCAGCCTCTTGGGTTGACCGAAGCCATATTCAC‑CCACCAGTTGCG‑GCGATCGCTCTGGGGATACCCCACGTTGCCCCGCCATGATCATCACCGGATCGCCCGGGa                                                                                                          <  2:702563/137‑1 (MQ=255)                                            cTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                      >  2:621327/1‑151 (MQ=255)                                              aGACCGCCATAAAAACCGCCCA‑AAATGAAGGGGAGATACTGCCGCCACGTCGTTTTATCTACGCCGAGTCCAG‑CCAGCAGCTGCGCGGGC‑CGCTGGGGGGAGATCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtg                                                                                                     <  1:943241/150‑1 (MQ=255)                                                           aacccccca‑AAATGAAATGGAGATACTGCTGCCACTTCGTTTTATCAACCCCGATTTCAC‑CCGGCAGCTGCGCGGGC‑CGTTGGGGGGAGATCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                                           <  1:247605/142‑1 (MQ=255)                                                               gCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                   >  1:724102/1‑151 (MQ=255)                                                                ccca‑AAATGTAATGGAGATCCTGCTCCCACATCGGTTTATCTAAGCGGAGTTCAG‑CCAGCAGCGGCGCGTGA‑CGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                                  <  2:391998/151‑1 (MQ=255)                                                                ccca‑AAATGAAAGGGAGATCCTGCTCCCACAGTGGTTTATATAAGCCGATTCCAG‑CCGGCAGCGGCGGGGGC‑CGTTCGGGGGAAACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAt                                                                                      <  1:656617/147‑1 (MQ=255)                                                                      aaTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                             >  1:313890/1‑151 (MQ=255)                                                                           aagggAAATACTGCCCCCACTTCGTTTTATCAAAGCCGATTCCAC‑CCAGCAGCGGCGCGGGC‑CGCTGGGGGGAGACCCCACGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                            <  1:457211/144‑1 (MQ=255)                                                                              ggAGATCCTGCTCCCACATCGTTTTATCAAAGCGGAGTTCGG‑CCAGCAGCTGCGGGTGC‑CGCTTGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAcc                                                                     <  2:559780/151‑1 (MQ=255)                                                                                agaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                    >  1:475756/1‑150 (MQ=255)                                                                                 gaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                  >  1:48935/1‑151 (MQ=255)                                                                                  aTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                  >  1:751329/1‑150 (MQ=255)                                                                                   tACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                                 >  1:916997/1‑150 (MQ=255)                                                                                   tACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                                 >  1:734133/1‑150 (MQ=255)                                                                                        ctgcAACATCGTTTTATCAAAGCCGAGTCCAG‑CCAGCGGCTGCGGGTGC‑CGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACg                                                             <  2:109562/149‑1 (MQ=255)                                                                                         tgcCACATCGTTTTATCAAAGCCGATTCCCG‑CCGGCAGCTGCGCGGGC‑CGCTCTGGGGAGACCCCCCGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                <  2:627425/145‑1 (MQ=255)                                                                                         tgcCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg                                                          >  1:915320/1‑151 (MQ=255)                                                                                           ccACATGGTTTTACCAAACCCGAGTCCCG‑CCGCCAGCGGCGCGGGC‑CGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                           <  1:323801/148‑1 (MQ=255)                                                                                              caTCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCCTGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                                     >  1:552880/1‑151 (MQ=255)                                                                                              caTCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                      >  1:317904/1‑150 (MQ=255)                                                                                                tCGTTTTATCTAAGCCGAGTTCAG‑CCAGCAGCGGCGCGTGA‑CGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                      <  2:857912/148‑1 (MQ=255)                                                                                                    tttATCTAAGCCGTGTCCAG‑CCAGCCGCTGCTGGGGG‑CGCTCTGGGGAGACCCCCCTTTCGACCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                                     <  2:773149/145‑1 (MQ=255)                                                                                                    tttATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGt                                               >  1:924853/1‑151 (MQ=255)                                                                                                    tttATCTAAGCCGAGTCCAG‑CCAGCGGCGGCGGGTGC‑CGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAg                                                <  2:35584/150‑1 (MQ=255)                                                                                                    tttATCAAAGCCGATTCCAG‑CGAGCGGCGGCGCGGGC‑CGCTCGGGGGAGACCCCGCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGAc                                                   <  2:637292/147‑1 (MQ=255)                                                                                                     ttATCAAACCCGATTACAC‑CCCGCGGCCGGGCGTGC‑CGCTGGGGGGAGCCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGa                                                    <  1:566516/145‑1 (MQ=255)                                                                                                           aaGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                        >  1:252030/1‑151 (MQ=255)                                                                                                              cgggttccag‑ccagcAGCGGCGGGGCC‑CGCGTGGGGGAGCCCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACg                                       <  1:877979/142‑1 (MQ=255)                                                                                                              ccGAGTCCAG‑CCAGCAGCGGCGCGTGA‑CGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                                     <  2:451068/151‑1 (MQ=255)                                                                                                              ccGAGTACAG‑CCACCAGCGGGGTGAGA‑GGCGTGGGGAACACCCCACGGCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCaaa                                         <  1:156964/147‑1 (MQ=255)                                                                                                               cGAGTTCAG‑CCGCCGGCTGCGGGGGC‑CGCTGGGGGGAGATCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                        <  1:243611/147‑1 (MQ=255)                                                                                                               cGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCg                                    >  2:212423/1‑151 (MQ=255)                                                                                                                gatttCAG‑CCAGCAGCTGGGCGTGC‑CGCTTGGGGGAGACCCCCCGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                   <  2:166894/148‑1 (MQ=255)                                                                                                                cgggtccgcccggcagcGGCGCGGGA‑CGTGGGGGGGAGCCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc                                             <  1:235484/130‑1 (MQ=255)                                                                                                                   ttCAG‑CCAGCAGCTGGGGGTGC‑CGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                   <  2:673348/148‑1 (MQ=255)                                                                                                                         ccagcAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaa                            >  2:432477/1‑150 (MQ=255)                                                                                                                         ccagcAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaa                            >  2:781698/1‑150 (MQ=255)                                                                                                                             cggcTGCGGGTGA‑CGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGc                       <  2:334762/149‑1 (MQ=255)                                                                                                                             cagcTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGc                       >  1:39778/1‑151 (MQ=255)                                                                                                                                 tGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc                    >  1:446591/1‑150 (MQ=255)                                                                                                                                  gcgcGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACAt                  >  1:345783/1‑151 (MQ=255)                                                                                                                                             ctctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAtt        >  1:167597/1‑151 (MQ=255)                                                                                                                                                 ggggAGACCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGaa                               <  1:35346/124‑1 (MQ=255)                                                                                                                                                   ggAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtat  >  1:383245/1‑151 (MQ=255)                                                                                                                                                     agaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtat  >  1:342175/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                  ACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCC‑AAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAG‑CCAGCAGCTGCGCGTGA‑CGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTAT  >  NC_000913/3705513‑3705806

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