Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,785,395	G→A	7.0%	S133S (TCG→TCA)	PP_5068 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,785,395	0	G	A	7.0%	167.0 / 2.4	57	S133S (TCG→TCA)	PP_5068	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base A (2/2);  total (29/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 = 7.67e-01

AACCCACGTCTGGGCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGCCATGTGC  >  NC_002947/5785252‑5785543                                                                                                                                                |                                                                                                                                                     aaCCCACGTCTGGGCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTcg                                                                                                                                                 >  4:275375/1‑149 (MQ=255)   cccACGTCTGGGCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTcgc                                                                                                                                                >  8:82811/1‑148 (MQ=255)    ccACGTCTGGGCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTcgc                                                                                                                                                >  8:260682/1‑147 (MQ=255)     cACGTCTGGGCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAg                                                                                                                                             <  5:72823/149‑1 (MQ=255)              gCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGcc                                                                                                                                    <  5:63925/149‑1 (MQ=255)                gACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTc                                                                                                                                  >  4:148299/1‑149 (MQ=255)                          gCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCt                                                                                                                         <  3:148299/148‑1 (MQ=255)                           ccTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTg                                                                                                                        >  4:275329/1‑148 (MQ=255)                             tgctgctgGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCAc                                                                                                                     >  7:254621/1‑149 (MQ=255)                              gctgctgGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACg                                                                                                                    >  5:28686/1‑149 (MQ=255)                              gctgctgGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACg                                                                                                                    >  4:74501/1‑149 (MQ=255)                              gctgctgGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACg                                                                                                                    >  8:286911/1‑149 (MQ=255)                                  ctgGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTAGTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCt                                                                                                                <  8:51838/149‑1 (MQ=255)                                    ggCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTgg                                                                                                              <  7:260682/149‑1 (MQ=255)                                     gCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTgg                                                                                                              <  1:131561/148‑1 (MQ=255)                                     gCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTgg                                                                                                              <  7:268040/148‑1 (MQ=255)                                      ccGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTggg                                                                                                             >  3:246827/1‑148 (MQ=255)                                          aGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATcg                                                                                                        >  8:44758/1‑149 (MQ=255)                                                  acgacgCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGc                                                                                                <  3:74501/149‑1 (MQ=255)                                                    gacgCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCg                                                                                              <  3:275375/149‑1 (MQ=255)                                                     acgCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCgg                                                                                             <  7:82811/149‑1 (MQ=255)                                                        cAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCa                                                                                          <  7:286911/149‑1 (MQ=255)                                                          ggCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACt                                                                                        >  5:35067/1‑149 (MQ=255)                                                                     ggACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAAt                                                                               <  4:246827/147‑1 (MQ=255)                                                                          tGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTc                                                                         >  5:144131/1‑148 (MQ=255)                                                                            cACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCa                                                                       <  8:254621/148‑1 (MQ=255)                                                                            cACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAg                                                                      <  6:35067/149‑1 (MQ=255)                                                                             aCCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGc                                                                     <  6:28686/149‑1 (MQ=255)                                                                               cTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCAATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGca                                                                    >  8:265144/1‑148 (MQ=255)                                                                               cTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGAGCTCGTCAATTCGCGAGTTTATCGACTACGAAGCCCTGCACGCCCTGGGTATCTAGACCAGACTGGCACTGTGCGTAATCGGTTCAAgaa                                                                    >  8:260978/1‑146 (MQ=255)                                                                                tCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGcaca                                                                  >  6:185691/1‑149 (MQ=255)                                                                                 cAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGATCCAGCacac                                                                 >  5:74915/1‑149 (MQ=255)                                                                                  aaGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACAcc                                                                >  7:263004/1‑149 (MQ=255)                                                                                    gggCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTtctc                                                                                                                                       <  8:265637/76‑1 (MQ=255)                                                                                    gggCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTtctc                                                                                                                                       >  7:265637/1‑76 (MQ=255)                                                                                    gggCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCgg                                                              >  7:68585/1‑149 (MQ=255)                                                                                            gCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGgcgc                                                      <  3:73481/149‑1 (MQ=255)                                                                                             ccAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTc                                                                         <  3:108501/129‑1 (MQ=255)                                                                                             ccAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTc                                                                         >  4:108501/1‑129 (MQ=255)                                                                                                 acacCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAAc                                                 <  7:44758/149‑1 (MQ=255)                                                                                                   acCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTc                                               >  7:146875/1‑149 (MQ=255)                                                                                                     ccTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGt                                             >  2:97211/1‑149 (MQ=255)                                                                                                             cGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAgcgc                                     >  8:242372/1‑149 (MQ=255)                                                                                                                aTGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGc                                  >  3:297054/1‑149 (MQ=255)                                                                                                                aTGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGc                                  <  8:263004/149‑1 (MQ=255)                                                                                                                  gggTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGgt                                                             >  7:97176/1‑120 (MQ=255)                                                                                                                  gggTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGgt                                                             <  8:97176/120‑1 (MQ=255)                                                                                                                   ggTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCa                                <  3:275329/148‑1 (MQ=255)                                                                                                                     tGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGc                             <  8:91479/149‑1 (MQ=255)                                                                                                                     tGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGc                             >  7:91479/1‑149 (MQ=255)                                                                                                                     tGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGc                             <  5:185691/149‑1 (MQ=255)                                                                                                                     tGACCGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCACTCGGTATCGCGACCAGCCGGGCACTGTGCGTAATCCGTTCCAGCACACCGGTGAGGCGCGAAACTCGTGAAAGCGCGGCCATGc                             >  5:119851/1‑149 (MQ=255)                                                                                                                         ggCCGTGCCGTGCTGCGCTCGTCGATTAGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCAGGCTGAc                         >  1:81180/1‑149 (MQ=255)                                                                                                                              tgccgtgcTGCGCTCGTCAATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTc                    <  7:265144/149‑1 (MQ=255)                                                                                                                              tgacgtgcTGCGCTCGTCAAATCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTc                    <  7:260978/146‑1 (MQ=255)                                                                                                                                  gtgcTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCAGGCTGACGCGTCTGgc                <  2:81180/149‑1 (MQ=255)                                                                                                                                         gcTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGc         <  8:146875/149‑1 (MQ=255)                                                                                                                                           tcgtcgATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGCCa       >  4:184791/1‑149 (MQ=255)                                                                                                                                             gtcgATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGCCAtg     >  7:165821/1‑149 (MQ=255)                                                                                                                                             gtcgATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGCCAt      <  6:74915/148‑1 (MQ=255)                                                                                                                                                gATTCGCGAGTTTCTCGCCTCCGACGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGCCATGTGc  <  8:68585/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     AACCCACGTCTGGGCGACGGCCGTGGCCTGCTGCTGGCCGAAGTGCTCAACGACGCAGGCGAACACTGGGACCTGCACCTCAAGGGCGCCGGCCAGACACCCTACTCCCGCATGGGTGACGGCCGTGCCGTGCTGCGCTCGTCGATTCGCGAGTTTCTCGCCTCCGAAGCCCTGCACGCCCTGGGTATCGCGACCAGCCGGGCACTGTGCGTAATCGGTTCCAGCACACCGGTGTGGCGCGAAACTCGTGAAAGCGCGGCCATGCTGACGCGTCTGGCGCAGAGCCATGTGC  >  NC_002947/5785252‑5785543

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