Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,571,036	0	T	G	42.2%	104.7 / 10.7	65	N239T (AAC→ACC)	eutE	aldehyde oxidoreductase, ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base T (8/29);  new base G (27/0);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 4.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACC  >  NC_000913/2570886‑2571185                                                                                                                                                       |                                                                                                                                                       aTCAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGt                                                                                                                                                        <  1:194468/151‑1 (MQ=255)   cAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGt                                                                                                                                                        <  2:75096/149‑1 (MQ=255)    aGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGc                                                                                                                                                     <  1:543777/151‑1 (MQ=255)         cATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCAGCCACCACTACCGGCGGGTTGCCAgcg                                                                                                                                                <  2:34132/151‑1 (MQ=255)               ttCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGcc                                                                                                                                              >  1:908035/1‑147 (MQ=255)                    cGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAAt                                                                                                                                     <  2:401414/151‑1 (MQ=255)                    cGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCCGCGGCaa                                                                                                                                      >  1:4420/1‑150 (MQ=255)                    cGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTTCgg                                                                                                                                         >  2:770016/1‑147 (MQ=255)                     ggCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATc                                                                                                                                    <  1:482838/151‑1 (MQ=255)                              aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCCGCGCCCGCGGGCATCAGACGgtt                                                                                                                            >  1:491307/1‑150 (MQ=255)                                    aacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCCATCCGACGTTTAtttggg                                                                                                                     >  2:412629/1‑147 (MQ=255)                                     acaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGgtgt                                                                                                                    <  1:311080/151‑1 (MQ=255)                                     acaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCCGCGCCCGCGGCCATCacccggttttttg                                                                                                                       >  2:583371/1‑136 (MQ=255)                                        aTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCCATAAGCCGATTTCTGgtgggt                                                                                                                  >  1:516876/1‑150 (MQ=255)                                        aTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAAATAGACGGTTATTGGGGTGtt                                                                                                                 >  1:267828/1‑151 (MQ=255)                                          cAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGttt                                                                                                                <  2:901053/150‑1 (MQ=255)                                                cTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGACGGTTATTGGGGTGTTTACGcgcc                                                                                                         >  1:354886/1‑151 (MQ=255)                                                 ttttCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTTCCCGCGCCTGCCGCGACCAGGGGTTTATTGgtgctttt                                                                                                                >  1:216278/1‑138 (MQ=255)                                                  tttCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGGGTGTTTACGcgccgc                                                                                                       >  1:715563/1‑151 (MQ=255)                                                  tttCGTCGGCACAAATGATGTTATTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGcgcc                                                                                                         <  1:429323/149‑1 (MQ=255)                                                     cgtcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGACGGCGGTTTCATCCACCACTACCGGACGGATGCCCGGGCCTTGCGCAAACAGACGTTTATTGGGGTGGTTTCGcgccg                                                                                                        >  1:941524/1‑147 (MQ=255)                                                      gtcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTct                                                                                                   <  1:277356/151‑1 (MQ=255)                                                       tcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCGTGCGGGAATCAGGCGTTTATTTGGGTGTTTACGcgcc                                                                                                         >  2:156517/1‑144 (MQ=255)                                                       tcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCAATCAGAAGGTTTTTTGGGGGTTTTCGcgccgcct                                                                                                     >  2:551566/1‑146 (MQ=255)                                                        cgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTctac                                                                                                 <  1:891472/151‑1 (MQ=255)                                                        cgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTGcta                                                                                                  <  2:363960/150‑1 (MQ=255)                                                          gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTctacta                                                                                               <  2:187558/151‑1 (MQ=255)                                                          gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTGTCCAGCGCCCGCGGCCATCCGACGTTTATTTGGGTGGTTTCGCGCCGCt                                                                                                      >  2:410169/1‑144 (MQ=255)                                                           cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGACGTTTTTTGGTGGGTTTACGCGCCGCTTctacta                                                                                               >  1:268208/1‑150 (MQ=255)                                                              aaaTGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGc                                                                                           <  2:476210/151‑1 (MQ=255)                                                              aaaTGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTCCGGCCATCAAACGGTTTTTTGGGTGTTTAAGCGCCCCTTctaata                                                                                               >  2:563401/1‑147 (MQ=255)                                                                 tgatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGGGCCAGCGCCCGCGGGGAACAGGCGGTTATTTGGGTGGTTTCGGGCCGCCTctacctcc                                                                                             >  2:21255/1‑142 (MQ=255)                                                                  gatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGAGCCTTGGGCAATCAGAAGGTTTTTTGGGGGTTTACGGGCCGCGTCTACGACCCCttt                                                                                        >  2:307164/1‑149 (MQ=255)                                                                       tgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGCGCCCGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCgccgcc                                                                                  >  1:785196/1‑151 (MQ=255)                                                                       tgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCCGCGCCTGCGGCAATCCGACGTTTTTTGGGGTGTTTACGCGCCGCTTCTACTACCGCTTCgccgcc                                                                                  >  1:570712/1‑151 (MQ=255)                                                                       tgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAACGCCTGCGGCAATCAGACGGTTATTGGGGTGTTTACGCCCCCCTTCTACTACCCCTTCcccccc                                                                                  >  1:733578/1‑151 (MQ=255)                                                                               aaaGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCACGCCTGCGGCAATCAGACGTTTTTTTGGGTGTTTACGCGCCGCTTCTACTACCGCCTCGCCGCC‑GGGTacc                                                                          >  2:672835/1‑150 (MQ=255)                                                                                    agCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGGTTATTGGGGTGTTTACGCGCCGCTTCTACTACCCCTTCGCCGCC‑GGGTaccacc                                                                       >  1:718269/1‑148 (MQ=255)                                                                                        ccTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAAACAGACGTTTATTGGGGGGGTTACGCGCCGCCTCTACTACCGCCTCGCCGCC‑GGTTaccacc                                                                       >  1:470671/1‑144 (MQ=255)                                                                                             gACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGGCATCAGACGGTTATTGGGGGGTTTTCGCGCCGCCTCTACCACCGCCTCCCCGCC‑GGGTACCACCACCAGGcc                                                               >  2:379853/1‑147 (MQ=255)                                                                                                gATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATAccc                                                        <  1:272126/151‑1 (MQ=255)                                                                                                gATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGGCGCGAACAGACGTTTTTTTGGGTGGTTACGCGCCGCTTCTTCTACCCCTTCCGCGCC‑GGGTTCCACCCACAGGCCGGTAccc                                                        >  2:66482/1‑151 (MQ=255)                                                                                                   ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGa                                                     <  1:663796/151‑1 (MQ=255)                                                                                                   ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGGGGGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGGTACCACCAGCAGGCCGATAccccg                                                      >  1:459520/1‑148 (MQ=255)                                                                                                    gACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGaa                                                    <  2:667091/151‑1 (MQ=255)                                                                                                     aCTGAGCGGCACGGGCGAGGTCGGCGGATTCATCCACCACTACCGGCGGGGTGCCAGGGCCTGGGGCCATCCGGCGTTTATTGGTGGGTTTACGCGCCCCCTCTCCTACCGCCTCGCCGCC‑GGTTAACCACAACAGGCCGATACCCCGaa                                                    >  2:722721/1‑150 (MQ=255)                                                                                                        gAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACtt                                                <  1:768196/151‑1 (MQ=255)                                                                                                        gAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACtt                                                <  1:554145/151‑1 (MQ=255)                                                                                                        gAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACtt                                                <  2:778961/151‑1 (MQ=255)                                                                                                         aGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCCATCACACGGTTTTTTGTGTGTTTTCGCGCCGCCTCCTCTCCCCCCTTCCCGCC‑GGTTTCCCCCCCCCAGGCGATACCCGGGAAAtt                                                >  2:468855/1‑150 (MQ=255)                                                                                                                cGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGc                                        <  1:115521/151‑1 (MQ=255)                                                                                                                cGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCGATCCGGCGGTTATTTGGGTGTTTACGCGCCGCCTCTAATACCGCCTCGCCGCC‑GGGTACCACCCGCAGGCCCCTACCCGGGAACTTTAACa                                           >  2:817204/1‑148 (MQ=255)                                                                                                                  ggCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGt                                      <  1:847415/151‑1 (MQ=255)                                                                                                                  ggCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGGCCGCGGGCAACCGACGGTTTTTGGTGGGTTTTCGCGCCGCCTCTACTACCGCCTTGCCCCCGGGTTACCA‑CAACAAGCCCATACCCCGGAACTTTAACAAGCg                                       >  1:475880/1‑150 (MQ=255)                                                                                                                      aGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCCGACGGTTTTTGGGGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAACGTTg                                    >  1:431234/1‑149 (MQ=255)                                                                                                                      aGGTCGGCGGTTTCATCCACCACTACCGGAGGGGTGCCAGCGCCTTTGGTAAACACACGTTTATTTGTGTGTTTAAGCCCCGCTTCTTCTACCCCGTCGCCCCC‑CGGTTCCACCAGCAGGCCCATAACCGGAAACTTGAACAAGCGt                                      >  1:628678/1‑147 (MQ=255)                                                                                                                            gcggTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGttt                            <  2:264612/151‑1 (MQ=255)                                                                                                                             cggTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTc                           <  1:266107/151‑1 (MQ=255)                                                                                                                                   cATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGatat                      <  1:766735/150‑1 (MQ=255)                                                                                                                                      ccaccaCTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCgg                  <  1:352022/151‑1 (MQ=255)                                                                                                                                         ccacTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTCCCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAtt               <  2:786280/151‑1 (MQ=255)                                                                                                                                         ccacTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAtt               <  1:222074/151‑1 (MQ=255)                                                                                                                                         ccacTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAtt               <  1:397187/151‑1 (MQ=255)                                                                                                                                          cacTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAttt              <  1:873336/151‑1 (MQ=255)                                                                                                                                            cTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGACGTTTATTTGGGGGTTTTCGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGGTACCACCCGCCGGGCGATAACCGGGAACTTGAAAAAGCGGTGCGCGGGTTCGGTATCCGGGTTTg             >  2:312244/1‑150 (MQ=255)                                                                                                                                                cggcggGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCcaca        <  2:338589/151‑1 (MQ=255)                                                                                                                                                    ggggTGCCAGCGCCGGCGGCAATCAGACGTTTATTGGGGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCGGGCCGATACCCGGAAACTTGAACAAGCGGTGCGCGGTTTCGAAATCCGGATTTGCCACAg       >  2:195450/1‑148 (MQ=255)                                                                                                                                                      gTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCACCAGCCCGATACCCGGCAACTGGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAAcc  >  2:891602/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       ATCAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCC‑GGTTACCACCAGCAGGCCGATACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACC  >  NC_000913/2570886‑2571185

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