Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,116,512	C→A	6.3%	P33Q (CCG→CAG)	PP_4502 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,116,512	0	C	A	6.3%	199.5 / 3.8	63	P33Q (CCG→CAG)	PP_4502	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base A (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

CCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATCGC  >  NC_002947/5116364‑5116649                                                                                                                                                     |                                                                                                                                          ccGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGGCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCgcccgccc                                                                                                                                           <  6:84594/149‑1 (MQ=255)  cGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGcccgcccg                                                                                                                                          >  3:324953/1‑149 (MQ=255)    tGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGgc                                                                                                                                        >  5:30438/1‑149 (MQ=255)        cTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGgcgc                                                                                                                                      <  4:15592/147‑1 (MQ=255)        cTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGgcgc                                                                                                                                      >  3:15592/1‑147 (MQ=255)            cTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGcaca                                                                                                                                   >  2:93943/1‑146 (MQ=255)            cTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGACAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGcaca                                                                                                                                   <  1:93943/146‑1 (MQ=255)               gtggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAAc                                                                                                                             >  5:212153/1‑149 (MQ=255)               gtggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAAc                                                                                                                             >  4:299273/1‑149 (MQ=255)                 ggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGc                                                                                                                           >  8:59467/1‑149 (MQ=255)                 ggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGc                                                                                                                           <  6:30438/149‑1 (MQ=255)                    gcAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTg                                                                                                                        >  7:157510/1‑149 (MQ=255)                                ccTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaag                                                                                                            <  8:157510/149‑1 (MQ=255)                                   gTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggcc                                                                                                         >  8:169898/1‑149 (MQ=255)                                    tCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggccg                                                                                                        <  5:226345/149‑1 (MQ=255)                                    tCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggcc                                                                                                         >  3:276366/1‑148 (MQ=255)                                      gcgcCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggccg                                                                                                        >  4:253061/1‑147 (MQ=255)                                             cTGCACGATGACCCTTCTACCTCCAAACACAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCcgcg                                                                                               >  1:16444/1‑149 (MQ=255)                                              tGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgt                                                                                              >  1:246582/1‑149 (MQ=255)                                              tGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgt                                                                                              >  2:45243/1‑149 (MQ=255)                                               gCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgtg                                                                                             >  5:277028/1‑149 (MQ=255)                                                cACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgtgg                                                                                            >  3:256103/1‑149 (MQ=255)                                                cACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgtg                                                                                             <  2:16444/148‑1 (MQ=255)                                                   gATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCgg                                                                                                                                         >  6:192962/1‑101 (MQ=255)                                                   gATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCgg                                                                                                                                         <  5:192962/101‑1 (MQ=255)                                                       aCCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggc                                                                                     <  1:45243/149‑1 (MQ=255)                                                             cTACCTCCACACCCAACGATTGACATGCTCTACGCCATTATCGCCCGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggccggc                                                                                 <  4:324953/147‑1 (MQ=255)                                                                  tCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGc                                                                             <  8:101458/146‑1 (MQ=255)                                                                  tCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGc                                                                             >  7:101458/1‑146 (MQ=255)                                                                     aaaCCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCgg                                                                       <  3:253061/149‑1 (MQ=255)                                                                             ggATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGAc                                                               <  6:212153/149‑1 (MQ=255)                                                                              gATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACa                                                              <  2:246582/149‑1 (MQ=255)                                                                                ttGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGc                                                            <  4:62923/149‑1 (MQ=255)                                                                                   aCATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAAc                                                         >  5:263570/1‑149 (MQ=255)                                                                                      tGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGAc                                                      >  4:134120/1‑149 (MQ=255)                                                                                             cGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgtggt                                                                                           >  2:235651/1‑105 (MQ=255)                                                                                             cGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgtggt                                                                                           <  1:235651/105‑1 (MQ=255)                                                                                                 aTTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGc                                           <  3:299273/149‑1 (MQ=255)                                                                                                   tATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCgg                                         >  5:4435/1‑149 (MQ=255)                                                                                                    aTCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCgg                                         <  4:276366/148‑1 (MQ=255)                                                                                                    aTCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCgg                                         <  6:263570/148‑1 (MQ=255)                                                                                                      cGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCgg                                         <  5:258384/146‑1 (MQ=255)                                                                                                      cGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCgg                                         >  6:258384/1‑146 (MQ=255)                                                                                                             gACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCgg                               <  7:59467/149‑1 (MQ=255)                                                                                                               cGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTa                             <  6:277028/149‑1 (MQ=255)                                                                                                                 tCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGc                           >  6:284547/1‑149 (MQ=255)                                                                                                                    cAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggc                                                                                                          >  6:281051/1‑67 (MQ=255)                                                                                                                    cAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggc                                                                                                          <  5:281051/67‑1 (MQ=255)                                                                                                                    cAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCt                         >  5:91020/1‑148 (MQ=255)                                                                                                                         tCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGt                   >  2:168658/1‑149 (MQ=255)                                                                                                                           ccTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggcc                                                                                    >  7:157866/1‑82 (MQ=255)                                                                                                                           ccTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggcc                                                                                    <  8:157866/82‑1 (MQ=255)                                                                                                                            cTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGc                                                                      <  7:4380/95‑1 (MQ=255)                                                                                                                            cTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGc                                                                      >  8:4380/1‑95 (MQ=255)                                                                                                                              ggAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGcc               >  7:16090/1‑148 (MQ=255)                                                                                                                              ggAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGcc               >  7:295323/1‑148 (MQ=255)                                                                                                                               gAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGACGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCg              >  8:299768/1‑148 (MQ=255)                                                                                                                                  aaGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAAtt          <  4:256103/149‑1 (MQ=255)                                                                                                                                   tgcgcCGGGCGGGCGGCCAGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGc                <  5:180696/141‑1 (MQ=255)                                                                                                                                   tgcgcCGGGCGGGCGGCCAGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGc                >  6:180696/2‑142 (MQ=255)                                                                                                                                    gcgcCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGc                                                            >  1:89453/1‑97 (MQ=255)                                                                                                                                    gcgcCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGc                                                            <  2:89453/97‑1 (MQ=255)                                                                                                                                    gcgcCGGGCGGGCGGCCAGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAAtt          >  1:194764/1‑147 (MQ=255)                                                                                                                                    gcgcCGGGCGGGCGGCCAGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAAtt          <  2:194764/147‑1 (MQ=255)                                                                                                                                          ggCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATCGc  >  1:280569/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                          CCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCATCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATCGC  >  NC_002947/5116364‑5116649

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