Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,917,966	T→C	8.2%	V122A (GTC→GCC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,917,966	0	T	C	8.2%	163.6 / 4.4	61	V122A (GTC→GCC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base C (2/3);  total (29/32)
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.53e-01

CAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTC  >  NC_002947/4917821‑4918114                                                                                                                                                   |                                                                                                                                                     cAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCgg                                                                                                                                                   >  4:54561/1‑149 (MQ=255) cAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCg                                                                                                                                                    <  1:290968/148‑1 (MQ=255)  aGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCgg                                                                                                                                                   <  5:304737/148‑1 (MQ=255)  aGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGc                                                                                                                                                  <  1:263063/149‑1 (MQ=255)    gCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGt                                                                                                                                                      <  3:159850/143‑1 (MQ=255)    gCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGt                                                                                                                                                      >  4:159850/1‑143 (MQ=255)    gCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCAc                                                                                                                                                >  6:206252/1‑149 (MQ=255)     ccAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCAcc                                                                                                                                               <  1:2253/149‑1 (MQ=255)     ccAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCAcc                                                                                                                                               >  5:197650/1‑149 (MQ=255)         cGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgc                                                                                                                                           <  3:28086/149‑1 (MQ=255)          gccagccaGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgct                                                                                                                                          >  1:65345/1‑149 (MQ=255)            cagccagcCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgctg                                                                                                                                         >  3:77618/1‑148 (MQ=255)             agccagccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgca                                                                                                                                       <  2:174005/149‑1 (MQ=255)                     gCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtg                                                                                                                                            >  1:241804/1‑136 (MQ=255)                     gCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtg                                                                                                                                            <  2:241804/136‑1 (MQ=255)                      cTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGcaaggcaagg                                                                                                                               <  1:84642/149‑1 (MQ=255)                         ccccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTc                                                                                                                             <  6:197650/147‑1 (MQ=255)                          cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                           >  7:195982/1‑148 (MQ=255)                          cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTg                                                                                                                          <  8:283708/149‑1 (MQ=255)                           ccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAATCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgctg                                                                                                                                         >  8:232818/1‑133 (MQ=255)                           ccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgctg                                                                                                                                         <  7:232818/133‑1 (MQ=255)                            cAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGac                                                                                                                        >  1:201883/1‑149 (MQ=255)                             aGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGacc                                                                                                                       <  7:106065/149‑1 (MQ=255)                                       gATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAgctgc                                                                                                             <  3:222648/149‑1 (MQ=255)                                           aCCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGc                                                                                                         >  7:70916/1‑149 (MQ=255)                                            ccGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGcc                                                                                                        >  1:13057/1‑149 (MQ=255)                                            ccGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGGCCTGACCAACCAGCTGCTGGcc                                                                                                        <  7:47908/149‑1 (MQ=255)                                                 aGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGTCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAgg                                                                                                    <  5:62377/148‑1 (MQ=255)                                                  gCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGac                                                                                                                        <  3:160809/127‑1 (MQ=255)                                                   ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAc                                                                                                 <  5:255313/149‑1 (MQ=255)                                                   ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAc                                                                                                 <  8:195982/149‑1 (MQ=255)                                                   ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAc                                                                                                 >  2:159271/1‑149 (MQ=255)                                                                aaTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGc                                                                                    >  4:228346/1‑149 (MQ=255)                                                                 aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGcc                                                                                   <  2:13057/149‑1 (MQ=255)                                                                 aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGcc                                                                                   >  4:51182/1‑149 (MQ=255)                                                                  tCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCg                                                                                  >  3:136123/1‑149 (MQ=255)                                                                  tCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCg                                                                                  >  6:53711/1‑149 (MQ=255)                                                                      cTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTgg                                                                                                          <  4:178919/121‑1 (MQ=255)                                                                      cTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTgg                                                                                                          >  3:178919/1‑121 (MQ=255)                                                                                tgcttcAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CCCCCAGCTGGACCCCCGCAAGGTGCCGGGCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                                    <  8:302500/149‑1 (MQ=255)                                                                                  ctgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTgg                                                                   <  8:70916/148‑1 (MQ=255)                                                                                     cAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCa                                                                <  4:77618/148‑1 (MQ=255)                                                                                      aaCAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGc                                                              <  6:41405/149‑1 (MQ=255)                                                                                          ggCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGc                                                          <  3:54561/149‑1 (MQ=255)                                                                                            ccagcgccagcCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTc                                                        <  1:159271/149‑1 (MQ=255)                                                                                                 gccagccaGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAg                                                   >  5:22527/1‑149 (MQ=255)                                                                                                 gccagacaGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCAGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGATGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGTGCGCTGGTAAGACTGCTCAACAg                                                   >  2:109718/1‑149 (MQ=255)                                                                                                  ccagccagGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGc                                                  <  4:136123/149‑1 (MQ=255)                                                                                                      ccagGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                               >  5:11098/1‑148 (MQ=255)                                                                                                      ccagGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                               >  6:33108/1‑148 (MQ=255)                                                                                                        agGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                    >  1:206445/1‑77 (MQ=255)                                                                                                        agGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                    <  2:206445/77‑1 (MQ=255)                                                                                                                 cccT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGc                                    >  8:323/1‑148 (MQ=255)                                                                                                                   cT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCcac                                 <  6:11098/149‑1 (MQ=255)                                                                                                                      cACCCAGCTGGACACCAGCAAGGTGCCGCTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCacac                               >  1:182891/1‑149 (MQ=255)                                                                                                                       aCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACt                              >  2:98710/1‑149 (MQ=255)                                                                                                                          cAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAg                           <  2:2509/149‑1 (MQ=255)                                                                                                                            gCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCa                         >  2:116699/1‑149 (MQ=255)                                                                                                                                 acacCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGAc                    >  5:47708/1‑149 (MQ=255)                                                                                                                                  cacCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACa                   <  5:206252/149‑1 (MQ=255)                                                                                                                                     cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcc                >  1:82506/1‑149 (MQ=255)                                                                                                                                                cGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTg     <  3:228346/149‑1 (MQ=255)                                                                                                                                                 ggTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCAGCCGTGTACCGGGCGCTGGTCAGCCTGCTCCACAGCACCCAGGCCGGTGCCACCCTGCGCATCGACAGCCCGCGCCCGCTgg    >  4:201490/1‑149 (MQ=255)                                                                                                                                                   tCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGt   >  2:83779/1‑148 (MQ=255)                                                                                                                                                   tCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCCACCGCCCGCGCCCGCTGGTc  <  2:182891/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCT‑CACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTC  >  NC_002947/4917821‑4918114

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