Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,777,550	C→A	6.6%	Q675H (CAG→CAT)	PP_5065 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,777,550	0	C	A	6.6%	225.6 / 3.7	70	Q675H (CAG→CAT)	PP_5065	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base A (2/3);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGCCTTGAAAAGGGCCCGATGGGCGTC  >  NC_002947/5777407‑5777694                                                                                                                                                 |                                                                                                                                                 gCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTgcgc                                                                                                                                             >  1:15445/1‑149 (MQ=255)     cTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGt                                                                                                                                         <  2:15445/149‑1 (MQ=255)      tttCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTg                                                                                                                                        >  2:277651/1‑149 (MQ=255)      tttCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTg                                                                                                                                        >  2:164576/1‑149 (MQ=255)        tCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGTGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCAGGCGCCAggc                                                                                                                                        <  3:23499/147‑3 (MQ=255)        tCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTCC‑CCGCCTGGCGCAGGCGCCAggcg                                                                                                                                       <  6:315143/148‑4 (MQ=255)             tttGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAg                                                                                                                                 <  7:20804/149‑1 (MQ=255)                 cTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGca                                                                                                                              >  4:52272/1‑148 (MQ=255)                      gCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCAGGCGCCAggcg                                                                                                                                       >  3:183828/1‑131 (MQ=255)                      gCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCAGGCGCCAggcg                                                                                                                                       <  4:183828/134‑4 (MQ=255)                       cGGGCAGCGCATAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCg                                                                                                                       <  7:192405/149‑1 (MQ=255)                       cGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCg                                                                                                                       <  1:277651/149‑1 (MQ=255)                         ggCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCgg                                                                                                                      >  4:20495/1‑148 (MQ=255)                           cAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTaa                                                                                                                   <  3:100460/149‑1 (MQ=255)                            aGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTaaa                                                                                                                  <  5:117281/149‑1 (MQ=255)                              cgcAAAAAATGTCCGCAACAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCACTGCCTTCAGGCACACCTCGGTAAAGc                                                                                                                >  6:223789/1‑149 (MQ=255)                                 aacaactgtgCGCATCAACGACGCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGAt                                                                                                             <  6:224954/143‑1 (MQ=255)                                 aaaaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTgc                                                                                                                                               >  8:99621/1‑115 (MQ=255)                                 aaaaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTgc                                                                                                                                               <  7:99621/115‑1 (MQ=255)                                 aaaaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGcaca                                                                                                                            <  1:317526/134‑1 (MQ=255)                                 aaaaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGcaca                                                                                                                            >  2:317526/1‑134 (MQ=255)                                 aaaaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGa                                                                                                              >  7:132258/1‑148 (MQ=255)                                 aaaaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGa                                                                                                              >  6:276840/1‑148 (MQ=255)                                   aaaaTGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATaa                                                                                                           <  5:151500/149‑1 (MQ=255)                                         tgCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCTCCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATtgt                                                                                                       <  1:250183/148‑1 (MQ=255)                                                 aaCGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAAc                                                                                             <  1:164576/149‑1 (MQ=255)                                                  aCGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAAtt                                                                                                         <  3:232637/136‑1 (MQ=255)                                                  aCGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAAtt                                                                                                         >  4:232637/1‑136 (MQ=255)                                                          tcCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCAc                                                                                    <  5:276840/149‑1 (MQ=255)                                                          tcCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCAc                                                                                    >  3:109746/1‑149 (MQ=255)                                                               gACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCttt                                                                                <  5:223789/148‑1 (MQ=255)                                                               gACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCttt                                                                                <  1:47946/148‑1 (MQ=255)                                                                     cAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACtt                                                                         >  3:6799/1‑149 (MQ=255)                                                                       gTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTg                                                                        >  4:66835/1‑148 (MQ=255)                                                                          gCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCa                                                                      >  7:160079/1‑147 (MQ=255)                                                                          gCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCa                                                                      <  8:160079/147‑1 (MQ=255)                                                                          gCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAg                                                                     >  6:93980/1‑148 (MQ=255)                                                                          gCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAg                                                                     <  5:93980/148‑1 (MQ=255)                                                                            cccGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCTCGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAgg                                                                    >  1:239887/1‑147 (MQ=255)                                                                            cccGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGcc                                                                                         <  1:156227/126‑1 (MQ=255)                                                                            cccGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGcc                                                                                         >  2:156227/1‑126 (MQ=255)                                                                             ccGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAg                                                                 >  7:159041/1‑149 (MQ=255)                                                                              cGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAgc                                                                <  8:132258/149‑1 (MQ=255)                                                                                     gcCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCggc                                                         >  8:100491/1‑148 (MQ=255)                                                                                      ccTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCAGTTGCAGGTAGCGCCCggg                                                         >  2:319377/1‑148 (MQ=255)                                                                                      ccTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGc                                                        <  4:109746/149‑1 (MQ=255)                                                                                        tGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGcat                                                      >  5:39289/1‑149 (MQ=255)                                                                                          gCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGcatc                                                     <  6:39289/148‑1 (MQ=255)                                                                                           cGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAgcg                                                               >  5:131206/1‑137 (MQ=255)                                                                                           cGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAgcg                                                               <  6:131206/137‑1 (MQ=255)                                                                                              tCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCAtcatca                                                 >  5:91219/1‑148 (MQ=255)                                                                                                gAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCaaa                                               <  6:91219/148‑1 (MQ=255)                                                                                                 aacaacTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCg                                                                                        >  8:138476/1‑106 (MQ=255)                                                                                                 aacaacTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCtttt                                                                               <  2:267982/115‑1 (MQ=255)                                                                                                 aacaacTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCtttt                                                                               >  1:267982/1‑115 (MQ=255)                                                                                                 aacaacTGTCCTGTGCGCATTTAACGCAGGGCCTTC‑CGGCCTGGCGTCTGCGTGAGTGCGTTCATGCTCACCTCGGTAATGCGATAATTGTTCAGTTGAACGGCCg                                                                                        <  7:138476/106‑1 (MQ=255)                                                                                                      cTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCg                                        >  7:254663/1‑149 (MQ=255)                                                                                                         tCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCt                                      <  3:52272/148‑1 (MQ=255)                                                                                                         tCCTGGGCGCATTTCACGCTCGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTACAGCGATAATTGTGCCCGTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCtt                                     <  4:247394/149‑1 (MQ=255)                                                                                                            tgtgCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTtcat                                  >  4:143154/1‑149 (MQ=255)                                                                                                              tgCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTtcatc                                 <  4:6799/148‑1 (MQ=255)                                                                                                                 gcATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGt                             >  7:111507/1‑149 (MQ=255)                                                                                                                 gcATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGt                             <  3:143154/149‑1 (MQ=255)                                                                                                                      tCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGcc                         <  8:254663/148‑1 (MQ=255)                                                                                                                            acgGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGCCTTGaaaa                  <  3:66835/149‑1 (MQ=255)                                                                                                                               gCCTTC‑CCGCCTGGCGCATGCGCCAGTGCCTTCAGGCACACATCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGCCTTGAAAAggg               >  7:123592/1‑149 (MQ=255)                                                                                                                                 cTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCa                                                                            <  6:320678/86‑1 (MQ=255)                                                                                                                                 cTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCa                                                                            >  5:320678/1‑86 (MQ=255)                                                                                                                                         ccTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGCCTTGAAAAGGGCCCGATggg      <  3:20495/149‑1 (MQ=255)                                                                                                                                             gcgcCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGCCTTGAAAAGGGCCCGATGGGCGTc  >  5:221669/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 GCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTC‑CCGCCTGGCGCCTGCGCCAGTGCCTTCAGGCACACCTCGGTAAAGCGATAATTGTGCCCTTGAACGGCCGTCACCTTTTCCACTTGCAGGTAGCGCCCGGGCATCATCAAAGACTCCGCTTCATCAGGTCGCCTTGAAAAGGGCCCGATGGGCGTC  >  NC_002947/5777407‑5777694

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