Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,854	G→T	8.3%	G115V (GGC→GTC)	PP_4562 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,854	0	G	T	8.3%	186.1 / 6.0	60	G115V (GGC→GTC)	PP_4562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base T (3/2);  total (32/28)
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.97e-01

CTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGC  >  NC_002947/5182709‑5183000                                                                                                                                                   |                                                                                                                                                   cTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGa                                                                                                                                                 <  7:77177/149‑1 (MQ=255)   gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGag                                                                                                                                                >  8:104922/1‑148 (MQ=255)    cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGagc                                                                                                                                               <  7:191680/148‑1 (MQ=255)      gcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGagcagc                                                                                                                                            >  1:167632/1‑149 (MQ=255)              aGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGGGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCGGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTc                                                                                                                                    <  7:36488/149‑1 (MQ=255)                    ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAAc                                                                                                                              >  2:185955/1‑149 (MQ=255)                      accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGT‑GTCGcgca                                                                                                                                              >  1:158695/1‑131 (MQ=255)                      accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGT‑GTCGcgca                                                                                                                                              <  2:158695/131‑1 (MQ=255)                       ccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTc                                                                                                                           >  5:46753/1‑149 (MQ=255)                          aTTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTa                                                                                                                        <  2:167632/149‑1 (MQ=255)                            ttCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTa                                                                                                                        >  8:275382/1‑147 (MQ=255)                                  aaGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                                 >  2:294276/1‑148 (MQ=255)                                  aaGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                                 >  3:88418/1‑148 (MQ=255)                                       tCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGcc                                                                                                           >  6:203671/1‑149 (MQ=255)                                         aCCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGt                                                                                                         >  2:157272/1‑149 (MQ=255)                                                gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCtt                                                                                                  >  3:75515/1‑149 (MQ=255)                                                gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCtt                                                                                                  >  5:209404/1‑149 (MQ=255)                                                gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCtt                                                                                                  >  7:177411/1‑149 (MQ=255)                                                   aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtct                                                                                               <  7:104922/149‑1 (MQ=255)                                                    aTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGTGAGCAGCTg                                                                                                                                          >  5:116603/1‑105 (MQ=255)                                                    aTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTg                                                                                                                                          <  6:116603/105‑1 (MQ=255)                                                     tGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCg                                                                                             <  1:317360/149‑1 (MQ=255)                                                         aGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGc                                                                                         >  4:114505/1‑149 (MQ=255)                                                           tGCGCAAGCCTGGCAATGGCCGCTCCACCATCAGCCTGGCCCTCAGCGGCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGGGa                                                                                       <  8:293604/149‑1 (MQ=255)                                                             cgcAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGAt                                                                                      <  4:88418/148‑1 (MQ=255)                                                                         aaGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCg                                                                                    <  2:9353/138‑1 (MQ=255)                                                                         aaGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCg                                                                                    >  1:9353/1‑138 (MQ=255)                                                                         aaGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGAc                                                                         >  2:122565/1‑149 (MQ=255)                                                                            cgccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGAt                                                                      <  2:182463/149‑1 (MQ=255)                                                                             gccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATCCATGACGATc                                                                     <  8:177411/149‑1 (MQ=255)                                                                              ccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGAt                                                                      >  8:100697/1‑147 (MQ=255)                                                                                      ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                                            <  1:6971/149‑1 (MQ=255)                                                                                         tCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCtt                                                          >  2:161999/1‑148 (MQ=255)                                                                                                   ccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCg                                                >  1:265578/1‑149 (MQ=255)                                                                                                   ccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                           <  4:224504/105‑1 (MQ=255)                                                                                                   ccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                           >  3:224504/1‑105 (MQ=255)                                                                                                         cgcCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATc                                                                                     >  3:37045/1‑105 (MQ=255)                                                                                                         cgcCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATc                                                                                     <  4:37045/105‑1 (MQ=255)                                                                                                           ccGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGcc                                       <  4:75515/149‑1 (MQ=255)                                                                                                                   gACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                               <  6:46753/149‑1 (MQ=255)                                                                                                                    aCTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATgg                                                                              <  5:117886/101‑1 (MQ=255)                                                                                                                    aCTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATgg                                                                              >  6:117886/1‑101 (MQ=255)                                                                                                                     cTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGt                             >  7:57930/1‑149 (MQ=255)                                                                                                                      tCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCGTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt                            >  6:196677/1‑149 (MQ=255)                                                                                                                      tCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt                            >  5:33375/1‑149 (MQ=255)                                                                                                                        gctgctGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                          >  3:88106/1‑149 (MQ=255)                                                                                                                        gctgctGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                          <  1:294276/149‑1 (MQ=255)                                                                                                                          agctgctCGCCACCCGGACGGGT‑GTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGccc                        >  5:228419/2‑149 (MQ=255)                                                                                                                             tgcGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAccc                                                 <  1:265725/121‑1 (MQ=255)                                                                                                                             tgcGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAccc                                                 >  2:265725/1‑121 (MQ=255)                                                                                                                               ctcgccacCCGGACGGGT‑GTCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAg                              >  2:147161/6‑138 (MQ=255)                                                                                                                               ctcgccacCCGGACGGGT‑GTCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAg                              <  1:147161/133‑1 (MQ=255)                                                                                                                               cgcgACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg                     <  3:114505/147‑1 (MQ=255)                                                                                                                                 cgACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGACGTTGCCCTTGCTGa                 >  6:243753/1‑149 (MQ=255)                                                                                                                                   acacCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACg               <  2:265578/149‑1 (MQ=255)                                                                                                                                     acCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGtt             <  1:185955/149‑1 (MQ=255)                                                                                                                                      cccGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTg            <  6:209404/149‑1 (MQ=255)                                                                                                                                        cGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtc                                                                                                >  6:340997/1‑63 (MQ=255)                                                                                                                                        cGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtc                                                                                                <  5:340997/63‑1 (MQ=255)                                                                                                                                                t‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc  <  1:122565/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGT‑GGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGC  >  NC_002947/5182709‑5183000

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