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

CAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTA  >  NC_000913/2570888‑2571182                                                                                                                                                     |                                                                                                                                                      cAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGtt                                                                                                                                                      <  1:784498/150‑1 (MQ=255)      gCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGcc                                                                                                                                                   >  1:785733/1‑148 (MQ=255)       caccAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAgcg                                                                                                                                               <  2:496218/148‑1 (MQ=255)          cAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGcc                                                                                                                                             >  1:328326/1‑150 (MQ=255)           aGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGcc                                                                                                                                             >  1:376225/1‑149 (MQ=255)               cATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGGAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCgg                                                                                                                                        <  1:456610/150‑1 (MQ=255)               cATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGc                                                                                                                                       <  2:1020862/151‑1 (MQ=255)                aTCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACATATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCgg                                                                                                                                        >  2:995647/1‑149 (MQ=255)                   ggCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATc                                                                                                                                   <  2:556979/151‑1 (MQ=255)                       aCGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAg                                                                                                                                 <  2:908973/149‑1 (MQ=255)                       aCGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCAGGGCCTGCGGCTAACAGa                                                                                                                                >  1:610266/1‑150 (MQ=255)                        cGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑g                                                                                                                             <  2:498002/151‑1 (MQ=255)                          cTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGAC‑ggt                                                                                                                           >  1:228195/1‑149 (MQ=255)                            aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCAATAAGAC‑Gttt                                                                                                                          >  2:464896/1‑150 (MQ=255)                               aacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACTGGCGGGGTGCCCGCGCCCGCGGGAATCAGAC‑GTTTAtt                                                                                                                       >  1:651453/1‑150 (MQ=255)                               aacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTg                                                                                                                      <  1:62915/151‑1 (MQ=255)                                  aacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCCAACAGAC‑GGTTATTGgtg                                                                                                                   >  1:657463/1‑151 (MQ=255)                                   acaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCAGCGCCCGCGGCAATCAGAC‑GTTTATTtgtg                                                                                                                   >  1:748900/1‑150 (MQ=255)                                    caaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCCACGCCTGCGGCAATCAGACGGTTTTTTggg                                                                                                                    >  2:780544/1‑148 (MQ=255)                                        cAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCACGCCCGCGGCAATCAGAC‑GTTTATTTGTGGGttt                                                                                                              >  1:461277/1‑150 (MQ=255)                                         aGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGGGGGGTGGCAGCGCCTGCGGCAATCACAC‑GGTTTTTTgggtgtgt                                                                                                              >  2:169544/1‑147 (MQ=255)                                         aGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTAc                                                                                                            <  1:99252/151‑1 (MQ=255)                                           tACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGGGGCAATCAGAG‑GTTTTTTgggggggtt                                                                                                              >  2:491091/1‑140 (MQ=255)                                            aCCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTAcg                                                                                                           <  2:15458/149‑1 (MQ=255)                                             ccTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACgcgc                                                                                                        <  1:750653/151‑1 (MQ=255)                                             ccTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGGGCCTGCGGGAATTAGGC‑GGTTATTTGGGTGttttaccgc                                                                                                        >  2:393693/1‑145 (MQ=255)                                               ttttCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCCGCGGGAATCAGAC‑GTTTATTGGTGGGTTTACGcgcc                                                                                                       >  1:1000490/1‑150 (MQ=255)                                                tttCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGcgccg                                                                                                      <  2:117273/150‑1 (MQ=255)                                                   cgtcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCCATCAGAC‑GGTTATTGGGGTGTTTTAGCGCCGCtt                                                                                                   >  2:704427/1‑150 (MQ=255)                                                         cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCTCCGCCATCAGGCGGTTTATAGGGGTGTTTaaccgcgcctcta                                                                                                 >  2:294330/1‑135 (MQ=255)                                                             aaTGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCt                                                                                        <  2:313474/151‑1 (MQ=255)                                                               tgatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGCGCCTGCGGCAATTAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGcct                                                                                       >  1:160472/1‑148 (MQ=255)                                                               tgatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCGATCAGAA‑GGTTATTGGGGGGTTTTCGCGCCGCCTCTACTAACGCt                                                                                        >  2:614443/1‑149 (MQ=255)                                                                   gttgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCAGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTcgcc                                                                                   <  2:431064/150‑1 (MQ=255)                                                                          tCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTGGGCCGCGCCCGCGGCAATCAGAC‑GATTATTTGGGTGTTTTCCCGCCGCTTCCAATAACGCCTTGCCC‑CCGGtt                                                                           >  1:178604/1‑150 (MQ=255)                                                                           cGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACCACCGGCGGGATACCAGCGCCCGCTGCGATCCGCA‑GTTTTTTGGTGTGTTTATCCGCCTCGTTTATGTCCACTTCCCCG‑CCggctcc                                                                         >  2:804605/1‑147 (MQ=255)                                                                            gAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTacc                                                                        <  1:919796/151‑1 (MQ=255)                                                                            gAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGAC‑GTTTTTTGGGGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTacc                                                                        >  1:717918/1‑151 (MQ=255)                                                                              aagaagCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGGGCCTGCGGCCAACAGAC‑GGTTTTTGGTGTGGTTTAGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGGTacccc                                                                      >  1:55881/1‑149 (MQ=255)                                                                               agaagCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCAATCAGAC‑GTTTTTTTGGGTGTTTTCGCGCCCCTTCTACTACCGCTTCGCCC‑CCGGGTaccacc                                                                     >  1:361725/1‑151 (MQ=255)                                                                                    cgcCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCACCCCCTGCGGGCATCACAC‑GGTTTTTTGGGGGTTTTCGCGCCGCTTCCTCTAACGGTTTGCCGCCCGGTTacacccac                                                                   >  2:814132/1‑143 (MQ=255)                                                                                      ccTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCAATCAGAC‑GGTTTTTGGGGTGTTTTCGCGCCGCTTCTACTACCGCCTCGCCG‑CCGGGTACCACCAGCAgg                                                               >  1:131739/1‑150 (MQ=255)                                                                                       cTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCCATCAGAC‑GTTTTTTGGGGGGTTTTCGCGCCGCTTCTTCTACCCCCTCGCCC‑CCCGGTTCCCCCCGCAgg                                                               >  2:584034/1‑149 (MQ=255)                                                                                            aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGTCAGCGCCCGCGGCCATCAGAC‑GTTTTTTTGGGTGTTTTCGCCCCCCTTTTTCTTCCCCCTCCCCG‑CCGGTTAcccccacgcag                                                               >  1:181378/1‑139 (MQ=255)                                                                                              gATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑Accc                                                     <  1:336072/151‑1 (MQ=255)                                                                                                 ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGCCAATCAGTC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCg                                                    <  2:485697/149‑1 (MQ=255)                                                                                                    cTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTTTTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGGTACCACCAGCAGGCCGAT‑ACCCGGaaa                                                >  1:995530/1‑150 (MQ=255)                                                                                                    cTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGGGCCTGCGGCGATCAGCC‑GGTTATTTGTGTGTTTTTGCCGCGCCTCCACTACCCCCTCCCCC‑GCCGTTACCACCCGCCGGcccgtttcccgaaa                                                 >  2:220206/1‑138 (MQ=255)                                                                                                     tGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGGGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGaaa                                                >  1:681520/1‑149 (MQ=255)                                                                                                     tGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCAATCAGAC‑GTTTTTTTGGGGGTTTACGCGCCGCCTCCACTACCCCCTCGCCG‑CCGGTTACCACCAGCAGGGCCATAACCCGGaa                                                 >  2:199768/1‑149 (MQ=255)                                                                                                         cGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCAAATAGAC‑GTTTTTTTGGGTGGTTACGCGCCGCCTCTTCTACCGCCTCCCCC‑CCCGGTACCACCCGCCGGGCGAT‑ACCCGGgaa                                                >  2:148240/1‑145 (MQ=255)                                                                                                              cGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAg                                      <  1:93197/150‑1 (MQ=255)                                                                                                              cGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCCATCCGAC‑CGTTATTGGTGTGTTTTCGCGCCGCTTCTACTTCCGCTTCGCCG‑CCGGTTACCCCCAGCAGGCCCAT‑ACCCGGGAACTTGAACAAGc                                     >  1:680860/1‑151 (MQ=255)                                                                                                                  cGAGGTCGGCGGTTTCATCCGCCACCACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGCC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGtt                                  <  2:400298/150‑1 (MQ=255)                                                                                                                    aGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTgcg                               <  1:97465/151‑1 (MQ=255)                                                                                                                           cggTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTc                        <  1:4975/151‑1 (MQ=255)                                                                                                                             gTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGa                      <  1:451877/151‑1 (MQ=255)                                                                                                                                tCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGata                    <  1:1042762/150‑1 (MQ=255)                                                                                                                                 cATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGatat                   <  1:826501/150‑1 (MQ=255)                                                                                                                                 cATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATc                  <  1:638976/151‑1 (MQ=255)                                                                                                                                  aTCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATcc                 >  2:546608/1‑151 (MQ=255)                                                                                                                                    ccaccaCTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCgg               <  1:192370/151‑1 (MQ=255)                                                                                                                                    ccaccaCTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCg                <  1:447320/150‑1 (MQ=255)                                                                                                                                       ccacTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAt             <  2:886962/150‑1 (MQ=255)                                                                                                                                         acTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTg          <  2:354889/151‑1 (MQ=255)                                                                                                                                               ggcggGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAg    <  1:364364/151‑1 (MQ=255)                                                                                                                                                gcggGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGt   <  1:176796/151‑1 (MQ=255)                                                                                                                                                 cggGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAg    >  2:729661/1‑149 (MQ=255)                                                                                                                                                 cggGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTa  <  2:473170/151‑1 (MQ=255)                                                                                                                                                  gggTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTa  >  1:888888/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CAGACGCATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGAC‑GTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCG‑CCGGTTACCACCAGCAGGCCGAT‑ACCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTA  >  NC_000913/2570888‑2571182

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