Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,601,784	0	A	G	39.6%	76.6 / 16.0	54	E195G (GAG→GGG)	hyfA	hydrogenase 4, 4Fe‑4S subunit
Reads supporting (aligned to +/- strand):  ref base A (6/26);  new base G (21/0);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATC  >  NC_000913/2601639‑2601933                                                                                                                                                   |                                                                                                                                                       gACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑Ggggggg                                                                                                                                                      >  2:478839/1‑148 (MQ=255)     tgtgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGATcccc                                                                                                                                                >  1:769734/1‑150 (MQ=255)      gtgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCt                                                                                                                                               <  1:123224/150‑1 (MQ=255)       tgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGGTCCCttt                                                                                                                                             >  1:335082/1‑151 (MQ=255)          gATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTccc                                                                                                                                          <  2:228268/151‑1 (MQ=255)             ttCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGGCCCCCTTccccc                                                                                                                                        >  1:613208/1‑148 (MQ=255)              tCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGGTCACCCTTcccctac                                                                                                                                      >  2:857681/1‑147 (MQ=255)              tCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCctca                                                                                                                                       <  2:278612/150‑1 (MQ=255)                 tGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCctcact                                                                                                                                     <  1:477534/149‑1 (MQ=255)                  gCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCctcactca                                                                                                                                   <  2:12138/150‑1 (MQ=255)                    cAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGa                                                                                                                                <  2:319163/151‑1 (MQ=255)                        aGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCTTTTCCTCCCTCAGGAg                                                                                                                               >  1:840314/1‑148 (MQ=255)                         ggTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCCTTTCctccccc                                                                                                                                    >  2:297689/1‑138 (MQ=255)                          gTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGGACCCCCTTCCCTCACCCAGGGGCaa                                                                                                                            >  1:306493/1‑149 (MQ=255)                           tCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGGGGGCGCCCCTTTCCttcaccacggg                                                                                                                                >  1:91220/1‑139 (MQ=255)                             cGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTa                                                                                                                      <  1:603075/151‑1 (MQ=255)                                ccTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATg                                                                                                                   <  1:27370/151‑1 (MQ=255)                                    cGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGAACCCCCTTTCCCCCCTCAGGGGCAAAG‑CTAAAGGa                                                                                                                 >  1:311832/1‑149 (MQ=255)                                        cgcgCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCTTTCCCCCACCCCGGAGCAACG‑CTAATGGATGccc                                                                                                            >  1:626481/1‑150 (MQ=255)                                           gcTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCa                                                                                                        <  1:834369/151‑1 (MQ=255)                                               gCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAatta                                                                                                    <  1:738066/151‑1 (MQ=255)                                                cccGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGATCCCCCTTTCCCCAATCAAGAACAACGCCTAATGGAGGCCCtcaaat                                                                                                      >  2:741310/1‑144 (MQ=255)                                                   gAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTaa                                                                                                <  2:806440/151‑1 (MQ=255)                                                         ggCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGt                                                                                          <  2:138543/151‑1 (MQ=255)                                                            gTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGGCCCCTTTCCCTCACTCACGAGGAACG‑CTCATGGACTCCCTGCCattttt                                                                                                  >  2:346393/1‑140 (MQ=255)                                                              tACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCt                                                                                      <  2:554000/150‑1 (MQ=255)                                                               aCGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCCTTCCCTCCCCTAGGGGCCACC‑CTAATAGATGGCCCGGAATTATTTACCt                                                                                             >  1:164325/1‑142 (MQ=255)                                                               aCGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGa                                                                                    <  1:576326/151‑1 (MQ=255)                                                                cGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGAACCCCTTTCCCCCACTCCGGAGCAAAG‑CTAATGGAAGCCCTGCCATTATTAACCTGGTCGCTg                                                                                     >  1:271208/1‑149 (MQ=255)                                                                 ggCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGAtt                                                                                  <  2:131852/151‑1 (MQ=255)                                                                    tGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCTTTTCCTCACTCAGGGGCAAAG‑CTAATGGGTGGCCCGCAATTATTAACCTGGTCGCTGAtaga                                                                                >  1:29992/1‑147 (MQ=255)                                                                    tGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGAt                                                                                   <  2:29988/147‑1 (MQ=255)                                                                      aTCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTa                                                                             <  2:259695/151‑1 (MQ=255)                                                                          ccGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTg                                                                         <  2:409964/151‑1 (MQ=255)                                                                           cGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGACCCCCTTCCCCCCCCCAGGAGGAACC‑CCAAAGGGTGGCCCGCCATTATTTAACTGGTCGCTGATTTTCTATc                                                                           >  1:4058/1‑148 (MQ=255)                                                                           cGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATACCCTTTCCCCCCCTCAGGGGCCACG‑CTCATGGATGGCCTGCAATTATTAACATCGGCGCTGAGTTTCTATCTGt                                                                        >  2:97815/1‑151 (MQ=255)                                                                                aTAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGCGGGACCCCCTTTCCTCACCCACGAGCAACC‑CTAAAGGATGGCCCCCAATTATTCACATCGTCGCTGAATTTCTATTTCTTTGCt                                                                   >  2:469013/1‑151 (MQ=255)                                                                                  aGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAg                                                                 <  2:541032/151‑1 (MQ=255)                                                                                     cTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGGGGGACCCCCTTTCCC‑CACTCAGGGGCAACC‑CTAATGGGGGCCCCGCAATTATTTACCTGGGCGCTGGTTCTCTCTCTGTTTGCtcattt                                                              >  2:691330/1‑146 (MQ=255)                                                                                       gCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGAACCCCCTTCCCCCACTCAGGAGCACCA‑CGAATGGATGCCCTGCAATTATTTACCTGGGCGCTGATTCTCTATTTGTTTTCTTGTCt                                                              >  2:243194/1‑149 (MQ=255)                                                                                            gTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCg                                                       <  2:220740/151‑1 (MQ=255)                                                                                              cAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCCTTTCCCCACTCAAGAGCAACG‑CTAATGGATGCCCCGCAATTATTAAACTGGTCGCTGGTTTTCTTTCTGTTTTCTAGTCTGGCtt                                                         >  1:59669/1‑147 (MQ=255)                                                                                                   gaaagaaaAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGAACCCCCTTCCCTCACTCCGGAGCAACG‑CTAATGGGTGCCCTGCAATTATTTAACCGGGCGCCGGTTTTCTATTTGTTTTCTAGTCTGGGTTCGCTGtttt                                                >  1:852135/1‑151 (MQ=255)                                                                                                         aaaaCAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCg                                          <  1:659554/151‑1 (MQ=255)                                                                                                            aCAGCGCCTGGCCGCCAGCTGGCTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCCCTCAGGAGCACCG‑CTAATGGCTGCCCCGCTATTATTAACCTGTTCGCCGAGTCTCTATCTGTTTGCTAGTTTGGCTTCGCTGTTTTTACTCGGt                                        <  2:122714/150‑1 (MQ=255)                                                                                                                cgcCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGTTCCCCCTTCCCCCCCTCCCGAGCAACA‑CCACAGGGTGCCCTGCAATTTTTTCCACTGTCGCCGGTTCTCTCTCTGTTTGCTTGTGTGGGTTCGCTGctttttttttttctggg                                   >  2:102951/1‑139 (MQ=255)                                                                                                                  ccTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACa                                 <  2:840053/151‑1 (MQ=255)                                                                                                                     tGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCCTTTCCTCACCCCGGGGCAAAG‑CCAATGGATGCCCTCCCCTTTTTTATCTTGTCGCTGATTTTCTATCTTTTTTCTTGTCTGGGTTCGCTGTTTTTAtttggtttgg                                    >  2:630973/1‑145 (MQ=255)                                                                                                                     tGCCGCAAGCTGGTTTGCCAATGGC‑GGGGGGGAACCCCTTTTCCTCACTCCGGCGGAACC‑CTAATGGGTGGCCTGCAATTATTTACCTGGGCGCTGCTTCTTTTTTTGGTTGTTAGTTTTGCTTCTCTGTTTTTTCTTGGTCTGGTCAg                                >  2:332587/1‑149 (MQ=255)                                                                                                                      gCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGAc                              <  2:134824/150‑1 (MQ=255)                                                                                                                         gCAAGCTGGTTTGCCAATGGC‑GGGGGGGAACCCCTTTTCCTCCCTCAGGAGCAAAG‑CCAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTTTCTTTTTTTTTGCTAGTCTGGCTTCGCTGTTTTTTCTCGGGCTgggaagaacgg                           >  1:828316/1‑141 (MQ=255)                                                                                                                           aaGCTGGTTTGCCAATGGC‑GGGGGGGGTCCCCTTTTCCTCCCTCAGGGGCAACG‑CTAATGGGTGCCCCGCAATTATTAACCTGGTCGCTGATTTTCTATCTGTTTGCTAGTCTGGCTTTGCTGTTTTTTCTCGGTCTGGGCAGACTgg                           >  1:184203/1‑148 (MQ=255)                                                                                                                           aaGCTGGTTTGCCAATGGC‑GGGGGGGGGCCCCTTTTCCTCAATCAAGAACAACG‑CCAATGGGTGCCCCTCGAATTTTAAACTGGTGGCGGTATTTCTCTCTGTTTGTTAGTCCGGCTTCGCTGTTTTTATCCGGTCTGGGCAGACCGGCt                         >  2:398546/1‑150 (MQ=255)                                                                                                                                   ttGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCt                 <  1:587782/150‑1 (MQ=255)                                                                                                                                       cAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCg            <  2:858205/151‑1 (MQ=255)                                                                                                                                                 ggaggaTCCCCTTTCCCTCACTCAGGAGCACCG‑GTAATGGATGCCCTGCAAGTATTAGCCTGCTCGGTGATTCTCTATCTGTTTGGTAGTCTGGCATCGCTGTTTTTACTCGGTCTGGACAGAGTTGCTGTTAAGCCTTCCGGCATcaca    <  2:409942/150‑1 (MQ=255)                                                                                                                                                  gaggaTCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATc  <  1:844157/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGC‑GGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACG‑CTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATC  >  NC_000913/2601639‑2601933

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