Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,604,205	C→G	8.2%	P12R (CCG→CGG) ‡	PP_4076 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,604,205	0	C	G	8.2%	225.1 / 9.4	73	P12R (CCG→CGG) ‡	PP_4076	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base G (3/3);  total (36/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 = 6.68e-01

GGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGA  >  NC_002947/4604069‑4604352                                                                                                                                         |                                                                                                                                                     gggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑Cg                                                                                                                                         >  8:237631/1‑149 (MQ=255)  ggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑Cgg                                                                                                                                        <  7:38162/149‑1 (MQ=255)   gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑Cgg                                                                                                                                        >  5:164379/1‑148 (MQ=255)    tGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGa                                                                                                                                      <  5:167642/149‑1 (MQ=255)       tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGc                                                                                                                                   <  6:165876/149‑1 (MQ=255)          gggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          <  5:276658/140‑3 (MQ=255)          gggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          >  6:276658/1‑138 (MQ=255)          gggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGtc                                                                                                                                >  1:303171/1‑149 (MQ=255)                gTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGcc                                                                                                                           <  2:303171/148‑1 (MQ=255)                 tACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAg                                                                                                                         >  6:116548/1‑149 (MQ=255)                   ccAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAgg                                                                                                                        >  7:203271/1‑148 (MQ=255)                       gTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          >  7:26176/1‑125 (MQ=255)                       gTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          <  8:26176/127‑3 (MQ=255)                          tcgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATc                                                                                                                                            <  4:332528/122‑1 (MQ=255)                          tcgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATc                                                                                                                                            >  3:332528/1‑122 (MQ=255)                          tcgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCt                                                                                                                >  4:85054/1‑149 (MQ=255)                            gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                              >  1:352817/1‑149 (MQ=255)                            gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                              <  3:78981/149‑1 (MQ=255)                             ctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCg                                                                                                             >  1:336448/1‑149 (MQ=255)                              tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCg                                                                                                             >  6:213717/1‑148 (MQ=255)                                         ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAGGTTcgcg                                                                                                 <  2:352817/149‑1 (MQ=255)                                           gCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGc                                                                                               >  2:44191/1‑149 (MQ=255)                                              cTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCAt                                                                                            >  1:353222/1‑149 (MQ=255)                                               ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGc                                                                                                                                   >  7:193343/1‑109 (MQ=255)                                               ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGc                                                                                                                                   <  8:193343/109‑1 (MQ=255)                                               ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCAt                                                                                            >  3:187459/1‑148 (MQ=255)                                                tCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          <  1:33360/102‑3 (MQ=255)                                                tCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          >  2:33360/1‑100 (MQ=255)                                                 cAAGTTCTTGGGCGCGAACTTCGAAATTTGCCTGGTTGGGCGGCTCGTCTTGCATTGAGTGGTTTCTTCGCTCGCTTGATGGCGGACCCGTCCGCCATC‑CGTGAGGCGTCTCGCCAGGATAGCGCTTCGCTCTTGTTCGGGGCGATc                                                                                           <  2:249357/147‑1 (MQ=255)                                                  aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAg                                                                                        >  2:265078/1‑149 (MQ=255)                                                    gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTg                                                                                      <  2:336448/149‑1 (MQ=255)                                                           ggCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTc                                                                               <  5:116548/149‑1 (MQ=255)                                                            gcgcGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTc                                                                               >  7:152261/1‑148 (MQ=255)                                                            gcgcGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCg                                                                              <  7:237631/149‑1 (MQ=255)                                                             cgcgAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGa                                                                             <  4:187459/149‑1 (MQ=255)                                                                gAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAAc                                                                          >  7:45878/1‑149 (MQ=255)                                                                     tCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                     >  8:279902/1‑149 (MQ=255)                                                                      cGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                      <  8:45878/147‑1 (MQ=255)                                                                      cGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑AGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                      <  8:152261/147‑1 (MQ=255)                                                                        aaaTTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTg                                                                   >  7:252172/1‑148 (MQ=255)                                                                          aTTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTGTTTTTCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAAc                                                                <  1:301052/149‑1 (MQ=255)                                                                           tttCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCCCCTCGAAACGTTGCTGAACg                                                               >  7:158794/1‑149 (MQ=255)                                                                                 tGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGc                                                          >  4:127684/1‑148 (MQ=255)                                                                                      ggCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                     >  5:204586/1‑132 (MQ=255)                                                                                      ggCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                     <  6:204586/132‑1 (MQ=255)                                                                                       gCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACaa                                                                                                        <  6:34528/96‑1 (MQ=255)                                                                                       gCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACaa                                                                                                        >  5:34528/1‑96 (MQ=255)                                                                                        ccGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTg                                                  >  3:48796/1‑149 (MQ=255)                                                                                               gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCACGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTcc                                              >  1:274371/1‑147 (MQ=255)                                                                                               gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCag                                           <  3:85054/149‑1 (MQ=255)                                                                                                 cacGCAATGATTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCagag                                         <  8:203271/149‑1 (MQ=255)                                                                                                    gCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagcag                                      >  6:229569/1‑149 (MQ=255)                                                                                                      aaTGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGcc                                    >  2:65194/1‑149 (MQ=255)                                                                                                           gTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                      <  5:78561/110‑1 (MQ=255)                                                                                                           gTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                      >  6:78561/1‑110 (MQ=255)                                                                                                           gTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAg                               <  5:213717/149‑1 (MQ=255)                                                                                                               tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                           >  4:133841/1‑149 (MQ=255)                                                                                                               tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                           <  3:133841/149‑1 (MQ=255)                                                                                                               tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                           <  3:127684/149‑1 (MQ=255)                                                                                                                ttcttcGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTggg                          >  2:296669/1‑149 (MQ=255)                                                                                                                  cttcGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAATCGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGc                         <  2:353222/148‑1 (MQ=255)                                                                                                                   ttcGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGccc                                   <  3:223205/137‑1 (MQ=255)                                                                                                                   tacGAACGCTTGATGGCGGACCCGTCAGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTACGCACAAGTTCGCGGCCATCAAGTGCCACCACGAAACGAAGCTGAACGCTCGCCAAGGTTGCTCCCAGAGAAGccc                                   >  4:223205/3‑137 (MQ=255)                                                                                                                      gATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGt                                                                                                                                 <  8:24787/40‑1 (MQ=255)                                                                                                                      gATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGt                                                                                                                                 >  7:24787/1‑40 (MQ=255)                                                                                                                      gATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTccc                                             >  4:148673/1‑124 (MQ=255)                                                                                                                      gATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTccc                                             <  3:148673/124‑1 (MQ=255)                                                                                                                      gATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTgc                     <  1:44191/148‑1 (MQ=255)                                                                                                                              gATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc             <  8:252172/148‑1 (MQ=255)                                                                                                                                    ggACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAGCGGGc      <  6:164379/149‑1 (MQ=255)                                                                                                                                    ggACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGc      <  4:48796/149‑1 (MQ=255)                                                                                                                                     gACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCa     <  7:279902/149‑1 (MQ=255)                                                                                                                                      aCCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAt    <  2:274371/149‑1 (MQ=255)                                                                                                                                        ccGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCGAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGa  >  5:300346/1‑149 (MQ=255)                                                                                                                                        ccGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCCCACAACTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCCGAGCAGCCCCGAGCTGGGCCAGCGCGACTTCAACTGGCATgc  >  6:286005/1‑148 (MQ=255)                                                                                                                                         |                                                                                                                                                     GGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATC‑CGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGA  >  NC_002947/4604069‑4604352

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