Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	936,213	T→C	6.0%	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	6.0%	195.8 / 2.0	66	N3137D (AAC→GAC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCA  >  NC_002947/936076‑936359                                                                                                                                          |                                                                                                                                                    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                                                                                                                                          <  1:79278/143‑1 (MQ=255)      gCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGcc                                                                                                                                          >  2:79278/1‑143 (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                                                                                                      >  6:331525/1‑130 (MQ=255)                                                       ggTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCg                                                                                                      <  5:331525/130‑1 (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                >  7:226765/1‑149 (MQ=255)                                                                                                                          tCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGcc                >  8:172337/1‑149 (MQ=255)                                                                                                                           ccAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCt                      <  2:7990/142‑1 (MQ=255)                                                                                                                           ccAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCt                      >  1:7990/1‑142 (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)                                                                                                                                          |                                                                                                                                                    GGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTCTCGCCGTTGGC‑CTGGGCCGGGGTCAGGCCGACGGTGAAGCTGCCATCCGGACCGGCCTGGGCGCTGCCCAGGACATTGCCCTGCGCGTCACGCACTTCCACCTGGGTATCTGCCGGTGCGGTACCGGAAACGCTGGAGCCGTCACCGGCCACCA  >  NC_002947/936076‑936359

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