Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	47.5%	85.9 / 29.6	61	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (30/2);  new base C (0/29);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.27e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGC  >  NC_000913/1376357‑1376646                                                                                                                                                          |                                                                                                                                              cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGCGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CAcc                                                                                                                                             >  1:704216/1‑151 (MQ=255)  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CAccc                                                                                                                                            >  2:867581/1‑151 (MQ=255)   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCt                                                                                                                                           >  2:926269/1‑151 (MQ=255)     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCg                                                                                                                                         <  1:352739/151‑1 (MQ=255)      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGc                                                                                                                                        >  1:812795/1‑151 (MQ=255)        ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCg                                                                                                                                      >  2:1013617/1‑151 (MQ=255)              aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGccc                                                                                                                                <  2:717891/151‑1 (MQ=255)                              tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                 >  2:495289/1‑150 (MQ=255)                                  cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TGTC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                            >  1:809134/1‑151 (MQ=255)                                   ggcccccGAAAAAGGGGGGGTGTTCCCCGGCGCTTACTGGAAAGGGTCGCCGCCCGAACCGGTCCCGTGGTTTACCCCGACCCGTT‑AAA‑CCGC‑TGTC‑AGGATCATAGGACCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                           <  2:1146627/146‑1 (MQ=255)                                   ggccaccaaaaaaTGGGGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCGTA‑AAC‑CCCC‑TATC‑AGGATCATATGACCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                           <  2:1098705/142‑1 (MQ=255)                                   ggCGCCCGAAAAAGGGTGAGTTATCCCCTGCGCGTACTGGAACAGGTCTCCGCGCGTACCGGACCGGTGGTTTATCTCAACCCGTT‑AAA‑CCGC‑TATC‑GGGATAATATGAACAA‑CCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAAc                                                                                                            <  2:167426/150‑1 (MQ=255)                                   ggCGACCAAAAAATGGTGAGTTTTCCCCTGCGCTTATTGAAACGGTCCCCCGCGGGCCCCGGACCGGGGGTTTATCTCAACCCGTA‑AAC‑CCTT‑TAAC‑AGGATCATATAACCAA‑CCCCCCCCCGGATGCCCGCCGTTACAtcgtcg                                                                                                                 <  1:281970/145‑1 (MQ=255)                                      gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                        >  2:136315/1‑151 (MQ=255)                                        ccGAAAAATGGTGAGTGTTTCCCTGCGCGTCCTGAAACGGGCCGCCGCGCTTACCGGACCCGTGGTGTATCCCGACCCGTT‑AAC‑CCGC‑TATC‑AGGATCATAGGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                       <  2:338581/150‑1 (MQ=255)                                            aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                  >  2:105595/1‑151 (MQ=255)                                               aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                >  2:319173/1‑150 (MQ=255)                                                 ggggagtgaTTCCCGCCGCGTACGGGAACAGTTCGCCGGGCGTCCGGGACGCGGGGTGTATCGCGACCGTTT‑AAA‑CGGC‑TACA‑GGGATCATATG‑TCAACCCCCCTCGCCAAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                              <  2:769903/147‑1 (MQ=255)                                                 ggTGATTGATCCCCTGCGCTTTCGGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGACCCGTA‑AAC‑CCCC‑TATC‑AGGATCATATGATCAC‑CCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                              <  2:370013/150‑1 (MQ=255)                                                  gttagtgaTCCCCTGCGCTTACGGAACCGGGTCGCCGCGCGTACCGAAACGGTGGTGTATCTCGACCCGTT‑AAA‑CCTC‑TACC‑AGAATCATATGACCAA‑CCCCCGCGCGGTGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                <  2:1188458/144‑1 (MQ=255)                                                    gagtgaTTCCCTGCGCTTACTGGAACAGTTCGCCGCGCGTACCGGACCCGTGTTTTATCTCGAACCGTT‑AAC‑CCGT‑TATC‑AGGATCATATGACCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                           <  2:657657/150‑1 (MQ=255)                                                    gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                            >  2:7662/1‑149 (MQ=255)                                                      tttttccccGGCGCGTCCGGAAACGGGTCCCCGCGCGTCCCGAACCGGTGGTTTCTC‑CCAACCCTT‑AAACCCGC‑TACA‑AGGATAATATAACCAC‑CCCCCTCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                         <  1:457884/144‑1 (MQ=255)                                                        gttcccgggggcttACTGGAAGAGGCCGCCGCCCGTACCGGAGCGGGGTTGTACCCCGAACCGTT‑AAC‑CCGT‑TATC‑AGGTTATTATAATCAC‑CCCCCGCCCGAACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                      <  1:928160/138‑1 (MQ=255)                                                            cccGGCGCGTCCGGGAACGGGGCCCCGCGCGACGGGGAACGGGGGTGTATCTCG‑ACCCTT‑AAA‑CCCCTTACC‑AGGATCATAAAACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                       <  1:327419/146‑1 (MQ=255)                                                            cccGGCGCGTCCGGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGACCCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAC‑CCCCCTCGTTGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                     <  1:446412/148‑1 (MQ=255)                                                                 cgcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                             >  2:909131/1‑151 (MQ=255)                                                                  gcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                            >  2:335129/1‑151 (MQ=255)                                                                    gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                          >  2:417685/1‑151 (MQ=255)                                                                    gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                          >  1:324638/1‑151 (MQ=255)                                                                    gTACTGGAACAGGTCGACGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                          >  2:40936/1‑151 (MQ=255)                                                                     tACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                          >  2:1170965/1‑150 (MQ=255)                                                                         ggAACAGGTCGCCGCGCGTCCCGAACCCGTGGTTTATCTCGACCCGTT‑AAC‑CCGC‑TATC‑AGGATCATATGATCAC‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                      <  2:17233/150‑1 (MQ=255)                                                                         gaaaCAGGTCCCCGCGCTTCCCGGAACCGGGGTTTATCTC‑AACCCTT‑AAACCCCC‑TATC‑AGAATCATATGACCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                 <  2:500706/137‑1 (MQ=255)                                                                                   gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                           >  2:671142/1‑151 (MQ=255)                                                                                    ccgcGCGTACCGGAACGGGGGTTTACCTCGACCCGTA‑AAC‑CCGT‑TATA‑GGGATAATATAACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                            <  1:1104509/149‑1 (MQ=255)                                                                                    ccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                           >  1:118311/1‑150 (MQ=255)                                                                                      gcgcgTACCGGAACCGTGGTGTTTCTCAACCCGTT‑AAC‑CCGT‑TATC‑AGAATCATATGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                         <  2:848070/150‑1 (MQ=255)                                                                                          gtcccGAAACCGTGGTGTACCCCAACCCGTA‑AAC‑CCGT‑TATC‑AGAATAATATAAACAC‑CCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                       <  1:454712/145‑1 (MQ=255)                                                                                            aCCGGACCCGTGGTGTACCTC‑AAACTTT‑AAA‑CCCG‑CATCAAGGATTATATAACCAC‑CCCCCTCTCCGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                    <  2:254832/149‑1 (MQ=255)                                                                                                 aaCCGGGGTGTATCTCGACCCGTA‑AAC‑CCTC‑TATA‑AGGATAATATGACAAA‑CCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                              <  2:625205/150‑1 (MQ=255)                                                                                                   ccGTGGTTTATCTC‑AACCTTTAAAA‑CCGC‑TATC‑AGAATCATATGACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                            <  1:849863/150‑1 (MQ=255)                                                                                                   ccGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                           >  2:257375/1‑151 (MQ=255)                                                                                                    cGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                          >  2:1154763/1‑151 (MQ=255)                                                                                                     gtgttttttctcGACCCTTA‑AAC‑CCGT‑TATC‑AGAATCATATAACCAA‑CCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                           <  1:130653/141‑1 (MQ=255)                                                                                                     gtgttttATCCCGACCCGTA‑AAC‑CCTC‑TATC‑AGGATAATATAACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                          <  1:816953/144‑1 (MQ=255)                                                                                                     ggggttTATCCCAACCCGTA‑AAC‑CCGT‑TATC‑AGAATCATATGAACAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                          <  1:856517/148‑1 (MQ=255)                                                                                                       ggttTACCTCGACCCTTT‑AAC‑CCGC‑TATC‑AGATTCATATGATCAC‑CCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                          <  1:989179/148‑1 (MQ=255)                                                                                                          tgtatcccaaCCCTT‑AAACCCGT‑TATA‑AGGAAAATATTACCAC‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                        <  1:453211/140‑1 (MQ=255)                                                                                                                    ctttt‑AAC‑CCGC‑TATC‑AGGATAATATGACCAA‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                               <  1:394539/143‑1 (MQ=255)                                                                                                                        t‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                      >  2:159426/1‑151 (MQ=255)                                                                                                                         taaacccGT‑TATC‑AGGATAATATGATCAA‑CCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                        <  2:67864/145‑1 (MQ=255)                                                                                                                          aaa‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                     >  1:652018/1‑151 (MQ=255)                                                                                                                            a‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                    >  2:1135267/1‑150 (MQ=255)                                                                                                                                    aTC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt             >  2:122849/1‑151 (MQ=255)                                                                                                                                     tC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            >  1:870544/1‑151 (MQ=255)                                                                                                                                      c‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg           >  1:8336/1‑151 (MQ=255)                                                                                                                                          gATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc        >  1:79021/1‑151 (MQ=255)                                                                                                                                           aTCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc        >  1:291886/1‑150 (MQ=255)                                                                                                                                                atGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc  >  1:1102803/1‑151 (MQ=255)                                                                                                                                                      aa‑CCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt      <  2:245507/141‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                              CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTT‑AAA‑CCGC‑TATC‑AGGATCATATGATCAA‑CACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGC  >  NC_000913/1376357‑1376646

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