Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	42.4%	82.2 / 30.7	66	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (30/8);  new base A (0/28);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.40e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5012‑5307                                                                                                                                                     |                                                                                                                                                     catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGaa                                                                                                                                                    <  1:232302/150‑1 (MQ=255)  atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAAt                                                                                                                                                   >  2:201572/1‑150 (MQ=255)   tcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATa                                                                                                                                                  <  2:264828/150‑1 (MQ=255)        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                            <  2:11787/151‑1 (MQ=255)        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                            <  2:57628/151‑1 (MQ=255)              tATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCCCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGtt                                                                                                                                      >  2:789172/1‑151 (MQ=255)                        caaaaTCAGGCTGGTTTTGTTTTTAGGCAGCCCTTTTTTTTTAAAAAAAAATTATGGAAAAAAAGGACAGGGAAAAAGAAAAAATTCTCAATAATGGCGTTAACTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAg                                                                                                                             <  2:874782/147‑1 (MQ=255)                           aaTCGGGCCGGGTTTGCTTTTATGAAACCCGGTTTTTTTATGAAGAAATTAGGAAAAAAACTGAAAGGAAAAAAGGAAAAATACCAAATAAAGGCGGTAACTT‑AGGTATTGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCg                                                                                                                               <  2:218141/145‑1 (MQ=255)                           aaTCAGGCGGGTTTTGTTTTTATGCAGCCGGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCACAAAAATGGCGGTAACTT‑AGAAATAGGGTTGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGt                                                                                                                            <  2:345291/148‑1 (MQ=255)                            aTCAGGCGGGTTTTGTTTTTATGCAGCCCGGTTTTTTAATGAAAAAATTAGGGAAAAAAATGAAGGGAAAAAAGGAAAAATTCAAAAAAAAGGCGTTATTTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                        <  2:156413/151‑1 (MQ=255)                                ggcgggttttttttttATGCAGCCCGGTTTTTTTAAAAAAAAATTATGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTA‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                     <  1:915459/138‑1 (MQ=255)                                ggCGGGTTTTGTTTTAATGCGGCCGGGTTTTTTTAAAAAAAAATTAGGAAAAAAAATGACAGAAAAAAAGAAAAAATTCTCAAAAAATGCGGTAACTT‑AGAGATAGGGTTGGGGGAAAATAACACCCGCCGTTCTCATCGAGTAATCTcc                                                                                                                     <  2:378221/150‑1 (MQ=255)                                ggCGGGTTTTGCTTTTATGCACCCCGGTTTTTTTATAAAAAAATTAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAAA‑GGGGTAAATTAAGAAATTAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                    <  1:195260/151‑1 (MQ=255)                                ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                    >  2:959974/1‑151 (MQ=255)                                 gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   >  2:336847/1‑151 (MQ=255)                                      gTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              >  2:972611/1‑151 (MQ=255)                                       tttGTTTTTATGCGCCCCGTTTTTTTTATGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGTTAACTT‑AGAGATTAGGATGGGGGAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  2:860471/145‑1 (MQ=255)                                           cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                         >  1:993964/1‑151 (MQ=255)                                            ttttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                         >  2:697975/1‑150 (MQ=255)                                                      cccGGTTTTTTTAGAAAAAAATTATGAAAAAAAATAACAGGAAAAAAGGAAAAATTCTCAAAAAATGGGATAATTA‑AGAGATTAGAATTGGGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                 <  1:788900/148‑1 (MQ=255)                                                         ggtttttttATGAAGAAATTATGGAGAAAAATGACAGGAAAAAAGGAGAAATTCTAAATAAATGCGGTACCTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                           <  1:485231/148‑1 (MQ=255)                                                          gCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          >  1:832358/1‑151 (MQ=255)                                                           cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          >  1:776187/1‑150 (MQ=255)                                                            ttttttAAAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAATAAGGGCGGTAATTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          <  2:757858/149‑1 (MQ=255)                                                             tttttAAAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATCCTAAAAAAATGCGTTAACTT‑AGAGATTAGGATGGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                        <  2:35863/150‑1 (MQ=255)                                                                ttATGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAg                                                                                     <  1:531521/150‑1 (MQ=255)                                                                 tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    >  2:971125/1‑150 (MQ=255)                                                                  ataaaaaaaTTAGGAAAAAAAATGAAAGGAAAAAAGGAGAAATTCTAAAAAAATGCGGTAACTT‑AGAGATTAGGTTGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAg                                                                                     <  2:529897/142‑1 (MQ=255)                                                                     aaaaaaTTATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTATCTT‑AGAGATTAGGATGGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                   <  2:393534/144‑1 (MQ=255)                                                                     aaaaaaTTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTACCTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                <  1:399469/147‑1 (MQ=255)                                                                     aaaaaaTTAGGAAAAAAAATGACGGGAAAAAAGGAGAAATTCTAAATAAAGGGGGTAATTT‑AGAAATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                <  1:217945/147‑1 (MQ=255)                                                                     aaaaaaTAAGGGAAAAAAATGACAGGAAAAAAGAAGAAATTCTAAAAAAAGGCGTTAACTT‑AGAGATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                <  2:506711/147‑1 (MQ=255)                                                                       gaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                             >  1:643426/1‑151 (MQ=255)                                                                        aatttAGGGAAAAAAAAGACAGGAAAAAAGGAGAAATTCAAAATAAATGGGGTAATTA‑GGAGGTTAGGTGTGCGGAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                                  <  1:122704/142‑1 (MQ=255)                                                                        aaaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                            >  1:189929/1‑151 (MQ=255)                                                                           ttATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         >  1:751004/1‑151 (MQ=255)                                                                            tATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         >  1:701368/1‑150 (MQ=255)                                                                               ggAAAAAAAGAAAAGGAAAAAAGGAAAAATTCTAAATAAATGGGTTAACTT‑AGAGATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        <  2:397219/148‑1 (MQ=255)                                                                                   aaaaaTGACAGGAAAAAAGAAAAAATTCTAAAAAAATGCGTTAACTT‑AGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                 <  1:514186/151‑1 (MQ=255)                                                                                      aaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               >  1:1012350/1‑150 (MQ=255)                                                                                      aaTGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGTTAACTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                              <  2:640813/151‑1 (MQ=255)                                                                                          aCAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  2:102633/1‑151 (MQ=255)                                                                                          aCAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  1:1008821/1‑151 (MQ=255)                                                                                          aCAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  2:493885/1‑151 (MQ=255)                                                                                           cAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  1:681953/1‑150 (MQ=255)                                                                                             ggaaaaaaGGAAAAATTCACAATAAATGCGGTAACTA‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                         <  1:923947/146‑1 (MQ=255)                                                                                              ggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                      >  1:596461/1‑151 (MQ=255)                                                                                                  aaaGGAAAAATTCTCAATAAATGCGGTAACTA‑AAAGATTAGGATGGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                  <  1:774656/151‑1 (MQ=255)                                                                                                         aaaTTCTAAAAAAAGGCGATAATTA‑AGAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                           <  2:412644/151‑1 (MQ=255)                                                                                                            ttCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                        >  1:482895/1‑151 (MQ=255)                                                                                                                cAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                    >  2:977376/1‑151 (MQ=255)                                                                                                                 aaaaaaGGCGGTAACTT‑AGAGATAAGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                      <  2:159129/145‑1 (MQ=255)                                                                                                                    aaaTGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                >  2:313033/1‑151 (MQ=255)                                                                                                                     aaTGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               >  2:756917/1‑151 (MQ=255)                                                                                                                     aaTGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               >  2:856954/1‑151 (MQ=255)                                                                                                                        gCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACGCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                             >  2:388820/1‑150 (MQ=255)                                                                                                                        gCGGTAACTT‑AGAGATTAGGATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                             <  2:529161/150‑1 (MQ=255)                                                                                                                            tAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTc                        >  1:941458/1‑151 (MQ=255)                                                                                                                                tt‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    >  2:368840/1‑151 (MQ=255)                                                                                                                                        ttAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgca                <  1:147938/148‑1 (MQ=255)                                                                                                                                         tAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac            >  2:309988/1‑151 (MQ=255)                                                                                                                                           ggtttGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAATCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg          <  1:232405/148‑1 (MQ=255)                                                                                                                                           ggATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc              <  1:756986/147‑1 (MQ=255)                                                                                                                                              ttGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc        <  1:441802/150‑1 (MQ=255)                                                                                                                                              ttGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt       >  2:24440/1‑151 (MQ=255)                                                                                                                                              tggCGGAGAAAATCATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg          <  2:337220/146‑1 (MQ=255)                                                                                                                                                    agaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  <  1:832207/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5012‑5307

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