Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,350,084	C→T	8.0%	A853V (GCC→GTC) ‡	PP_4704 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,350,084	0	C	T	8.0%	179.7 / 4.1	63	A853V (GCC→GTC) ‡	PP_4704	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

TCAGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAG  >  NC_002947/5349937‑5350229                                                                                                                                                    |                                                                                                                                                   tCAGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCACCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGc                                                                                                                                                    >  7:252525/1‑148 (MQ=255)   aGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCACCAGCGCCTTGAGCGTGAGTACCAGGCCAc                                                                                                                                                 >  2:112155/1‑149 (MQ=255)             cccAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActca                                                                                                                                       >  7:70241/1‑146 (MQ=255)              ccAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActcac                                                                                                                                      <  8:70241/149‑5 (MQ=255)                aGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActcacgc                                                                                                                                    <  6:236377/149‑7 (MQ=255)                     tcCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCACCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCa                                                                                                                               <  8:252525/149‑1 (MQ=255)                       cTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGc                                                                                                                             >  7:105145/1‑149 (MQ=255)                             cAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACggg                                                                                                                       <  8:105145/149‑1 (MQ=255)                              aaGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACggg‑g                                                                                                                     <  3:352229/149‑1 (MQ=255)                                gCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑Gt                                                                                                                    >  4:287025/1‑148 (MQ=255)                                  tGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCAACGGG‑GTGa                                                                                                                  >  5:313864/1‑148 (MQ=255)                                    gCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGAc                                                                                                                 <  5:342503/147‑1 (MQ=255)                                     cATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAg                                                                                                              <  3:287025/149‑1 (MQ=255)                                       tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGCGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGt                                                                                                            >  1:187576/1‑149 (MQ=255)                                         cgcgCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGc                                                                                                          >  4:234086/1‑149 (MQ=255)                                                 aaGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCCACGGG‑GTGACAAGGTGCCGGCCAAt                                                                                                  <  6:313864/149‑1 (MQ=255)                                                     gcgcCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCt                                                                                              <  8:148663/149‑1 (MQ=255)                                                       gcCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTc                                                                                             >  5:268782/1‑148 (MQ=255)                                                              aCTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtgg                                                                                     <  2:187576/149‑1 (MQ=255)                                                                  cgcCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggt                                                                                 <  1:272387/149‑1 (MQ=255)                                                                    ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtg                                                                                >  3:197703/1‑148 (MQ=255)                                                                    ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑CTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                               >  2:135929/1‑149 (MQ=255)                                                                     cTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGAAAGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGga                                                                              >  1:297675/1‑149 (MQ=255)                                                                       ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                        >  2:182855/1‑106 (MQ=255)                                                                       ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                        <  1:182855/106‑1 (MQ=255)                                                                       ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGgac                                                                             >  5:286901/1‑148 (MQ=255)                                                                           aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtg                                                                                      <  4:296300/135‑1 (MQ=255)                                                                           aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtg                                                                                      >  3:296300/1‑135 (MQ=255)                                                                           aTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAg                                                                        <  6:286901/149‑1 (MQ=255)                                                                                    aCCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGAc                                                                    >  2:284814/1‑144 (MQ=255)                                                                                    aCCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGAc                                                                    <  1:284814/144‑1 (MQ=255)                                                                                     ccAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCaa                                                                                                   >  1:306699/1‑112 (MQ=255)                                                                                     ccAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCaa                                                                                                   <  2:306699/112‑1 (MQ=255)                                                                                        gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcg                                                             >  4:159427/1‑147 (MQ=255)                                                                                        gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcg                                                             <  3:159427/147‑1 (MQ=255)                                                                                             cAGGGCACCGTCGGGCCGCTGATTGCCCAACAGCGCCTTGAGCGTTATTTCCACGGCACCTTGTAAAATATCCTGCATCGTCGGGTGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGc                                                       <  1:106800/149‑1 (MQ=255)                                                                                                 ccAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑CTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGaa                                                  <  1:135929/149‑1 (MQ=255)                                                                                                   gtcgctggCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAgg                                                <  1:12861/143‑1 (MQ=255)                                                                                                      gACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAt                                             <  2:297675/149‑1 (MQ=255)                                                                                                       aCGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAtt                                            >  2:176238/1‑149 (MQ=255)                                                                                                                      cccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTAcc                                        <  4:265938/138‑1 (MQ=255)                                                                                                                      cccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTAcc                                        >  3:265938/1‑138 (MQ=255)                                                                                                                       ccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTCACCGTACCCGCCTg                             >  8:43274/1‑148 (MQ=255)                                                                                                                        cAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                     >  5:161007/1‑107 (MQ=255)                                                                                                                        cAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                     <  6:161007/107‑1 (MQ=255)                                                                                                                         aaCAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                          <  6:268782/149‑1 (MQ=255)                                                                                                                          aCAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                          >  6:99179/1‑148 (MQ=255)                                                                                                                             gcgcCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTTTATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagc                      <  1:112155/149‑1 (MQ=255)                                                                                                                              cgcCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagca                     <  5:99179/149‑1 (MQ=255)                                                                                                                                 cTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAgcc                  >  5:102305/1‑149 (MQ=255)                                                                                                                                        gTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGcc           >  5:146704/1‑149 (MQ=255)                                                                                                                                           tgTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGaa        >  4:170463/2‑149 (MQ=255)                                                                                                                                            gTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt       >  1:225603/1‑149 (MQ=255)                                                                                                                                            gTACCAGGCCACCTGGTACAAGCTCCTGCAGCGACGGG‑GTGACAACGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt       <  8:260976/149‑1 (MQ=255)                                                                                                                                            gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                          <  4:227441/130‑1 (MQ=255)                                                                                                                                            gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                          >  3:227441/1‑130 (MQ=255)                                                                                                                                            gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt       <  6:146704/149‑1 (MQ=255)                                                                                                                                             tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg      >  1:236915/1‑149 (MQ=255)                                                                                                                                             tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg      <  1:176238/149‑1 (MQ=255)                                                                                                                                               ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATgg     <  4:197703/148‑1 (MQ=255)                                                                                                                                                 aGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAg  <  2:236915/149‑1 (MQ=255)                                                                                                                                                  ggCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTCCCGTCCCCGCCTGAACCAGCAGCCGCCAGCCGAATGGGAg  >  7:312541/1‑148 (MQ=255)                                                                                                                                                   gCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                               <  4:108577/70‑1 (MQ=255)                                                                                                                                                   gCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                               >  3:108577/1‑70 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCAGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGG‑GTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAG  >  NC_002947/5349937‑5350229

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