Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,079,182	G→C	6.9%	A28A (GCC→GCG)	PP_2688 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,079,182	0	G	C	6.9%	179.0 / 3.6	58	A28A (GCC→GCG)	PP_2688	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base C (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GCGCGGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCGGTTCATCTGCGTCTTTCAGGGTGCTGC  >  NC_002947/3079036‑3079330                                                                                                                                                   |                                                                                                                                                     gcgcGGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTg                                                                                                                                                      <  8:39793/147‑1 (MQ=255) gcgcGGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTg                                                                                                                                                      <  8:33946/147‑1 (MQ=255) gcgcGGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGATGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGgc                                                                                                                                                    <  3:77147/149‑1 (MQ=255)   gcgGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGgc                                                                                                                                                    >  4:36559/1‑147 (MQ=255)   gcgGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTCgccc                                                                                                                                                  >  5:353711/1‑148 (MQ=255)   gcgGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTCgccc                                                                                                                                                  >  5:348126/1‑148 (MQ=255)     ggATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGgccagc                                                                                                                                                <  5:333310/149‑1 (MQ=255)           tgcGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCgg                                                                                                                                          >  7:124472/1‑149 (MQ=255)           tgcGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCgg                                                                                                                                          >  4:551/1‑149 (MQ=255)            gcgcGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTg                                                                                                                                                      <  4:186205/136‑1 (MQ=255)            gcgcGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTg                                                                                                                                                      >  3:186205/1‑136 (MQ=255)                      atatCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTCGCCCGCTGGCTGGCCAACAccc                                                                                                                                <  6:348126/148‑2 (MQ=255)                      atatCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTCGCCCGCTGGCTGGCCAACAccc                                                                                                                                <  6:353711/148‑2 (MQ=255)                       tatCGAAGGGGATCAGGCGCTCGGTATGCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGa                                                                                                                              >  2:45143/1‑149 (MQ=255)                         tCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGAt                                                                                                                             >  2:308113/1‑148 (MQ=255)                            aaGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGCCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATg                                                                                                                         >  8:79949/1‑149 (MQ=255)                            aaGGGGATCAGGCGCTCGGTATCCTGCATGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGATCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGAGAACTCCTGATAATg                                                                                                                         >  6:9283/1‑149 (MQ=255)                              ggggATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATgcg                                                                                                                       <  7:79949/149‑1 (MQ=255)                                  aTCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATgcgggcg                                                                                                                   >  2:162771/1‑149 (MQ=255)                                     aGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATgcgggcgg                                                                                                                  >  2:19840/1‑147 (MQ=255)                                      tgcgcTCGGGTTCCTGCTTTTCGCCATACAAGGGGAAGGTGTTTCCTGCGCGGTGAAACAGGAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAg                                                                                                               <  7:231438/148‑1 (MQ=255)                                          atCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGc                                                                                                           >  6:340245/2‑149 (MQ=255)                                                      cTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGt                                                                                               <  3:551/149‑1 (MQ=255)                                                           cGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTc                                                                                                 >  6:223781/1‑142 (MQ=255)                                                           cGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTc                                                                                                 <  5:223781/142‑1 (MQ=255)                                                            gCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACAtt                                                                                          >  6:69815/1‑148 (MQ=255)                                                             ccATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCgg                                                                                                                                          <  3:236872/99‑1 (MQ=255)                                                             ccATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTc                                                                                        <  3:226544/149‑1 (MQ=255)                                                             ccATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGAGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCgg                                                                                                                                          >  4:236872/1‑99 (MQ=255)                                                              cATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCa                                                                                       <  1:45143/149‑1 (MQ=255)                                                               aTACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCAt                                                                                      >  1:148274/1‑149 (MQ=255)                                                                 aCAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCATCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTg                                                                                    <  7:179295/149‑1 (MQ=255)                                                                 aCAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTCTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCAtt                                                                                     >  5:190854/1‑148 (MQ=255)                                                                     ggtgaaggtgaTGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCa                                                                                       >  4:296894/1‑142 (MQ=255)                                                                     ggtgaaggtgaTGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCa                                                                                       <  3:296894/142‑1 (MQ=255)                                                                      gtgaaggtgaTGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCTACTCCAGCGGGGTGTTGTCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGaaaaa                                                                               <  1:162771/149‑1 (MQ=255)                                                                       tgaaggtgaTGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGaaaa                                                                                >  4:221522/1‑147 (MQ=255)                                                                         aaggtgaTGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGc                                                                            >  7:143004/1‑149 (MQ=255)                                                                            gtgaTGCCTGCGGGGTGTATCAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGGGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTGATTGAAAATAGCCCg                                                                         <  4:317429/149‑1 (MQ=255)                                                                                      cgcgGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATAcc                                                                <  3:36559/148‑1 (MQ=255)                                                                                        cgGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCAc                                                             <  5:69815/149‑1 (MQ=255)                                                                                           tGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGa                                                                      <  1:147935/137‑1 (MQ=255)                                                                                           tGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTcc                                                          >  5:189626/1‑149 (MQ=255)                                                                                              aaCAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTg                                                       >  2:278302/1‑149 (MQ=255)                                                                                               aCAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTg                                                       >  6:268814/1‑148 (MQ=255)                                                                                               aCAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTg                                                       <  5:268814/148‑1 (MQ=255)                                                                                                   cagGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTgccgcc                                                  <  1:278302/149‑1 (MQ=255)                                                                                                      gTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGtt                                               >  8:83191/1‑149 (MQ=255)                                                                                                      gTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGtt                                               >  7:286391/1‑149 (MQ=255)                                                                                                      gTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGtt                                               >  8:333090/1‑149 (MQ=255)                                                                                                           ccTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCgg                                          <  4:358166/149‑1 (MQ=255)                                                                                                            cTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCAGGGACATACCCACTCCTTGCCGCCGTTGCCgg                                          >  1:328596/1‑148 (MQ=255)                                                                                                               gcgcCGGCGCTGTTCCAGCAACTCCAGCGGGGTGATGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAg                                       <  8:124472/148‑1 (MQ=255)                                                                                                                cgcCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCt                                     <  2:37035/149‑1 (MQ=255)                                                                                                                      cgcTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTg                               <  5:340245/149‑1 (MQ=255)                                                                                                                        cTGTTTCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTg                               <  7:324154/147‑1 (MQ=255)                                                                                                                         tGTTCCAGCAACTCCAGCTGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATCCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAACCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCgg                            <  2:113131/149‑1 (MQ=255)                                                                                                                           ttCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCgg                            <  8:143004/147‑1 (MQ=255)                                                                                                                               aGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCCCATCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCGGTTCATc                      >  1:288175/1‑149 (MQ=255)                                                                                                                                    cTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCgg                            >  3:74506/1‑138 (MQ=255)                                                                                                                                    cTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGCCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCgg                            <  4:74506/138‑1 (MQ=255)                                                                                                                                                 ttGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCGGTTCATCTGCGTCTTTCAGGGtgct    >  8:283859/1‑149 (MQ=255)                                                                                                                                                   ggCCAGCCAGCGCGCGAACTCATGATAATCAGGGCGGCAGTTGCCGTCCGCTTCCTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGACGTTGCCGGAAGCATCATTGCGGTTCATCTGCGACTGTCAGGGtgatcc  >  8:288161/1‑145 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCGCGGATGCTGCGCGGGATGATATCGAAGGGGATCAGGCGCTCGGTATCCTGCTTGTCGCCATACAAGGTGAAGGTGATGCCTGCGCGGTGAAACAGCAGGTCGGCCTCGCGCCGGCGCTGTTCCAGCAACTCCAGCGGGGTGTTGGCCAGCCAGCGGGCGAACTCCTGATAATGCGGGCGGCAGTTGCCGTCCGCTTCGTACATTTCATTGAAAAAAGCCCGGGACATACCCACTCCTTGCCGCCGTTGCCGGCAGCTTCATTGCGGTTCATCTGCGTCTTTCAGGGTGCTGC  >  NC_002947/3079036‑3079330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑