Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	48.3%	61.8 / 33.1	59	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (26/4);  new base C (0/28);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.25e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376356‑1376651                                                                                                                                                         |                                                                                                                                                   gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAacac                                                                                                                                                   >  1:953680/1‑151 (MQ=255) gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAacac                                                                                                                                                   >  1:207168/1‑151 (MQ=255)    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCt                                                                                                                                                >  2:62334/1‑151 (MQ=255)     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCt                                                                                                                                                <  2:544734/150‑1 (MQ=255)      ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCg                                                                                                                                              <  2:1065156/151‑1 (MQ=255)        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                            >  1:883371/1‑151 (MQ=255)        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                            >  1:792178/1‑151 (MQ=255)                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                     <  1:314070/150‑1 (MQ=255)                      cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                  <  1:565001/147‑1 (MQ=255)                           aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                         >  2:65214/1‑151 (MQ=255)                                ggACGGCGCCCGAAAAATGGTGAGTGATTCCCGGCGCGTACTGAAACAGGTCGCCGCGCGTACCGAACCG‑GTGGTGTATCTCGACCCTTAAAA‑CCGC‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                       <  2:619702/148‑1 (MQ=255)                                   cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                 >  2:473129/1‑151 (MQ=255)                                    ggCGCCCGCAAAATGGTGAGTGATCCCCTGCGCGTACGGGAACGGGTCGCCGCGGGTCCGGGACCG‑GGGGTTTATCCCGACCCTTAAAC‑CCGC‑TACCAGAATAATATGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                <  2:917612/151‑1 (MQ=255)                                    ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                >  2:514612/1‑151 (MQ=255)                                    ggCGACCGAAAAAGGGTGAGTGATCCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGACCCGTAAAC‑CCGC‑TATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                <  2:766343/151‑1 (MQ=255)                                      cgcccGCAAAATGGGGGGTGATCCCCTGCGCTTATTGGAACAGGCCCCCGCGCGTCCCGGACCC‑GTGGTGTATCTCGAACCGTTAAC‑CCGC‑TATCAGGATCATATGATCACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGa                                                                                                               <  2:440857/147‑1 (MQ=255)                                      cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                              >  1:313041/1‑151 (MQ=255)                                        cccaaaaaggggTGATTGATCCCCGGCGCGTGATGGAACGGTTCGCCGCGCGTCCCGGACCC‑GTGTTTTACCCCGAACCTTTAAC‑CCCC‑TATCAGGATCTTATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                            <  2:99503/141‑1 (MQ=255)                                         ccaaaaaaTGGTGAGTGATCCCCTGCGCTTCCTGGAAAGGGTCGCCGCGCGTACGGGACCG‑GGGTTGTATCCCGACCCGTAAAC‑CCGC‑TATAGGTATCATATAACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                <  1:737018/142‑1 (MQ=255)                                         ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                           >  1:958446/1‑151 (MQ=255)                                         ccGAAAAAGGGTGAGTGTTCCCCTGCGCGTCGTGGAACGGGTCGCCGCGCGT‑CCGGAACCGGGGGTGTATCCCGACCCTTAAAC‑CCGT‑TATAGGTATCATATAACCAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGa                                                                                                               <  1:816239/147‑1 (MQ=255)                                          cGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                           >  1:1175407/1‑150 (MQ=255)                                             aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                       >  1:1050801/1‑151 (MQ=255)                                                aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAT‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                    >  2:260085/1‑151 (MQ=255)                                                      agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                 >  1:328613/1‑148 (MQ=255)                                                             cccGGCGCGTACGGGAACAGGTCGCCGCGCGACCCGGACCG‑GTGGTGTTTCTCGACCCTTAAAC‑CCGC‑TATAAGGATAATATGATCAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                              <  1:473109/144‑1 (MQ=255)                                                                   gcttACGGGAACAGGCCGCCGCGCGTCCCGAACCG‑GTGTTGTATCCCGACCTGTTAAC‑CCCT‑TATCAGAATAATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                 <  1:792147/148‑1 (MQ=255)                                                                    cgTTCTGAAACGGGCCCCGGCGCGTCCGGGAACG‑GTGGTGTTTCTCGAACCGTTAAA‑CCGC‑TATCAGGACCATTTGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                  <  2:873578/149‑1 (MQ=255)                                                                      tccggaaaCAGGTCGCCGCGCGTACGGGACCG‑GTGGTGTATCTC‑AACCCTTAAACCCTC‑TATCAGGATCATATAACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                <  2:909242/143‑1 (MQ=255)                                                                       aCTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                             >  1:458921/1‑151 (MQ=255)                                                                         tGGAACAGCTCGCTCCGCGT‑CCAGAACCGGTGGTTTCTCCCTAGCCTTTAAC‑CCGT‑TAAGATCATTTTATGATCCCCTCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                             <  2:840813/149‑1 (MQ=255)                                                                            aacgggcccccgcGCGTCCCGAAACG‑GGGGTTTATCCCAAACCGTAAAC‑CCGC‑TACAAGAATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                        <  1:318305/142‑1 (MQ=255)                                                                            aaCAGGTCCCCGCGCGTCCCGGACCG‑GTGGTTTATCCCAACCCGTAAAC‑CCCT‑TACCAGTATAATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                        <  1:991925/151‑1 (MQ=255)                                                                             acggcgcgccgcGCGTCCCGGACCC‑GTGGTTTATCTC‑AACCCTTAAACCCGC‑TATAAGAATTATAAGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                             <  1:1093108/139‑1 (MQ=255)                                                                             aCAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                        >  1:202348/1‑150 (MQ=255)                                                                                    gccgcGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                >  1:504626/1‑151 (MQ=255)                                                                                     ccgcgctccccGGAACC‑GTGTTGTACCTC‑AACCCTTAAA‑CCCCTTATAAGTATAATATAATCAACCCCCCCCCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                   <  1:370466/138‑1 (MQ=255)                                                                                     ccgcGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                >  1:98830/1‑150 (MQ=255)                                                                                           gtcccGGACCC‑GTGGTTTATCTCGAACCGTTAAC‑CCGC‑TATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                          <  2:313055/147‑1 (MQ=255)                                                                                           gtcccGGAACG‑GTGGTGTACCTCGAACCGTTAAC‑CCGC‑TATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                             <  1:1135088/144‑1 (MQ=255)                                                                                           gTACCGGACCC‑GTGGTGTATCTCGACCCGTTAAA‑CCGC‑TATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                          <  2:404328/150‑1 (MQ=255)                                                                                                 ggaccgggggtttctctcG‑ACCTTTAAACCCGC‑TACAAGAATAATATAACCACCCCCCCCCCCAGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                     <  1:627086/136‑1 (MQ=255)                                                                                                  aaCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                  >  1:383150/1‑151 (MQ=255)                                                                                                         gttgtATCCCAAACCGTTAAC‑CCGC‑TATCAGGATCATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                <  2:178429/145‑1 (MQ=255)                                                                                                          gtgtATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                           >  1:856090/1‑151 (MQ=255)                                                                                                            gtATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAGCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                         >  1:936980/1‑151 (MQ=255)                                                                                                                  cGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                   >  2:185404/1‑151 (MQ=255)                                                                                                                            aacccGC‑TATCAGGACAATTGGACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                            <  2:198124/146‑1 (MQ=255)                                                                                                                               ccGC‑TATCAGAATCATATAACCACCCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                       <  2:917845/151‑1 (MQ=255)                                                                                                                                cGC‑TATCAGAATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                      <  2:718380/151‑1 (MQ=255)                                                                                                                                     accAGAATAATATGACCAACCCCCCCGCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                    <  2:856119/147‑1 (MQ=255)                                                                                                                                           aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc             >  1:67387/1‑150 (MQ=255)                                                                                                                                            tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt           >  1:360785/1‑151 (MQ=255)                                                                                                                                            tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt           >  1:730680/1‑151 (MQ=255)                                                                                                                                             cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg          >  1:641139/1‑151 (MQ=255)                                                                                                                                              atatAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa            <  1:854734/148‑1 (MQ=255)                                                                                                                                              atatAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt           <  1:1084419/149‑1 (MQ=255)                                                                                                                                              aaataaTCCACCCCCCCGCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg          <  1:1112711/145‑1 (MQ=255)                                                                                                                                                 tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc       >  2:1079597/1‑150 (MQ=255)                                                                                                                                                  acaaaaCCCCCCCGCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat     <  1:936351/147‑1 (MQ=255)                                                                                                                                                     cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt  >  2:92134/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                   GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376356‑1376651

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