Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,615	0	T	G	48.2%	83.1 / 13.2	56	T15P (ACC→CCC)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base T (7/22);  new base G (27/0);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 7.66e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGG  >  NC_000913/4137465‑4137764                                                                                                                                                       |                                                                                                                                                      cGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGt                                                                                                                                                       <  2:142257/151‑1 (MQ=255)     aGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGGTTTTAAGGAGTGGCTGGGaaa                                                                                                                                                    >  1:640191/1‑150 (MQ=255)            gtaagCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTa                                                                                                                                            <  1:126837/151‑1 (MQ=255)             taagCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGaaaatgcggt                                                                                                                                             >  2:18988/1‑142 (MQ=255)                  tGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGATTTTTAGGAGTGGGTGGGAAAAGGCGGTAGCAgg                                                                                                                                       >  1:118258/1‑150 (MQ=255)                         aCGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGGTTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGGTGag                                                                                                                                  >  2:767190/1‑148 (MQ=255)                         aCGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGagca                                                                                                                                <  2:664904/150‑1 (MQ=255)                          cGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGagcagc                                                                                                                              <  2:540953/151‑1 (MQ=255)                          cGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGCTGag                                                                                                                                  >  2:215399/1‑147 (MQ=255)                           gCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAgcg                                                                                                                             <  1:26980/151‑1 (MQ=255)                           gCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAgcg                                                                                                                             <  1:549078/151‑1 (MQ=255)                               gATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGGGGCTGGGGAAAGCCGTTAGCAGGGTTGGCAgcggc                                                                                                                           >  1:479265/1‑147 (MQ=255)                               gATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAg                                                                                                                         <  1:846925/151‑1 (MQ=255)                                aTCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGt                                                                                                                        <  2:258284/151‑1 (MQ=255)                                 tCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGGTAGCAGGGTGAGGAGCGcc                                                                                                                           >  1:155122/1‑147 (MQ=255)                                    tAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGGTTTTTAGGAGTGGGTGGGAAATGGCGTTAGCGGGGTGAGCCGCGCCCGt                                                                                                                        >  2:99890/1‑147 (MQ=255)                                     aaaaTAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAg                                                                                                                    <  1:161080/150‑1 (MQ=255)                                                ctgcCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGCTGAGCAGCGCCCGTGAAGGTTTCATCTgg                                                                                                         >  2:657694/1‑150 (MQ=255)                                                    cAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAAAGCCGTTAGCAGGCGGGGGAGCGCCAGGGAAGCTTTTATCTgggt                                                                                                       >  1:861111/1‑146 (MQ=255)                                                      aTTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAAAGCCGTTTGCAGGGTGAGCAGCGCCAGTGAAGCTTTTATCTGGGTTGTcc                                                                                                  >  2:446037/1‑151 (MQ=255)                                                      aTTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTcc                                                                                                  <  1:651191/151‑1 (MQ=255)                                                       tttGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAAAGCCGTTAGCAGGGTGAGCAGCGCCAGGGAAGCTTTCATCTGGtttggcc                                                                                                  >  2:612685/1‑145 (MQ=255)                                                        ttGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAAAGCCGGTTGGAGGGTGGGGAGCGCCCGGGAAGCttttcttttggt                                                                                                       >  1:644300/1‑135 (MQ=255)                                                        ttGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGCTGAGCAGCGCCCGTGAAGCTTTTATCTTGTTCGTcc                                                                                                  >  2:251120/1‑149 (MQ=255)                                                         tGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGGTGGGAAAAGGCGTTAGCAGGGTGAGGAGGGCCCAGGGAGGTTTTATCTCGTTCGGCCCTg                                                                                               >  1:435802/1‑151 (MQ=255)                                                            tGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAAAGCCGGTTGCAGGGTGAGGAGCGCCAGTGAAGCGTTTATCttggt                                                                                                       >  1:144799/1‑136 (MQ=255)                                                            tGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTg                                                                                               <  1:29780/148‑1 (MQ=255)                                                            tGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAAAGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTTGTTCGTcc                                                                                                  >  2:29780/1‑145 (MQ=255)                                                             gTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGAATTTTAGGCGGGGCTGGGAAAATCCGCTAGGAGGGCGGGCAGGCGCAGAGGAGCTTTTATCTGGtttctctcgggaa                                                                                            >  1:555175/1‑139 (MQ=255)                                                                tGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAAAGCCGGTAGCAGGCTGAGCAGCGCCCGGGAAGCTTTCATCTGGTTCGTCCGTGgaaaaa                                                                                         >  2:23118/1‑148 (MQ=255)                                                                         cTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGATTTTAAGGAGTGGCTGGGAAATGCCGGTTGCAGGCTGCGGAGCGCCCCTGAAGGTTTCCTCTGGTTTGTCCGGGGAAATAAGAGGcaca                                                                                >  1:636081/1‑150 (MQ=255)                                                                           gCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGGTGAGCAGCGCCCGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAAAGTca                                                                                  >  2:169648/1‑146 (MQ=255)                                                                             gATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAAAGCCGTTAGCAGGGGGAGCAGCGGCAGTGAAGCTTTCATCTGGGTTGTCCGGGAAAAAAAGAGGCACAGAtt                                                                            >  2:34115/1‑150 (MQ=255)                                                                                             ttGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGGTAGCCGGGTGAGCAGCGCCAGTGAAGCTTTTATCTGGTTCGTCCGGGAAAATAAAAGGCACAGATTCGACCTTCCGGGCTaa                                                            >  2:209083/1‑150 (MQ=255)                                                                                                gCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGAt                                                        <  1:374074/151‑1 (MQ=255)                                                                                                    ggATTTGTACCACGGTTGGTGGAACGGCGGGGGATTTTTAGGGGTGGCTGGGAAATACCGTTTTCAGGCTGAGCAGCGCCAGTGAAGATTTTATCTGGTTCGTCCGGGGAAATAAAAGGCACAGATTTTACCTTCCGGGCTAACGATTGtt                                                    >  2:499479/1‑151 (MQ=255)                                                                                                      aTTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAAAGCCGTTTGCAGGGGGAGCAGCGCCCGTGAAGGTTTTATTTGGGTCGTCCGTGGAAATAAGAGGCACAGAATCGCCCCTCCCGGGTAACGAtt                                                       >  1:331352/1‑146 (MQ=255)                                                                                                        ttGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTTATCTGGGTCGTCCGTGAAAATAAGAGGCACAGATTCGACCTTCCGGGCTAACGATTGTTTTGt                                                >  2:726942/1‑151 (MQ=255)                                                                                                               aCGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAAAGCCGGTAGCAAGGTGGGGCGCGCCAGGGGAGCTTTTATCTGGTTCGTCCGGGGAAAAAAGAGTCCCAGAATCGACCCTCCGGGCTAACGATTGtttt                                                  >  1:639001/1‑140 (MQ=255)                                                                                                                  gTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATg                                      <  2:75778/151‑1 (MQ=255)                                                                                                                        ggAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGa                                <  2:861209/151‑1 (MQ=255)                                                                                                                         gAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGaa                               <  2:388835/151‑1 (MQ=255)                                                                                                                          aaCGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAg                              <  1:344775/151‑1 (MQ=255)                                                                                                                                gggAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGAt                         >  1:480030/1‑150 (MQ=255)                                                                                                                                 ggAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATa                        >  1:619442/1‑150 (MQ=255)                                                                                                                                 ggAGATTTTAAGGAGTGGCTGGTAAATGCCGTCAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATa                        >  1:433980/1‑150 (MQ=255)                                                                                                                                  gagaTTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATaaa                      <  2:751324/151‑1 (MQ=255)                                                                                                                                  gagaTTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATaa                       >  1:105200/1‑150 (MQ=255)                                                                                                                                     aTTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATaaaaaa                   <  2:118266/151‑1 (MQ=255)                                                                                                                                          aaGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCt              <  2:815952/151‑1 (MQ=255)                                                                                                                                          aaGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCt              <  2:98666/151‑1 (MQ=255)                                                                                                                                            ggAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTg             >  2:424059/1‑150 (MQ=255)                                                                                                                                             gAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGaa           <  1:660761/151‑1 (MQ=255)                                                                                                                                                   cTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGt     >  1:548371/1‑151 (MQ=255)                                                                                                                                                     ggTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTg    <  1:451284/150‑1 (MQ=255)                                                                                                                                                      gTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTggg  <  2:161364/151‑1 (MQ=255)                                                                                                                                                      gTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTggg  >  1:690586/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGTGAAGCTTTCATCTGGTTCGTCCGTGAAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGG  >  NC_000913/4137465‑4137764

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