Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	936,213	T→C	5.4%	N3137D (AAC→GAC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	936,213	0	T	C	5.4%	223.9 / 2.2	74	N3137D (AAC→GAC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (2/2);  total (35/39)
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.30e-01

CGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACC  >  NC_002947/936066‑936358                                                                                                                                                    |                                                                                                                                                  cgccgcGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATTTCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTg                                                                                                                                                  <  8:275953/149‑1 (MQ=255)  gccgcGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTg                                                                                                                                                  <  8:19833/148‑1 (MQ=255)  gccgcGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTg                                                                                                                                                  <  2:36241/148‑1 (MQ=255)     gcgACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCt                                                                                                                                              <  8:17877/149‑1 (MQ=255)     gcgACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCt                                                                                                                                              >  8:1371/1‑149 (MQ=255)      cgACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTg                                                                                                                                             <  1:201378/149‑1 (MQ=255)            gTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCgg                                                                                                                                       >  3:282039/1‑149 (MQ=255)             tGAGCATGCTGCCATCGGCACTCAGTTTCAGTTCGCTGGCTTCTTTGGGCGGGGTGATATCGGGGGGGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTTGCCTGGGCCggg                                                                                                                                      <  8:83870/149‑1 (MQ=255)             tGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCggg                                                                                                                                      <  6:213418/149‑1 (MQ=255)                    gCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGc                                                                                                                               >  8:247205/1‑149 (MQ=255)                         cATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACGGCCACGACATCCAGGGTCTCGCCGTTGGCCTg                                                                                                                                             >  7:174295/1‑130 (MQ=255)                         cATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACGGCCACGACATCCAGGGTCTCGCCGTTGGCCTg                                                                                                                                             <  8:174295/130‑1 (MQ=255)                         cATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACg                                                                                                                          >  4:265271/1‑149 (MQ=255)                         cATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACg                                                                                                                          >  1:106866/1‑149 (MQ=255)                         cATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTCGGCCTGACCCCGGCCCAGGCCAACg                                                                                                                          <  8:63575/149‑1 (MQ=255)                                     gTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCggg                                                                                                                                      <  6:277998/125‑1 (MQ=255)                                     gTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCggg                                                                                                                                      >  5:277998/1‑125 (MQ=255)                                      ttGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCAt                                                                                                             >  3:52472/1‑149 (MQ=255)                                         cAGTTCGCTGGCTTCATTGGGCGGGGTTATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCg                                                                                                          <  8:55964/149‑1 (MQ=255)                                         cAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCg                                                                                                          <  8:60602/149‑1 (MQ=255)                                         cAGTTCGCTGGCTTCATTGGGCGGGGTGATATAGGGGGCGGTCATCTCCGTAGGCAGCCAACTGTTGCCAGCGCCGTCAACCGCCACCACATCCAGGGGCTCGACGTTGGCCTGGGCCGGGGACAGGCCGACGGTGAAGCTGCCATCCg                                                                                                          >  1:325988/1‑149 (MQ=255)                                           gTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGa                                                                                                        <  7:1371/149‑1 (MQ=255)                                                 tGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCggg                                                                                                                                      <  2:242728/113‑1 (MQ=255)                                                 tGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGCCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCggg                                                                                                                                      >  1:242728/1‑113 (MQ=255)                                                   gCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTg                                                                                                >  3:290415/1‑149 (MQ=255)                                                    cTTCATTGGGCGGGGTGATATGGTGGGCGGTGATCTGGTTTGGCAACGAACTGTTGCCAGCGCCGTCAACCGCAACGACATCCAGGGTCTCGCGGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTgg                                                                                               <  7:139958/149‑1 (MQ=255)                                                     ttCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTggg                                                                                              >  5:22827/1‑149 (MQ=255)                                                       cATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGgc                                                                                             <  4:52472/148‑1 (MQ=255)                                                              cGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACg                                                                                                                          >  6:254058/1‑112 (MQ=255)                                                              cGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCa                                                                                     <  2:106866/149‑1 (MQ=255)                                                              cGGGGTGATAACGGGGGCGGTGATCTGGGTAGGCAGCTAACTGTTGCCAGCGCATTCAACCTCAACGACATCCAGGTTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCa                                                                                     <  2:192969/149‑1 (MQ=255)                                                               ggggTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAg                                                                                    >  5:35332/1‑149 (MQ=255)                                                                 ggTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCAGCGCCCAggc                                                                                  <  4:129446/149‑2 (MQ=255)                                                                       atCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGa                                                                                                        >  6:351403/1‑121 (MQ=255)                                                                       atCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGa                                                                                                        <  5:351403/121‑1 (MQ=255)                                                                              gCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCACTGCCCAGGACATTGCCCTgcgc                                                                     >  7:257144/1‑149 (MQ=255)                                                                                ggTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGt                                                                   <  3:265271/149‑1 (MQ=255)                                                                                ggTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGt                                                                   >  5:95877/1‑149 (MQ=255)                                                                                   gATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGATGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCAc                                                                >  5:41888/1‑149 (MQ=255)                                                                                         ggTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCACTGCCCAGGACATTGCCCTGCGCGTCACGCACt                                                           <  8:257144/148‑1 (MQ=255)                                                                                                 gCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCa                                                                                     >  1:23214/1‑114 (MQ=255)                                                                                                 gCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCa                                                                                     <  2:23214/114‑1 (MQ=255)                                                                                                  cGAACTTTTGCCTGCTGGTTCAACCGCCCCGACATCCAGGGGGTCGCCGTTGGCCTGGGCCGGGGTCAGTCCGAGGGTGATGCTGCCATCCGGACCGGCCTGGGGTCTGCCCAGGACATTGCCCTGAGCGTGACGCACTTCCACCTgg                                                  <  2:325988/148‑1 (MQ=255)                                                                                                  cGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTggg                                                 >  4:165501/1‑149 (MQ=255)                                                                                                  cGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTggg                                                 <  4:282039/149‑1 (MQ=255)                                                                                                   ctgcTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCGCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGt                                                <  6:339394/146‑1 (MQ=255)                                                                                                       tGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCt                                            >  8:54615/1‑149 (MQ=255)                                                                                                           gCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCg                                        >  6:232460/1‑149 (MQ=255)                                                                                                           gCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCg                                        >  4:155108/1‑149 (MQ=255)                                                                                                               gcgcCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCAc                                                      >  3:184076/1‑131 (MQ=255)                                                                                                               gcgcCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCAc                                                      <  4:184076/131‑1 (MQ=255)                                                                                                               gcgcCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGATGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGc                                    <  6:41888/149‑1 (MQ=255)                                                                                                                  ccGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTgc                                                                       <  3:302910/111‑1 (MQ=255)                                                                                                                  ccGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTgc                                                                       >  4:302910/1‑111 (MQ=255)                                                                                                                   cGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTa                                <  6:22827/149‑1 (MQ=255)                                                                                                                     tCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTAcc                              <  7:54615/149‑1 (MQ=255)                                                                                                                     tCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTAcc                              <  5:232460/149‑1 (MQ=255)                                                                                                                       aaCCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCg                             >  2:332571/1‑148 (MQ=255)                                                                                                                        aCCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGCTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCgg                            >  5:228160/1‑148 (MQ=255)                                                                                                                             cACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAAc                        <  3:111438/147‑1 (MQ=255)                                                                                                                             cACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAAc                        >  4:111438/1‑147 (MQ=255)                                                                                                                               cGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTg                    <  6:35332/149‑1 (MQ=255)                                                                                                                                   aTCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGc                >  8:120326/1‑149 (MQ=255)                                                                                                                                     ccAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc               >  7:258070/1‑148 (MQ=255)                                                                                                                                     ccAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc               <  8:258070/148‑1 (MQ=255)                                                                                                                                      cAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCg              <  6:95877/148‑1 (MQ=255)                                                                                                                                       aGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGt             >  1:83228/1‑148 (MQ=255)                                                                                                                                       aGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTc            >  7:47057/1‑149 (MQ=255)                                                                                                                                         ggTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTc            <  1:332571/147‑1 (MQ=255)                                                                                                                                            ctcGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCgg       >  6:289127/1‑149 (MQ=255)                                                                                                                                             tcaccGTCGGCCTGACCCCGGCCCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGt                                 <  6:304053/119‑1 (MQ=38)                                                                                                                                             tcaccGTCGGCCTGACCCCGGCCCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGt                                 >  5:304053/4‑122 (MQ=38)                                                                                                                                               gCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGc                                                                                                   <  6:70293/54‑1 (MQ=255)                                                                                                                                                 cGTCGGCCTGACCCCGGCCCAGGCCAACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGccacc  >  3:194545/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGCCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACC  >  NC_002947/936066‑936358

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