Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,987,760	A→C	5.9%	E93A (GAG→GCG)	oplA →	5‑oxoprolinase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,987,760	0	A	C	5.9%	211.4 / 2.5	67	E93A (GAG→GCG)	oplA	5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base C (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGCCGGCCTGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTGGCGATGGCCGTGAGCA  >  NC_002947/3987613‑3987907                                                                                                                                                     |                                                                                                                                                    gcgcCGGCCTGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAg                                                                                                                                                    <  6:160372/149‑1 (MQ=255) gcgcCGGCCTGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAg                                                                                                                                                    >  2:219982/1‑149 (MQ=255)  cgcCGGCCTGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCACCCGGGGTCAAGGTGGGCCTGCTGTGCACCGAg                                                                                                                                                    <  1:219982/149‑1 (MQ=255)  cgcCGGCCTGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAgg                                                                                                                                                   <  7:191853/149‑1 (MQ=255)          tGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAg                                                                                                                                           >  6:14763/1‑149 (MQ=255)              acaGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACa                                                                                                                                       <  4:143177/149‑1 (MQ=255)                aGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCCCCTCGGTGCAcagc                                                                                                                                           <  6:173953/143‑4 (MQ=255)                aGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCCCCTCGGTGCAcagc                                                                                                                                           >  5:173953/1‑140 (MQ=255)                aGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGc                                                                                                                                     <  3:352983/149‑1 (MQ=255)                      ccGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCCCCTCGGTGCAcagcag                                                                                                                                         <  4:146548/139‑6 (MQ=255)                      ccGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCCCCTCGGTGCAcagcag                                                                                                                                         >  3:146548/1‑134 (MQ=255)                        gACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAAt                                                                                                                             >  7:24265/1‑149 (MQ=255)                           gCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCg                                                                                                                          >  3:309755/1‑149 (MQ=255)                                 gggCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTggg                                                                                                                    <  2:117428/149‑1 (MQ=255)                                  ggCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTgggg                                                                                                                   <  5:14763/149‑1 (MQ=255)                                   gCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGc                                                                                                                  >  8:83104/1‑149 (MQ=255)                                          ccGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAgg                                                                                                            >  2:249490/1‑148 (MQ=255)                                             gCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGgaaga                                                                                                        <  8:24265/149‑1 (MQ=255)                                              ccGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGgaagaa                                                                                                       <  5:232967/149‑1 (MQ=255)                                                gAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGa                                                                                                                                            <  1:195460/110‑1 (MQ=255)                                                gAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGa                                                                                                                                            >  2:195460/1‑110 (MQ=255)                                                          ccAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTac                                                                                           >  4:78189/1‑149 (MQ=255)                                                          ccAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTac                                                                                           >  6:273879/1‑149 (MQ=255)                                                                 tgACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGcc                                                                                     <  4:309755/148‑1 (MQ=255)                                                                 tgACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGcc                                                                                     >  3:269150/1‑148 (MQ=255)                                                                    ccTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCAtt                                                                                 <  8:291014/149‑1 (MQ=255)                                                                            tCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCt                                                                                                                       >  8:215425/1‑103 (MQ=255)                                                                            tCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCt                                                                                                                       <  7:215425/103‑1 (MQ=255)                                                                            tCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACccgccg                                                                         >  7:269718/1‑149 (MQ=255)                                                                               aCGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGca                                                                      >  7:317829/1‑149 (MQ=255)                                                                                     ccaccGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCAc                                                                                                                                              <  8:81461/71‑1 (MQ=255)                                                                                     ccaccGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCAc                                                                                                                                              >  7:81461/1‑71 (MQ=255)                                                                                      caccGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTgg                                                                                                                                 <  3:200746/83‑1 (MQ=255)                                                                                      caccGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTgg                                                                                                                                 >  4:200746/1‑83 (MQ=255)                                                                                       accGTGGCACTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGAACAAGGAAGAAGCACACCGCTACGACGCCCATTACCCGCCGCCACACATGCt                                                              >  3:270283/1‑149 (MQ=255)                                                                                        ccGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCATCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCtt                                                             >  4:277509/1‑149 (MQ=255)                                                                                           tGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGcc                                                          <  5:273879/149‑1 (MQ=255)                                                                                               ccTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCAc                                                      >  7:29883/1‑149 (MQ=255)                                                                                                   aaCGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAg                                                                                                                                           <  3:39665/60‑1 (MQ=255)                                                                                                   aaCGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAg                                                                                                                                           >  4:39665/1‑60 (MQ=255)                                                                                                            aTCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTggg                                                                                                                    <  5:176276/74‑1 (MQ=255)                                                                                                            aTCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTggg                                                                                                                    >  6:176276/1‑74 (MQ=255)                                                                                                               gAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGccc                                      >  1:213708/1‑149 (MQ=255)                                                                                                                   gcA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCa                                                       >  4:86359/1‑128 (MQ=255)                                                                                                                   gcA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCa                                                       <  3:86359/128‑1 (MQ=255)                                                                                                                    cA‑CCGGGGTCAAGGTGGGCCTGCTGTGCATCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCg                                 <  3:277509/149‑1 (MQ=255)                                                                                                                        cGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCtt                                                             <  4:207597/118‑1 (MQ=255)                                                                                                                        cGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCtt                                                             >  3:207597/1‑118 (MQ=255)                                                                                                                        cGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTg                              >  8:153533/1‑149 (MQ=255)                                                                                                                            gTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGcc                           >  4:234798/1‑148 (MQ=255)                                                                                                                              cAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGaa                        <  2:146764/149‑1 (MQ=255)                                                                                                                               aaGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGc                            >  1:29301/1‑144 (MQ=255)                                                                                                                               aaGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGc                            <  2:29301/144‑1 (MQ=255)                                                                                                                                aGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATc                      >  2:126936/1‑149 (MQ=255)                                                                                                                                aGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATc                      <  4:269150/149‑1 (MQ=255)                                                                                                                                  gTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCt                    >  3:154279/1‑149 (MQ=255)                                                                                                                                     ggCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTgg                 <  3:78189/149‑1 (MQ=255)                                                                                                                                     ggCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTgg                 <  7:83104/149‑1 (MQ=255)                                                                                                                                     ggCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTgg                 <  7:153533/149‑1 (MQ=255)                                                                                                                                     ggCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTg                  <  8:29883/148‑1 (MQ=255)                                                                                                                                     ggCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTg                  <  2:213708/148‑1 (MQ=255)                                                                                                                                        ctgctgTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTGGCGa              <  6:86432/149‑1 (MQ=255)                                                                                                                                          gctgTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCtt                                                             >  2:176509/1‑100 (MQ=255)                                                                                                                                          gctgTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCtt                                                             <  1:176509/100‑1 (MQ=255)                                                                                                                                           ctgTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTGGCGAt             <  8:317829/147‑1 (MQ=255)                                                                                                                                                 aCCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTGGCGATGGCCGTGa     <  3:234798/149‑1 (MQ=255)                                                                                                                                                    gAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTGGCGATGGCCGTGAGca  <  1:249490/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCGCCGGCCTGGCACAGATCGCCGACGCCACCGGGCGTAGCCCGGCCGAACTGATCGCCAACTGTGACCTGTGCATCAACGGCACCACCGTGGCCCTCAACGCGCTGATCGAGCGCA‑CCGGGGTCAAGGTGGGCCTGCTGTGCACCGAGGGGCACGAGGACAGCCTGGAAATCCGCCTGGGGCACAAGGAAGAAGGACACCGCTACGACGCCCATTACCCGCCGGCACACATGCTTGCCCCACGCCACCTGCGCCGCCCGATCGGTGGCCGAATCCTTGGCGATGGCCGTGAGCA  >  NC_002947/3987613‑3987907

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