Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,139,896	G→T	6.3%	P76Q (CCG→CAG)	PP_5481 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,139,896	0	G	T	6.3%	193.8 / 2.4	63	P76Q (CCG→CAG)	PP_5481	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base T (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

GGTCGGTGATCCACTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTGCAGCAACT  >  NC_002947/2139753‑2140043                                                                                                                                                |                                                                                                                                                    ggTCGGTGATCCACTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAgg                                                                                                                                                >  6:65666/1‑149 (MQ=255)     ggTGATCCACTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAATACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCACGCTGACGGCGACCGGCAGGTTGc                                                                                                                                            >  3:119938/1‑149 (MQ=255)     ggTGATCCACTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGc                                                                                                                                            >  8:132921/1‑149 (MQ=255)           ccACTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCgcgc                                                                                                                                       <  6:73201/148‑1 (MQ=255)             aCTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCAcc                                                                                                                                    >  6:94460/1‑149 (MQ=255)             aCTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCAcc                                                                                                                                    >  1:316392/1‑149 (MQ=255)               tGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGAACCgg                                                                                                                                  >  3:189367/1‑149 (MQ=255)                  aaGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGAt                                                                                                                               >  8:190447/1‑149 (MQ=255)                     gTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGc                                                                                                                            >  4:207556/1‑149 (MQ=255)                         gCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAgg                                                                                                                        <  4:302704/149‑1 (MQ=255)                         gCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAgg                                                                                                                        <  3:207556/149‑1 (MQ=255)                           tcgtcgGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGgc                                                                                                                      >  6:104232/1‑149 (MQ=255)                             gtcgGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCggg                                                                                                                                 >  4:2746/1‑136 (MQ=255)                             gtcgGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCggg                                                                                                                                 <  3:2746/136‑1 (MQ=255)                              tcgGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGgccgc                                                                                                                   <  6:197078/149‑1 (MQ=255)                                 ggCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCAAGTAGGTGACGAAGAACACCGCCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCTCACCGGCATCGCCAGTgccgccgc                                                                                                                >  8:243451/1‑149 (MQ=255)                                  gCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGgccgccgc                                                                                                                >  4:262658/1‑148 (MQ=255)                                       gTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCa                                                                                                                          <  2:110789/133‑1 (MQ=255)                                       gTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCa                                                                                                                          >  1:110789/1‑133 (MQ=255)                                        tGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCcagcag                                                                                                          >  2:328806/1‑148 (MQ=255)                                            cggcggTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCt                                                                                                     >  1:83035/1‑149 (MQ=255)                                              gcggTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGc                                                                                                   <  1:328806/149‑1 (MQ=255)                                                ggTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCAt                                                                                                 <  6:291989/149‑1 (MQ=255)                                                        cATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGAt                                                                                           <  5:94460/147‑1 (MQ=255)                                                          tCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGa                                                                                       <  7:243451/149‑1 (MQ=255)                                                          tCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGa                                                                                       <  5:65666/149‑1 (MQ=255)                                                             gCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCg                                                                                                                                   >  3:251400/1‑102 (MQ=255)                                                             gCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCg                                                                                                                                   <  4:251400/102‑1 (MQ=255)                                                                aGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAgc                                                                                 >  3:347074/1‑149 (MQ=255)                                                                 ggCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAgcg                                                                                >  3:53030/1‑149 (MQ=255)                                                                    aCCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGcgccg                                                                             >  4:210073/1‑149 (MQ=255)                                                                        gTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGaa                                                                         <  3:262658/149‑1 (MQ=255)                                                                           aTGTAGGTGACGAAGAATACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACa                                                                       <  4:119938/148‑1 (MQ=255)                                                                                ggTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCca                                                                 <  5:104232/149‑1 (MQ=255)                                                                                  tGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCcaca                                                               >  2:73637/1‑149 (MQ=255)                                                                                      gaagaaCACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCa                                                           >  1:221358/1‑149 (MQ=255)                                                                                        agaaCACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCcagc                                                                                                            >  5:310798/1‑98 (MQ=255)                                                                                        agaaCACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCcagc                                                                                                            <  6:310798/98‑1 (MQ=255)                                                                                                ggcggcATGGGCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGc                                                 <  8:261348/149‑1 (MQ=255)                                                                                                   ggcATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATCGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCAc                                              <  1:73637/149‑1 (MQ=255)                                                                                                               gggTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGtt                                  >  7:70158/1‑149 (MQ=255)                                                                                                               gggTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGtt                                  <  2:83035/149‑1 (MQ=255)                                                                                                                 gTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCcaca                                                               >  8:21213/1‑118 (MQ=255)                                                                                                                 gTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGtt                                  >  3:24031/1‑147 (MQ=255)                                                                                                                      tCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGc                           >  6:239622/1‑149 (MQ=255)                                                                                                                        aGCCACACCAGGCTGACGGCGACCGGGAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCa                         <  6:32030/149‑1 (MQ=255)                                                                                                                        aGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGaa                                                                         >  3:249146/1‑101 (MQ=255)                                                                                                                        aGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGaa                                                                         <  4:249146/101‑1 (MQ=255)                                                                                                                          ccACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCa                         >  4:84332/1‑147 (MQ=255)                                                                                                                          ccACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCa                         >  5:247314/1‑147 (MQ=255)                                                                                                                           cacaCCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCagag                      <  7:190447/149‑1 (MQ=255)                                                                                                                                 aGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTggg                <  4:347074/149‑1 (MQ=255)                                                                                                                                   gtgcgcGGCAACCTGCCGGTCGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCa                         <  8:313969/133‑1 (MQ=255)                                                                                                                                   gtgcgcGGCAACCTGCCGGTCGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTATGGCCACCGAATGGCGGTTCAGGTGCCa                         >  7:313969/6‑138 (MQ=255)                                                                                                                                   gCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTc                                 >  6:160495/1‑130 (MQ=255)                                                                                                                                       gcGGCAACCTGCCGGTCGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACATCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGcc                          <  7:337634/132‑1 (MQ=255)                                                                                                                                       gcGGCAACCTGCCGGTCGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCAGCGAACAGCCCCACATCCATGGCTCTGGCCAACCAATGGCGTTTCATGTGcc                          >  8:337634/2‑133 (MQ=255)                                                                                                                                       aCGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGCGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTg          >  8:159917/1‑149 (MQ=255)                                                                                                                                        cGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTgc         <  6:247314/149‑1 (MQ=255)                                                                                                                                        cGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTgc         <  2:221358/149‑1 (MQ=255)                                                                                                                                        cGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGCGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTgc         <  7:159917/149‑1 (MQ=255)                                                                                                                                              ccGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTGCAGCAAc   <  5:239622/149‑1 (MQ=255)                                                                                                                                               cGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTGCAGCAACt  >  3:189181/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GGTCGGTGATCCACTGGAAGGTCAGCTCGTCGGGCAGGGTGCGCGGCGGTACCTGCATCAGCCAGGCACCGGTCATGTAGGTGACGAAGAACACCGGCGGCATGGTCAGCGGGTTGGTCAGCCACACCAGGCTGACGGCGACCGGCAGGTTGCCGCGCACCGGGATCGCCAGGGCCGCCGCCAGCAGCATCTGCATCGGGATGGGAATAAGCGCCGCGAACAGGCCCACACCCATGGCTCTGGCCACCGAATGGCGGTTCAGGTGCCAGAGGTTGGGATCGTGCAGCAACT  >  NC_002947/2139753‑2140043

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