Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	44.3%	82.4 / 20.9	61	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (5/29);  new base G (27/0);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 2.02e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508508‑3508805                                                                                                                                                      |                                                                                                                                                       cAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaa                                                                                                                                                       <  2:579214/151‑1 (MQ=255)   gTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaaa                                                                                                                                                     <  1:264797/151‑1 (MQ=255)    tCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGgaaa                                                                                                                                                     >  2:918654/1‑150 (MQ=255)      cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaaaa                                                                                                                                                    >  2:55382/1‑149 (MQ=255)         gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGgaaaaat                                                                                                                                                  >  2:1071673/1‑146 (MQ=255)             tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATcaca                                                                                                                                              >  2:596045/1‑148 (MQ=255)                  aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                       <  2:931730/150‑1 (MQ=255)                  aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                      <  2:1076782/151‑1 (MQ=255)                  aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                      <  2:668704/151‑1 (MQ=255)                      ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGGAGcccc                                                                                                                                  >  2:1006305/1‑151 (MQ=255)                          agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcaca                                                                                                                              <  1:1002243/151‑1 (MQ=255)                              aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGAGCGCGGGAAAATTTAAAAACCGTAGACCCcccacc                                                                                                                            >  2:1125039/1‑147 (MQ=255)                                cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATC‑AGa                                                                                                                        >  2:600070/1‑150 (MQ=255)                                 aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACAAACGGTAAccccccc                                                                                                                               >  1:925507/1‑141 (MQ=255)                                 aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATC‑AGaa                                                                                                                       >  2:850074/1‑150 (MQ=255)                                        aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCgtggt                                                                                                               <  2:502961/151‑1 (MQ=255)                                         gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCgtggtg                                                                                                              <  1:307072/151‑1 (MQ=255)                                                aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATC‑AGAAAAAGTGGGGCGCacac                                                                                                       >  2:339543/1‑151 (MQ=255)                                                aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGTAGCCCCCCCATC‑AAAAAACGTGGTgcg                                                                                                            >  2:484537/1‑146 (MQ=255)                                                 tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAACCCCCCCCTC‑AAAAAACGGGGTGGGCACAcc                                                                                                      >  1:938313/1‑151 (MQ=255)                                                   ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACGGTAGCCCCCACCTC‑AGAAAGCGGGGGGGGCCCAcc                                                                                                      >  2:250571/1‑149 (MQ=255)                                                   ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACCGGAGCCCCCCCATC‑AGAAAGCGGGGGGGGCACAcc                                                                                                      >  1:264151/1‑149 (MQ=255)                                                      gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGAAAAATTCACACACGGTAGCCCCCCCATC‑AGAAAACGGGGGGGGCACACCCAGCt                                                                                                 >  2:484076/1‑151 (MQ=255)                                                            gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACAGGAGACCCCCCATC‑CGAAAACGTGGGGGGCACACCGCGCCGCGt                                                                                             >  2:505005/1‑149 (MQ=255)                                                            gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCATC‑AGAAAACGGGGGGCGCACACCCAGCTg                                                                                                >  2:674524/1‑146 (MQ=255)                                                              aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCATC‑AGAAAGCGTGGGGGGCACACCCAGCTGCGTa                                                                                            >  2:359852/1‑148 (MQ=255)                                                              aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                         <  2:681072/151‑1 (MQ=255)                                                              aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                         <  2:357443/151‑1 (MQ=255)                                                                 ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                      <  2:196935/151‑1 (MQ=255)                                                                    tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATC‑AGAAAGCGGGGTGCGCACACCCAGCTGCGTAGCGatc‑a                                                                                    >  2:1068123/1‑149 (MQ=255)                                                                      ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACAAGGTAGCCCCCCCATC‑AGAAAACGTGGGGCGCACACCGCGCTGCGTAGCGatc‑at                                                                                   >  2:532714/1‑148 (MQ=255)                                                                        aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCtt                                                                              <  2:938384/151‑1 (MQ=255)                                                                          cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTg                                                                             <  1:377161/150‑1 (MQ=255)                                                                             aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATC‑AGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATC‑ATccc                                                                                >  2:923281/1‑144 (MQ=255)                                                                             aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGGTGa                                                                         <  1:305638/151‑1 (MQ=255)                                                                              aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGGAGCCCCCCCATC‑AGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATC‑ATCCCCTTCTGAt                                                                        >  2:511647/1‑151 (MQ=255)                                                                              aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCCTC‑CGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATC‑CTCCCCTTCTGAt                                                                        >  2:269051/1‑151 (MQ=255)                                                                                tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAAGGTAGCCCCCCCATC‑AGAAAAAGTGGTGCGCACACCCCGCTGCGTTGCGGTC‑ATTCCctgctgtt                                                                        >  1:459028/1‑147 (MQ=255)                                                                                 cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATaaa                                                                     <  2:1032729/151‑1 (MQ=255)                                                                                 cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATaa                                                                      <  2:1089053/150‑1 (MQ=255)                                                                                    aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCa                                                                  <  2:82258/151‑1 (MQ=255)                                                                                    aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCa                                                                  <  2:481893/151‑1 (MQ=255)                                                                                      ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGc                                                                <  1:6359/151‑1 (MQ=255)                                                                                                 cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGGGGGGGAAAAATTACACACCGCAGGCCCCCCATC‑AGGAAACGGTGGGCGCACACCGGGCCCTGTTGGCAAC‑ATCCCTTGCTGATaaaacccaactgctacc                                                       >  2:258719/1‑135 (MQ=255)                                                                                                  aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATC‑AGAAAGCGTGGGGGGCACACCCAGCTGCGGAGCGATC‑ATCCCTTGCTGATAAAACAGCACTGGccccac                                                     >  1:915504/1‑150 (MQ=255)                                                                                                      cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                  <  2:455469/149‑1 (MQ=255)                                                                                                       ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACCGTAGCCCCCACATC‑AGAAAGAGTGGGGCGCACACCGCGCTGCGTAGCGAGC‑ATCCCTTGCTGATAAAACAACACCgg                                                           >  2:483812/1‑139 (MQ=255)                                                                                                           ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGGAAAAATTACAAACAGGAGCCCCCCCATCAAAAAAGCG‑GGGGCGCACACCCAACCGCGTAGCGGTCAATCCCTGCCTTTTAAACCA‑CACTGGCACCACCTCCCggggt                                           >  1:455749/1‑149 (MQ=255)                                                                                                               gaaagaaaAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACGGGAGCCCCCACAAC‑AGAAAACGGGGGGCGCACACCCAGCCGCGGAGCGATC‑AACCCTTTCTGATAAAACAACACCGGCACCCCTTCCCGGGGAtgg                                        >  2:525520/1‑150 (MQ=255)                                                                                                               gaaagaaaAGAGCGACCATTTTGTGGCGGAGGGGGGCGGGAAAAATCACAAACGGGAGGCCCCCCATC‑AGAAAAAGGGGGGGGCACAACGAAGTGGGGAGCGGTC‑AAACCCTTCTGATAAAACAGGACCGGGACCCCTTCCCGGGGAtgg                                        >  2:560295/1‑150 (MQ=255)                                                                                                                aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                      <  2:1042585/151‑1 (MQ=255)                                                                                                                  agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACGGGAGCCCCCACATC‑AGAAAACGTGGGGGGCACACCCAGCTGCGGAGCGATC‑ATCCCTTGCTGATAAAACAGCACTGGCCCCACCTCGCGGGGAtggggg                                     >  1:646288/1‑150 (MQ=255)                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                   <  2:421486/151‑1 (MQ=255)                                                                                                                     aaagaGCGACCATTTTGTGGCGGTTGCGCGCGGAAAAATTCACACACGGGAGCCCCCCCATC‑AGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGATC‑ATCCCTTTCTGCTAAAACAGCACTGGCACCACCTCGCGGGTAtggggg                                     >  1:54061/1‑147 (MQ=255)                                                                                                                      aagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAACCCCCACATC‑AGAAAGAGTGGTGCGCACACCCAGCTGCGTAACGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACCTCGCGGGTATGGGGGgtgt                                 >  2:904174/1‑150 (MQ=255)                                                                                                                      aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                <  1:703202/151‑1 (MQ=255)                                                                                                                                aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                      <  1:1068027/151‑1 (MQ=255)                                                                                                                                     gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                  <  2:697933/150‑1 (MQ=255)                                                                                                                                     gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                 <  1:320787/151‑1 (MQ=255)                                                                                                                                       ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                <  2:870842/150‑1 (MQ=255)                                                                                                                                           gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc           <  1:46036/151‑1 (MQ=255)                                                                                                                                             ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt         <  2:645454/151‑1 (MQ=255)                                                                                                                                               cgcgcgGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg        >  2:812952/1‑150 (MQ=255)                                                                                                                                                gcgcgGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt      <  2:164651/151‑1 (MQ=255)                                                                                                                                                gcgcgGAAAAAGTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt      <  1:527435/151‑1 (MQ=255)                                                                                                                                                    ggAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc  >  2:230578/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATC‑AGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATC‑ATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508508‑3508805

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