Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,943,632	G→T	6.0%	R87R (CGC→CGA)	PP_3478 ←	secretion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,943,632	0	G	T	6.0%	216.0 / 3.3	67	R87R (CGC→CGA)	PP_3478	secretion protein
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base T (2/2);  total (31/36)
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.27e-01

GCCGCCTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGATGCCGGCCTCGTACTGGCCCTCTT  >  NC_002947/3943495‑3943780                                                                                                                                          |                                                                                                                                                      gccgccTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTc                                                                                                                                            <  7:321639/149‑1 (MQ=255)   cgccTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCGCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑c                                                                                                                                         >  1:66250/1‑149 (MQ=255)    gccTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGc                                                                                                                                                      <  8:158852/136‑1 (MQ=255)    gccTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGc                                                                                                                                                      >  7:158852/1‑136 (MQ=255)      cTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGt                                                                                                                                      >  3:263828/1‑149 (MQ=255)          ccAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGc                                                                                                                                  <  2:66250/149‑1 (MQ=255)          ccAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCGCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGc                                                                                                                                  <  1:212609/149‑1 (MQ=255)             aGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTg                                                                                                                               <  1:160945/149‑1 (MQ=255)               cgcTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGc                                                                                                                             <  8:8817/149‑1 (MQ=255)                    tGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAAt                                                                                                                         <  4:263828/148‑1 (MQ=255)                     gCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAAt                                                                                                                         >  5:35392/1‑147 (MQ=255)                     gCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAAt                                                                                                                         <  6:35392/147‑1 (MQ=255)                       gTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAgg                                                                                                                     >  2:267565/1‑149 (MQ=255)                       gTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCAGACGCGGTCG‑CTGTCGGCCTGGCCAATCAgg                                                                                                                     <  6:142003/149‑1 (MQ=255)                             ccAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTc                                                                                                                >  4:254900/1‑148 (MQ=255)                                  gTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGgggg                                                                                                                                                    <  2:221084/108‑1 (MQ=255)                                  gTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTAGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCgg                                                                                                                                                    >  1:221084/1‑108 (MQ=255)                                         tCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGAccgcc                                                                                                         <  8:40662/143‑1 (MQ=255)                                         tCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGAccgcc                                                                                                         >  7:40662/1‑143 (MQ=255)                                           aGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTggcg                                                                                                 <  7:173308/149‑1 (MQ=255)                                                  gCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTgg                                                                                          <  2:208769/149‑1 (MQ=255)                                                   cTGGCGAATGCCGTCCTGGTCGCGGCGGTTGATCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGTGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGTGCCTCGACCGCCGCCTGGCGGCTGTgg                                                                                          <  5:4800/148‑1 (MQ=255)                                                         aaTGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAggg                                                                                   <  2:111747/149‑1 (MQ=255)                                                           tGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCGCCTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGc                                                                                  <  1:267565/149‑1 (MQ=255)                                                           tGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCg                                                                                 <  8:203456/149‑1 (MQ=255)                                                             ccTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGa                                                                                >  5:161346/1‑148 (MQ=255)                                                              cTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGata                                                                              >  7:219771/1‑149 (MQ=255)                                                               ttCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGata                                                                              >  3:127873/1‑148 (MQ=255)                                                                tCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAgg                                                                                                                     >  3:12899/1‑108 (MQ=255)                                                                tCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAgg                                                                                                                     <  4:12899/108‑1 (MQ=255)                                                                tCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCg                                                                                 >  6:143647/1‑144 (MQ=255)                                                                tCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCg                                                                                 <  5:143647/144‑1 (MQ=255)                                                                 ccTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGata                                                                              <  4:127873/146‑1 (MQ=255)                                                                  cTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTggcg                                                                                                 >  6:145122/1‑126 (MQ=255)                                                                  cTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTggcg                                                                                                 <  5:145122/126‑1 (MQ=255)                                                                     ggCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGc                                                                        <  6:161346/148‑1 (MQ=255)                                                                      gcgcGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGctc                                                                      >  5:272264/1‑149 (MQ=255)                                                                        gcggcgGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTg                                                                    <  6:272264/149‑1 (MQ=255)                                                                          ggcggTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTa                                                                  <  3:254900/149‑1 (MQ=255)                                                                                 gttCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCAcgcg                                                           >  6:226595/1‑149 (MQ=255)                                                                                   tCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGa                                                         >  6:232705/1‑149 (MQ=255)                                                                                    cGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGAt                                                        >  4:126454/1‑149 (MQ=255)                                                                                        tCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCAGTATCACGCGGAtcgtc                                                    >  6:178897/1‑149 (MQ=255)                                                                                          cAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTgg                                                                                          <  4:313001/109‑1 (MQ=255)                                                                                          cAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTgg                                                                                          >  3:313001/1‑109 (MQ=255)                                                                                             ggTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATc                                                                <  4:60958/132‑1 (MQ=255)                                                                                             ggTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATc                                                                >  3:60958/1‑132 (MQ=255)                                                                                                gAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGAtc                                                       >  5:161237/1‑138 (MQ=255)                                                                                                gAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGAtc                                                       <  6:161237/138‑1 (MQ=255)                                                                                                  gCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGcc                                          >  7:251832/1‑149 (MQ=255)                                                                                                   cacaCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCct                                         >  7:120010/1‑149 (MQ=255)                                                                                                   cACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCct                                         >  1:150418/1‑149 (MQ=255)                                                                                                   cACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCct                                         >  3:131277/1‑149 (MQ=255)                                                                                                        ggCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTcc                                    <  8:219771/149‑1 (MQ=255)                                                                                                            gTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGgc                                <  8:120010/149‑1 (MQ=255)                                                                                                                cccATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGgcgagc                            >  5:101422/1‑149 (MQ=255)                                                                                                                 ccATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTc                                     <  4:322939/139‑1 (MQ=255)                                                                                                                 ccATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTc                                     >  3:322939/1‑139 (MQ=255)                                                                                                                   ggatgggCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGAt                         <  7:102783/145‑1 (MQ=255)                                                                                                                   aTCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGAtcgtc                                                    <  2:324940/122‑1 (MQ=255)                                                                                                                       gggCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGt                                                                   <  5:345137/103‑1 (MQ=255)                                                                                                                       gggCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGt                                                                   >  6:345137/1‑103 (MQ=255)                                                                                                                       gggCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGATGCCg                     <  5:178897/149‑1 (MQ=255)                                                                                                                                     tCATGTCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGgcga                              <  4:228249/126‑1 (MQ=255)                                                                                                                                     tCATGTCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGgcga                              <  4:232670/126‑1 (MQ=255)                                                                                                                                     tCATGTCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGgcga                              >  3:232670/1‑126 (MQ=255)                                                                                                                                     tCATGTCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGgcga                              >  3:228249/1‑126 (MQ=255)                                                                                                                                     tCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGATGCCGGCCTCGTACTGGcc       <  6:101422/149‑1 (MQ=255)                                                                                                                                         ggCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGATGCCGGCCTCGTACTGGCCctct   >  5:317345/1‑149 (MQ=255)                                                                                                                                          gCGGGCGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGATGCCGGCCTCGTACTGGCCCTCtt  >  3:196806/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                      GCCGCCTGGCCAAGCGCTATGCGTTCGTCCAGGGTGCGGATCAGCTCCAGCTGGCGAATGCCTTCCTGGGCGCGGCGGTTGTTCGGCTCCAGGGTGAGCACCCGGCCGTAGCCCATCCGGGCATTGGCGAAGTCATGGCGGACGCGGTCG‑CTGTCGGCCTGGCCAATCAGGGCCTCGACCGCCGCCTGGCGGCTGTGGGCCAGGGCGATATTCAGCTCTGTATCACGCGGATCGTCGCGCAGGGCCTCTTCCAGGCGAGCGATGCCGGCCTCGTACTGGCCCTCTT  >  NC_002947/3943495‑3943780

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