Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,917,689	A→G	10.9%	I30V (ATC→GTC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,917,689	0	A	G	10.9%	192.2 / 12.1	73	I30V (ATC→GTC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base G (5/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCCA  >  NC_002947/4917543‑4917833                                                                                                                                                   |                                                                                                                                                 cATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATc                                                                                                                                                <  6:283321/149‑1 (MQ=255) cATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATc                                                                                                                                                <  6:380014/149‑1 (MQ=255)   ttattaCACGCGCCGCCGCGCTGCCAATAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCa                                                                                                                                              <  7:15224/149‑1 (MQ=255)        acacGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGgcg                                                                                                                                          <  4:296994/148‑1 (MQ=255)                cgccgcGCTGCCAAGTTGCCGGGCTCGACAGGTTGTTTACAATCATGTCTGAAATCAATAGTCTCGGCGGACAATCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTggg                                                                                                                                 <  5:320784/149‑1 (MQ=255)                   ctcgcTGCCATGAGGCCGGGCGGGTCTGGTTGTATTGATTCATGACTGATCTCAATAGGCTCTGCGCAGAAACCTCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGACCCAGGCTGCGCCTggatca                                                                                                                              <  8:127429/147‑5 (MQ=255)                    gcgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGc                                                                                                                             >  3:232734/1‑149 (MQ=255)                      gcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGc                                                                                                                             >  8:118278/1‑147 (MQ=255)                        tGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGc                                                                                                                         >  5:37287/1‑149 (MQ=255)                          ccAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTc                                                                                                                       >  1:326308/1‑149 (MQ=255)                            aaGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                       >  1:355258/1‑131 (MQ=255)                            aaGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                       <  2:355258/131‑1 (MQ=255)                            aaGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCt                                                                                                                     >  8:318067/1‑149 (MQ=255)                                       cTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACcgc                                                                                                          <  7:318067/149‑1 (MQ=255)                                        tCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACcgc                                                                                                          >  3:300779/1‑148 (MQ=255)                                          gACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACcgc                                                                                                          <  4:300779/146‑1 (MQ=255)                                               gTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGcc                                                                                                   >  1:250326/1‑148 (MQ=255)                                                ttGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGa                                                                                                 <  1:348993/149‑1 (MQ=255)                                                ttGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGa                                                                                                 <  8:158308/149‑1 (MQ=255)                                                  gTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAgg                                                                                               >  1:368783/1‑149 (MQ=255)                                                    ataCAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGaa                                                                                                               <  7:174814/131‑1 (MQ=255)                                                    ataCAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGaa                                                                                                               >  8:174814/1‑131 (MQ=255)                                                           cATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGt                                                                                      >  6:37025/1‑149 (MQ=255)                                                            aTGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGtt                                                                                     <  8:308961/149‑1 (MQ=255)                                                             tGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTg                                                                                    >  7:313740/1‑149 (MQ=255)                                                              gACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCACGCCGAGGTTTTATCGTTGc                                                                                   >  2:290738/1‑149 (MQ=255)                                                               aCTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGc                                                                                   <  2:250326/148‑1 (MQ=255)                                                                cTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACAGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGt                                                                                 >  5:125755/1‑149 (MQ=255)                                                                 tGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGa                                                                                                                <  5:131484/117‑1 (MQ=255)                                                                    aaTCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGa                                                                             >  6:280565/1‑149 (MQ=255)                                                                        aaTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACggg                                                                         >  2:357442/1‑149 (MQ=255)                                                                         aTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGc                                                                        <  3:203112/149‑1 (MQ=255)                                                                          tAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGgc                                                                                                                                           >  6:219134/1‑81 (MQ=255)                                                                          tAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGgc                                                                                                                                           <  5:219134/81‑1 (MQ=255)                                                                                cGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACt                                                                 >  2:232043/1‑149 (MQ=255)                                                                                 ggCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACcgc                                                                                                          <  2:157226/107‑1 (MQ=255)                                                                                 ggCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGa                                                                   <  5:37025/146‑1 (MQ=255)                                                                                 ggCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACtt                                                                <  7:118278/149‑1 (MQ=255)                                                                                     cacaAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAg                                                            >  7:180409/1‑149 (MQ=255)                                                                                      acaAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGc                                                           <  2:368783/149‑1 (MQ=255)                                                                                        aaaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGACGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTg                                                         >  5:304511/1‑149 (MQ=255)                                                                                         aaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCa                                                                      <  4:44944/135‑1 (MQ=255)                                                                                         aaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCa                                                                      >  3:44944/1‑135 (MQ=255)                                                                                          aCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGt                                                       <  8:180409/149‑1 (MQ=255)                                                                                            cgcgATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGATGGTGc                                                     >  5:238741/1‑149 (MQ=255)                                                                                                   aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGc                                                                        <  6:122233/123‑1 (MQ=255)                                                                                                   aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGc                                                                        >  5:122233/1‑123 (MQ=255)                                                                                                   aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGc                                               <  1:357442/148‑1 (MQ=255)                                                                                                   aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCt                                              >  7:232532/1‑149 (MQ=255)                                                                                                    aCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTg                                             <  2:326308/149‑1 (MQ=255)                                                                                                     cAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTTGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGa                                            >  7:353097/1‑149 (MQ=255)                                                                                                        ccaagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATg                                          >  1:100016/1‑148 (MQ=255)                                                                                                          aagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAg                                       <  8:313740/149‑1 (MQ=255)                                                                                                           agccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAgg                                      >  3:288621/1‑149 (MQ=255)                                                                                                            gccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGc                                     >  6:203907/1‑149 (MQ=255)                                                                                                             ccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGc                                     >  5:47411/1‑148 (MQ=255)                                                                                                                   aGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGc                              <  6:37287/149‑1 (MQ=255)                                                                                                                    ggCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCg                             >  7:195918/1‑149 (MQ=255)                                                                                                                     gCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGt                            >  6:15840/1‑149 (MQ=255)                                                                                                                        ggCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCa                         <  5:280565/149‑1 (MQ=255)                                                                                                                                cGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGc                 <  4:288621/149‑1 (MQ=255)                                                                                                                                    cccaggctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCtt               >  1:249818/6‑147 (MQ=255)                                                                                                                                         gctgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTg                                                                                    <  8:361362/73‑1 (MQ=255)                                                                                                                                         gctgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTg                                                                                    >  7:361362/1‑73 (MQ=255)                                                                                                                                         gctgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAgg                                      >  7:30254/1‑119 (MQ=255)                                                                                                                                         gctgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAgg                                      <  8:30254/119‑1 (MQ=255)                                                                                                                                          ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAgg           <  2:363430/145‑1 (MQ=255)                                                                                                                                          ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAgg           >  1:363430/1‑145 (MQ=255)                                                                                                                                          ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCaa       >  4:251266/1‑149 (MQ=255)                                                                                                                                          ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACCTTCAGCTGGTGCTACGGCTGAt                                           >  6:283074/1‑113 (MQ=255)                                                                                                                                          ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACCTTCAGCTGGTGCTACGGCTGAt                                           <  5:283074/113‑1 (MQ=255)                                                                                                                                               ccTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgcca  <  5:203907/149‑1 (MQ=255)                                                                                                                                               ccTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgc    >  5:164904/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCCA  >  NC_002947/4917543‑4917833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑