Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	936,221	G→C	6.4%	A3134G (GCC→GGC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	936,221	0	G	C	6.4%	228.1 / 2.4	71	A3134G (GCC→GGC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base C (3/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCACCACGC  >  NC_002947/936076‑936365                                                                                                                                                   |                                                                                                                                                 ggTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCg                                                                                                                                               <  5:79064/149‑1 (MQ=255)  gggAGCACGCGGCCATCGCCACTCAGTTCCAGTGCGCTCGCCTCAATGGGCGGGGTGATATCGGGGGCGGTGATCTGCGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACCACATCCAGGGTCTCGCCGCTGGC‑CTGGGCCgg                                                                                                                                              <  8:169936/147‑1 (MQ=255)  gTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCTACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCgg                                                                                                                                              >  5:319882/1‑149 (MQ=255)      gCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGcc                                                                                                                                                >  2:79278/1‑143 (MQ=255)      gCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGcc                                                                                                                                                <  1:79278/143‑1 (MQ=255)              ccATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGa                                                                                                                                   <  2:350559/148‑1 (MQ=255)              ccATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGAc                                                                                                                                  <  1:155870/149‑1 (MQ=255)              ccATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGAc                                                                                                                                  >  5:75760/1‑149 (MQ=255)               cATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACg                                                                                                                                 <  2:127577/149‑1 (MQ=255)               cATCGCCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTCGGC‑CTGACCCCGGCCCAGGCCAACg                                                                                                                                 <  2:204773/149‑1 (MQ=255)                aTCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACgg                                                                                                                                >  1:320549/1‑149 (MQ=255)                       cTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCt                                                                                                                         <  6:319882/149‑1 (MQ=255)                        tCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTg                                                                                                                        >  1:264804/1‑149 (MQ=255)                           gTTGCAGTTCGCTGTCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCACCGTCGGC‑CTGACCCCGGCCCAGGCCAACggc                                                                                                                               <  8:314030/139‑2 (MQ=38)                           gTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCACCGTCGGC‑CTGACCCCGGCCCAGGCCAACggc                                                                                                                               >  7:314030/1‑138 (MQ=38)                              gCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATcc                                                                                                                  >  4:240639/1‑149 (MQ=255)                                 gTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGa                                                                                                               <  7:289336/149‑1 (MQ=255)                                         gCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCt                                                                                                        <  8:1983/148‑1 (MQ=255)                                         gCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTg                                                                                                       >  7:234165/1‑149 (MQ=255)                                               ttGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCt                                                                                                 >  7:70628/1‑149 (MQ=255)                                                tGGGGGGGTTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGTACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCt                                                                                                 <  2:264804/148‑1 (MQ=255)                                                    cGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCa                                                                                            <  3:240639/149‑1 (MQ=255)                                                       ggTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCg                                                                                                            <  5:331525/130‑1 (MQ=255)                                                       ggTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCg                                                                                                            >  6:331525/1‑130 (MQ=255)                                                        gTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGa                                                                                         >  3:170786/1‑148 (MQ=255)                                                        gTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGAc                                                                                        <  8:70628/149‑1 (MQ=255)                                                            tatCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTg                                                                                    <  6:75760/149‑1 (MQ=255)                                                             atCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGc                                                                                   >  3:51084/1‑149 (MQ=255)                                                                    gCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTgcgc                                                                            <  2:320549/149‑1 (MQ=255)                                                                    gCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGGTGGC‑CTGGGCCGGGGGCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACAtt                                                                                     <  3:340198/140‑1 (MQ=255)                                                                      ggTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGt                                                                          >  8:44984/1‑149 (MQ=255)                                                                                  aGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCa                                                              >  5:245815/1‑149 (MQ=255)                                                                                       gCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTgg                                                         <  8:234165/149‑1 (MQ=255)                                                                                               ttGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCt                                                   <  4:170786/147‑1 (MQ=255)                                                                                               ttGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGc                                                 <  4:51084/149‑1 (MQ=255)                                                                                                  ccAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGCCGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAg                                                                                           >  6:122509/1‑104 (MQ=255)                                                                                                  ccAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAg                                                                                           <  5:122509/104‑1 (MQ=255)                                                                                                  ccAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCgg                                              >  5:366957/1‑149 (MQ=255)                                                                                                  ccAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CGGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCgg                                              >  1:77231/1‑149 (MQ=255)                                                                                                  ccAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGA‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCgg                                              >  1:256886/1‑149 (MQ=255)                                                                                                   cAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGt                                             >  3:168754/1‑149 (MQ=255)                                                                                                     gcgcCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGCCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTGCCACCTGGGAATCTGCCGGtcc                                           >  3:244603/1‑147 (MQ=255)                                                                                                          gTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTAc                                      >  4:69827/1‑149 (MQ=255)                                                                                                              aCCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCgg                                   >  3:233645/1‑148 (MQ=255)                                                                                                               ccGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGa                                  >  6:160716/1‑148 (MQ=255)                                                                                                                cGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGaaa                                >  5:1378/1‑149 (MQ=255)                                                                                                                  ccACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACg                              >  3:178646/1‑149 (MQ=255)                                                                                                                  ccACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCACTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACg                              >  6:297875/1‑149 (MQ=255)                                                                                                                    aCGACATCCAGGGTCTCGCCGTTGGCGCTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGc                             >  1:192539/1‑149 (MQ=255)                                                                                                                     cGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCACTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTg                           <  5:297875/149‑1 (MQ=255)                                                                                                                       aCATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGa                         <  6:245815/149‑1 (MQ=255)                                                                                                                       aCATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGa                         <  6:1378/149‑1 (MQ=255)                                                                                                                        cATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAg                        <  4:233645/149‑1 (MQ=255)                                                                                                                          tCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc                      >  8:172337/1‑149 (MQ=255)                                                                                                                          tCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc                      >  7:226765/1‑149 (MQ=255)                                                                                                                           ccAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCt                            >  1:7990/1‑142 (MQ=255)                                                                                                                           ccAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCt                            <  2:7990/142‑1 (MQ=255)                                                                                                                            cAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGt                    <  7:172337/149‑1 (MQ=255)                                                                                                                             aGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTgg                                                         <  2:116209/111‑1 (MQ=255)                                                                                                                             aGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTc                   <  7:44984/149‑1 (MQ=255)                                                                                                                                gcttcaccGTCGGC‑CTGACCCCGGCCCAGGCCAACGGCGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCAcc                >  7:304817/7‑149 (MQ=37)                                                                                                                                     gCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGCAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACAGTGCCCTGCGCCTCACGCACTTCCACCAGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGcca           >  5:125478/1‑149 (MQ=255)                                                                                                                                     gCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGcca           <  3:249953/149‑1 (MQ=255)                                                                                                                                       cGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc                      >  6:55539/1‑136 (MQ=255)                                                                                                                                       cGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc                      <  5:55539/136‑1 (MQ=255)                                                                                                                                        gTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGccacca        >  5:365598/1‑149 (MQ=255)                                                                                                                                             c‑cTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCACCAcg   >  5:190468/1‑149 (MQ=255)                                                                                                                                               cTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTgg                                                         >  4:219301/1‑94 (MQ=255)                                                                                                                                               cTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTgg                                                         <  3:219301/94‑1 (MQ=255)                                                                                                                                               cTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCACCAcgc  <  4:178646/149‑1 (MQ=255)                                                                                                                                               cTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCACCAcgc  >  1:312782/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCACCACGC  >  NC_002947/936076‑936365

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