Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	40.4%	81.1 / 19.5	57	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (10/24);  new base A (23/0);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 2.86e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508520‑3508809                                                                                                                                                |                                                                                                                                                     tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATt‑acaca                                                                                                                                              >  1:396347/1‑150 (MQ=255)  ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑acacac                                                                                                                                             <  1:772779/150‑1 (MQ=255)   tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACg                                                                                                                                            <  1:42448/150‑1 (MQ=255)    ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATC‑ACACACgg                                                                                                                                           >  1:909877/1‑150 (MQ=255)      aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAAttaaa                                                                                                                                                  >  1:819816/1‑140 (MQ=255)       aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGc                                                                                                                                       <  2:674093/151‑1 (MQ=255)        aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑GGGAAAAATC‑AAAAACAGTAGcc                                                                                                                                      >  2:76287/1‑151 (MQ=255)        aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGAAAAAATC‑ACACACGGGAg                                                                                                                                        >  2:646812/1‑149 (MQ=255)        aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTA‑AAAAACggaacc                                                                                                                                       >  2:92727/1‑147 (MQ=255)         gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGG‑GGGAAAAATC‑ACACAACGTAGacc                                                                                                                                     >  1:567696/1‑151 (MQ=255)              agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCcaca                                                                                                                                <  1:664959/151‑1 (MQ=255)                     aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATC‑AAACACCGGAGCCCCCCCACCCGaa                                                                                                                          >  2:924704/1‑150 (MQ=255)                      atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGCGCGC‑GGGAAAAATC‑ACCCACAGTAGCCCCCCCAt                                                                                                                               >  1:147950/1‑144 (MQ=255)                        atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                       <  1:415142/150‑1 (MQ=255)                         tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATC‑ACACACCGGAGCCCCCCCATCAGaa                                                                                                                          >  1:147772/1‑146 (MQ=255)                            aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGG‑GGGAAAAATT‑ACACACGGTAGCCCCCCCAACCGAAAACgcgg                                                                                                                   >  2:819102/1‑150 (MQ=255)                            aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAGAAATTA‑ACAAACGGGAACCCCCCCCTCGGaagga                                                                                                                       >  2:84910/1‑143 (MQ=255)                              aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCTACAAACGG‑AGACCCCCCAACAGAAAACggcgtg                                                                                                                 >  2:490819/1‑145 (MQ=255)                                      gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                                        <  2:767046/151‑1 (MQ=255)                                       ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                       <  1:65362/151‑1 (MQ=255)                                         cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATC‑ACACACGGTAGCCCCCACATCAGAAAACGTGGGGGGCACAcc                                                                                                         >  1:82253/1‑147 (MQ=255)                                           aggtgggaaCCCGTCGGATCGAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGAAGAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GTAAAAATTC‑ACACACGGTGGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                                      <  2:523647/141‑1 (MQ=255)                                           aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGC‑GGAAAAATTT‑ACACACGATAGGCCCCACAACATAAAGCGTGGTGCGCACAACCAGCTg                                                                                                   >  1:218748/1‑151 (MQ=255)                                           aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                   <  2:353891/151‑1 (MQ=255)                                           aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGAGCGC‑GGAAAAAATT‑CCACACAGGAACCCCCCCAACAGAAAACgcgggggggccacc                                                                                                         >  2:502301/1‑135 (MQ=255)                                            gggTTTAACCCGTCGGATAAAAACTCAGAGGGTGTCCCCGGAAAACCGGAGGAAGGAAAGAAAAGACGCACCATTTTGTGGCGGGGGCGCGC‑GAAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGGGTGGGGCGCACACCGaatt                                                                                                    <  2:159459/149‑4 (MQ=255)                                               tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATC‑ACACACCCAAGGCCCCCCCTCCGAAAGCGTGGGGGGGCCACCGCGCTGCggg                                                                                               >  2:532146/1‑149 (MQ=255)                                               tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATT‑ACACACGGTAGCCCCCCCCTCAGAAAGCGGGGGGGGCACACCGAGCTGCg                                                                                                 >  1:204894/1‑149 (MQ=255)                                                gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGCAAAATTC‑CCACACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACAcc                                                                                                         >  1:160142/1‑140 (MQ=255)                                                  aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATT‑ACACAAGGTCGCCCCCCCATCCGAAAGCGGGGTGCGCACACCCAGCTGCg                                                                                                 >  2:565603/1‑146 (MQ=255)                                                    cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                          <  2:155161/151‑1 (MQ=255)                                                              cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                                <  1:101170/151‑1 (MQ=255)                                                                   cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGAAAAAATC‑ACACACGGGAGCCCCCACAACAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGGTCATCCCCTGCTGAt                                                                            >  1:577924/1‑150 (MQ=255)                                                                      aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGAAAAAATC‑ACAAACGGGAGCCCCCACCTCAGAAAGAGGGGTGGGCACACCGCACTGCGTAGTGGTCATCCCCTGCTGGTaaa                                                                         >  2:881678/1‑150 (MQ=255)                                                                         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATT‑ACACACGGGAGACCCCCCACCAGACAGAGAGGGGGGCGCACCCAGCTGCGGAGCGAGTATCCCTCGCTGTTaaaac                                                                       >  2:541741/1‑147 (MQ=255)                                                                           gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGCGG‑GGGAAAATTT‑ACACAACGGAGGCCCCCCACCAAGAAACGAGGGGGGCACACCGCGCCGCGGCGGGGTCAACCCCTTCTGGTAAAACAGc                                                                    >  1:384128/1‑150 (MQ=255)                                                                             gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                                 <  2:336918/151‑1 (MQ=255)                                                                                   tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGAGCGG‑GGGAAAAATA‑AAAAACCGGAAACCCCCCACCAGAAAGCGGGGGGGGCACACCGCGCTGCGGAGCGATCATCCCCTGCTGATAAAACAACACCGGcacc                                                           >  2:266205/1‑151 (MQ=255)                                                                                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                                          <  1:464029/151‑1 (MQ=255)                                                                                         ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg                                                     <  2:398387/151‑1 (MQ=255)                                                                                           ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGC‑GGGAAAATTC‑ACACCCGGTAGACCCCCCATCAGAACGCGGGGGGCGGACACCGCGCTGCGGAGCGATCCTCCCCTTCTGATAAAACAGCACAGGCACCACCTcgcg                                                   >  2:212849/1‑151 (MQ=255)                                                                                           ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATC‑ACACACGGTAGCCCCCCCATCAGAAAACGTGGGGGGCACACCGAGCCGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACTTcgcg                                                   >  1:268388/1‑151 (MQ=255)                                                                                             agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAATTT‑ACACACCGGAGCCCCCCCCTCAGAAAGCGGGGGGCGCCCAACCGGCCGCGTAGCGATCATCCCCTGCTGATAAACCAGCACTGGCACCACCTCCCggg                                                 >  2:875681/1‑151 (MQ=255)                                                                                                 aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtg                                             <  1:21741/151‑1 (MQ=255)                                                                                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATC‑ACACACGGTAGGCCCCCCATCAGAAAACGGGGGGCGCACCCCGCGCCGCGTAGCGATCATCCCTTTCTGATAAAACAACACCGGCCCCCCCTCTCGGGGAtgggg                                          >  2:937443/1‑149 (MQ=255)                                                                                                          aagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATT‑AAACACCGTAACCCCCACATCAGAAAGCGGGGGGCGCACACCGAGCTGCGTGGCGATCATCCCTTTCTGTTAAACCAGCACTGGCACCCCTTTGCGGGGATggtgggtg                                      >  1:247705/1‑149 (MQ=255)                                                                                                                 aCCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                             <  1:153239/151‑1 (MQ=255)                                                                                                                  ccATTTTGTGGCGGTGGCGCGC‑CGGAAAAATC‑ACAAACCGGAGCCCCCCCATCAGAAAGCGTGGTGGGCACACCGAGCCGCGGCGAGGTCATCCCCTGCTGATAAAACAACACTGGCACCACCTCGCGGGGATGGGGGGGGTGACCAAAgg                            >  2:882837/1‑151 (MQ=255)                                                                                                                         gtggcggtggcgCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                      <  1:490091/150‑1 (MQ=255)                                                                                                                          tggcggtggcgCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                    >  1:878800/1‑151 (MQ=255)                                                                                                                           ggcggtggcgCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                    <  2:540744/150‑1 (MQ=255)                                                                                                                              ggtggcgCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc                >  1:284179/1‑151 (MQ=255)                                                                                                                              ggtggcgCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc                >  2:355865/1‑151 (MQ=255)                                                                                                                               gtggcgCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc               <  2:564196/151‑1 (MQ=255)                                                                                                                                   cgcgc‑gGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGa           <  1:239758/151‑1 (MQ=255)                                                                                                                                      gc‑gGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc        <  1:496878/151‑1 (MQ=255)                                                                                                                                      gc‑gGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc        <  1:608495/151‑1 (MQ=255)                                                                                                                                             aaaTTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCTTTGCCGTGATCCGCCATc  <  2:789979/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     TTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTC‑ACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508520‑3508809

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