Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	942,008	A→C	5.0%	F1205C (TTC→TGC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	942,008	0	A	C	5.0%	241.7 / 2.4	78	F1205C (TTC→TGC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (31/43);  new base C (2/2);  total (33/45)
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.94e-01

AGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCACGGCAGGA  >  NC_002947/941861‑942151                                                                                                                                                     |                                                                                                                                                aGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGaa                                                                                                                                                <  2:358439/149‑1 (MQ=255) aGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGaa                                                                                                                                                <  7:131906/149‑1 (MQ=255)    cGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGATAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGt                                                                                                                                             <  7:99579/149‑1 (MQ=255)       aCGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGat                                                                                                                                          <  8:47073/149‑2 (MQ=255)        cGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTTCACCGGTCGCGGCGAGCCCGGCGCTACC‑GTGCAAGTgcgc                                                                                                                                         <  7:95057/149‑4 (MQ=17)         gTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGa                                                                                                                                           >  4:20353/1‑146 (MQ=255)         gTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGa                                                                                                                                           <  3:20353/146‑1 (MQ=255)                  ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                                                                                >  1:266813/1‑148 (MQ=255)                  ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                                                                                >  1:35917/1‑148 (MQ=255)                    cACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCaa                                                                                                                              <  7:19238/148‑1 (MQ=255)                    cACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGCCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCaa                                                                                                                              <  5:312221/148‑1 (MQ=255)                     aCCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAgg                                                                                                                            >  3:243771/1‑149 (MQ=255)                       cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGc                                                                                                                           >  4:132872/1‑148 (MQ=255)                        aTTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCa                                                                                                                         <  8:184399/149‑1 (MQ=255)                         ttGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAg                                                                                                                        <  6:202515/149‑1 (MQ=255)                             ccTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCa                                                                                                                    <  1:66162/149‑1 (MQ=255)                             ccTGGGCCGGCGACAGGGTGAGACTGAACAGGCCGTCGGCACCGACGACGCCGGTGCAAATTAGGGTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGCCCGGTGCAAGTGACGCCCTCGGCGACGGCCACGt                                                                                                                      >  7:226225/1‑147 (MQ=21)                                    cGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGtt                                                                                                             <  1:194004/149‑1 (MQ=255)                                     ggCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTtg                                                                                                            <  3:132872/149‑1 (MQ=255)                                        gACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGc                                                                                                         >  7:115581/1‑149 (MQ=255)                                             ggTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGtt                                                                                                             <  3:19445/140‑1 (MQ=255)                                             ggTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGtt                                                                                                             >  4:19445/1‑140 (MQ=255)                                             ggTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAggcg                                                                                                    >  8:42084/1‑149 (MQ=255)                                             ggTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGCGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAggcg                                                                                                    >  4:159494/1‑149 (MQ=255)                                              gTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGTCCAAGGCCAGGTCAGTCGGTTGTGCAggcgg                                                                                                   <  6:219946/149‑1 (MQ=255)                                                 aGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGgcag                                                                                                >  4:300683/1‑149 (MQ=255)                                                   gCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGgcagc                                                                                               >  4:178699/1‑148 (MQ=255)                                                    cTGAACAGGCCGTCGGCACCGGCGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTGGGCGCGACCGGTGAAAGTGACGCCATCGGCCAGGGCCAGGTCCGTGGGTTGTGCAggcggcggc                                                                                               >  8:55448/1‑147 (MQ=17)                                                     tGAACAGGCCGTCGGCACCGACGACGCCGCTGCCAATTACGTTGCCGGCACCATCGCGCACTTGCA‑CGCTACCGCCGGGCTCGACGCGACAGCAGAAAGTGCCTCCATCGGCCAAGGCCAGGTCAGTAGGATGTACAGGAGtcagatag                                                                                            >  4:257860/1‑145 (MQ=21)                                                      gAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGt                                                                                           >  5:19364/1‑149 (MQ=255)                                                         cAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCgg                                                                                        <  7:236348/149‑1 (MQ=255)                                                              cGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAg                                                                                   <  8:208174/149‑1 (MQ=255)                                                                  ggCACCGACGACGCCGGTGCCAATTAGGTGGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCTGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGa                                                                               <  8:226225/149‑1 (MQ=255)                                                                   gCACCGGCGTCGTCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGAc                                                                              >  1:211248/1‑149 (MQ=255)                                                                     aCCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGAcc                                                                             >  5:199357/1‑148 (MQ=255)                                                                       cgacgaCGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAg                                                                          >  7:339366/1‑149 (MQ=255)                                                                       cgacgaCGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAg                                                                          >  8:47778/1‑149 (MQ=255)                                                                        gacgacGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGc                                                                         <  5:320399/149‑1 (MQ=255)                                                                        gacgacGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGc                                                                         <  7:47778/149‑1 (MQ=255)                                                                            acgCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCTAGGTGACCGAGCCCgg                                                                     <  2:35917/149‑1 (MQ=255)                                                                                cGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTgg                                                                 <  2:210414/149‑1 (MQ=255)                                                                                  gTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGa                                                               <  2:266813/149‑1 (MQ=255)                                                                                   tGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGAt                                                              <  3:178699/149‑1 (MQ=255)                                                                                      cAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACg                                                           >  1:335575/1‑149 (MQ=255)                                                                                       aaTTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATAc                                                            >  1:116869/1‑147 (MQ=255)                                                                                       aaTTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATAc                                                            <  2:116869/147‑1 (MQ=255)                                                                                            ggTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGcc                                                      <  6:199357/148‑1 (MQ=255)                                                                                             gTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCcgg                                                    >  3:102263/1‑149 (MQ=255)                                                                                                gCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCcggcg                                                  <  2:335575/148‑1 (MQ=255)                                                                                                  cGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCa                                               >  6:36242/1‑149 (MQ=255)                                                                                                    gcagcaTCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTggg                                                                >  3:145960/1‑130 (MQ=255)                                                                                                    gcagcaTCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTggg                                                                <  4:145960/130‑1 (MQ=255)                                                                                                         aTCGCGCACTTGCA‑CGCTAGCGCCGGTCTCGCCGCGACCGGTGAAAGTGACGCCATCCGCCAATTCCAGGTCAGTCGGTTGTGCAGGCGGCTGCAATTCGGCCGAGGTGACCGAGACCGGTTGGGATACGTTGCCGGCGGCATCGGTGa                                        <  8:97911/149‑1 (MQ=17)                                                                                                          tCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAg                                       >  6:245564/1‑149 (MQ=255)                                                                                                          tCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAg                                       <  7:89648/149‑1 (MQ=255)                                                                                                            gcgcACTTGCA‑CGGTAGCGCCGGGCTCGCCGGGACCGGTGATAGTGACGCCATCTGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCTTCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGTGATACGTTGCCGGCGGCATCGGTGAgg                                      <  3:257860/148‑1 (MQ=17)                                                                                                            gcgcACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAgg                                      <  8:115581/148‑1 (MQ=255)                                                                                                             gtcACTTTCACCGGTCGCGGCGAGCCCGGCGCTACC‑GTGCAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGt                                                                    >  5:348609/3‑117 (MQ=14)                                                                                                             gtcACTTTCACCGGTCGCGGCGAGCCCGGCGCTACC‑GTGCAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGt                                                                    <  6:348609/115‑1 (MQ=14)                                                                                                                aCTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGAc                                 <  2:211248/149‑1 (MQ=255)                                                                                                                    gCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCAt                                              >  1:345334/1‑132 (MQ=255)                                                                                                                    gCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCAt                                              <  2:345334/132‑1 (MQ=255)                                                                                                                    gCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTg                              <  5:36242/148‑1 (MQ=255)                                                                                                                      a‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGc                             >  2:10684/1‑147 (MQ=255)                                                                                                                        cGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAg                          <  5:245564/149‑1 (MQ=255)                                                                                                                           tAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCGATCGGCCAAGGCCAGGGCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGCGGTGACCGAGCCCGGTTGTGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtg                       <  7:62530/149‑1 (MQ=255)                                                                                                                           tAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtg                       <  1:10684/149‑1 (MQ=255)                                                                                                                               gcCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtg                       <  5:346356/145‑1 (MQ=255)                                                                                                                               gcCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtg                       >  6:346356/1‑145 (MQ=255)                                                                                                                                       cgccgcGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCATCGTTc           <  4:102263/149‑1 (MQ=255)                                                                                                                                         ccgcGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGtt                                                         >  4:62724/1‑101 (MQ=255)                                                                                                                                         ccgcGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGtt                                                         <  3:62724/101‑1 (MQ=255)                                                                                                                                          cgcgACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCACg        <  3:159494/149‑1 (MQ=255)                                                                                                                                             gACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCACGGCa     <  6:19364/149‑1 (MQ=255)                                                                                                                                             gACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCACGGCa     <  8:339366/149‑1 (MQ=255)                                                                                                                                                cGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtg                       <  3:239941/128‑1 (MQ=255)                                                                                                                                                cGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtg                       >  4:239941/1‑128 (MQ=255)                                                                                                                                                cGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCACGGCAGGa  >  8:192443/1‑149 (MQ=255)                                                                                                                                                   tGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCAc         >  8:164899/1‑139 (MQ=255)                                                                                                                                                   tGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCAc         <  7:164899/139‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                AGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCA‑CGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTGCCGTCGTTCACGGCAGGA  >  NC_002947/941861‑942151

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