Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,305,608	T→C	5.8%	L149L (CTA→CTG)	PP_3777 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,305,608	0	T	C	5.8%	261.6 / 4.1	86	L149L (CTA→CTG)	PP_3777	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/48);  new base C (3/2);  total (36/50)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTT  >  NC_002947/4305461‑4305745                                                                                                                                                    |                                                                                                                                          ggCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                          <  1:226445/149‑1 (MQ=255)  gCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAg                                                                                                                                         >  4:243012/1‑149 (MQ=255)  gCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAg                                                                                                                                         >  2:208030/1‑149 (MQ=255)   cAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGt                                                                                                                                        <  5:177549/149‑1 (MQ=255)   cAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGt                                                                                                                                        <  6:254388/149‑1 (MQ=255)      tGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     >  5:158067/1‑149 (MQ=255)      tGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     <  2:122668/149‑1 (MQ=255)        ccGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATc                                                                                                                                   >  4:110458/1‑149 (MQ=255)          gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGa                                                                                                                                 <  1:208030/149‑1 (MQ=255)          gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGa                                                                                                                                 <  8:36439/149‑1 (MQ=255)          gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                 <  3:110458/149‑1 (MQ=255)              aaCGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGAGTCCCGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGCggca                                                                                                                             >  3:156154/1‑147 (MQ=255)                    ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     <  2:141429/135‑1 (MQ=255)                    ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     >  1:141429/1‑135 (MQ=255)                        gCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAg                                                                                                                   >  7:181066/1‑149 (MQ=255)                         cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGTTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTTTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                  <  8:227887/149‑1 (MQ=255)                           cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                >  2:226422/1‑149 (MQ=255)                           cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                >  7:114967/1‑149 (MQ=255)                                atCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCt                                                                                                            <  3:248400/148‑1 (MQ=255)                                 tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAgg                                                                                                                               >  6:206334/1‑128 (MQ=255)                                 tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAgg                                                                                                                               <  5:206334/128‑1 (MQ=255)                                 tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGaa                                                                                                                             <  4:206581/130‑1 (MQ=255)                                 tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGaa                                                                                                                             >  3:206581/1‑130 (MQ=255)                                 tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGc                                                                                                          <  3:243012/149‑1 (MQ=255)                                   aGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCa                                                                                                        >  5:179031/1‑149 (MQ=255)                                    gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCAc                                                                                                       >  4:108657/1‑149 (MQ=255)                                     cTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCAcc                                                                                                      >  3:51380/1‑149 (MQ=255)                                        gCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTg                                                                                                   >  1:238869/1‑149 (MQ=255)                                          cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGTCAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                 >  7:222715/1‑117 (MQ=255)                                          cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGTCAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                 <  8:222715/117‑1 (MQ=255)                                          cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  3:108657/149‑1 (MQ=255)                                          cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  5:64318/149‑1 (MQ=255)                                          cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  6:243160/149‑1 (MQ=255)                                           tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                >  4:76771/1‑149 (MQ=255)                                           tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                >  8:240028/1‑149 (MQ=255)                                              gggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGACGTCTGCCACCCTGGCGGAt                                                                                             <  3:74607/149‑1 (MQ=255)                                               ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            <  4:238290/149‑1 (MQ=255)                                               ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            <  3:215019/149‑1 (MQ=255)                                                  ttCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATACAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  2:62614/1‑149 (MQ=255)                                                    cAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         <  4:51380/147‑1 (MQ=255)                                                    cAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAAt                                                                                        <  6:179031/148‑1 (MQ=255)                                                     aGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTNAATa                                                                                       <  5:3743/148‑1 (MQ=255)                                                       gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  2:54119/1‑144 (MQ=255)                                                       gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         <  1:54119/144‑1 (MQ=255)                                                       gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                    <  1:226422/149‑1 (MQ=255)                                                       gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                    <  3:205536/149‑1 (MQ=255)                                                       gTAGCGCGGCACCCGTTCGAAACACCGGCGCACGGTCCGCGCGCCATATAACAGCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                    <  4:156154/149‑1 (MQ=255)                                                           cTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATca                                                                                 >  2:112429/1‑148 (MQ=255)                                                              ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATc                                                                                                                                   >  4:22279/1‑95 (MQ=255)                                                              ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATc                                                                                                                                   <  3:22279/95‑1 (MQ=255)                                                              ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGAGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATcaca                                                                               <  2:148622/147‑1 (MQ=255)                                                                cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            >  2:61490/1‑132 (MQ=255)                                                                cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            <  1:61490/132‑1 (MQ=255)                                                                cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCg                                                                            >  5:230136/1‑148 (MQ=255)                                                                    cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTa                                                                          >  8:13777/1‑146 (MQ=255)                                                                    cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTa                                                                          <  7:13777/146‑1 (MQ=255)                                                                          aaaaaCCTGCGCAAGGTCGGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCctc                                                                 <  2:14952/149‑1 (MQ=255)                                                                                      aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAgg                                                                                     >  8:178548/1‑117 (MQ=255)                                                                                      aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAgg                                                                                     <  7:178548/117‑1 (MQ=255)                                                                                      aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAt                                                     >  7:85237/1‑149 (MQ=255)                                                                                       aGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  8:114967/149‑1 (MQ=255)                                                                                       aGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  3:76771/149‑1 (MQ=255)                                                                                       aGGTCCGCTCGCCATACAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  1:62614/149‑1 (MQ=255)                                                                                          tCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCg                                                 >  6:58637/1‑149 (MQ=255)                                                                                             gatagCCATAGAACATCGCTTACAGCGACAGCCCTTGGGCAGCCAAGATGCCCTGTAGATCATCGAGGAACCCATCCGAGCGCGTCTGCCACCCTCGCGTATTCAATAGGATCAAATCGAATACGACCTCTGCCCAAGGCATAGccag                                               >  4:82074/5‑145 (MQ=255)                                                                                             gCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGc                                                  <  3:82074/145‑1 (MQ=255)                                                                                                 gCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    >  7:69860/1‑139 (MQ=255)                                                                                                 gCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  8:69860/139‑1 (MQ=255)                                                                                                   cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATc                                                                             <  4:106454/112‑1 (MQ=255)                                                                                                   cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATc                                                                             >  3:106454/1‑112 (MQ=255)                                                                                                         aCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTccc                                  <  6:158067/149‑1 (MQ=255)                                                                                                          cATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCt                                 <  7:234851/149‑1 (MQ=255)                                                                                                          cATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCt                                 >  3:59260/1‑149 (MQ=255)                                                                                                                    aGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTc                       >  5:228510/1‑149 (MQ=255)                                                                                                                    aGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTc                       <  8:85237/149‑1 (MQ=255)                                                                                                                     gCGACAGCCATTGGTCAGCCTCGATGCGCTTTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCg                      <  2:238869/149‑1 (MQ=255)                                                                                                                      cGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGa                     <  7:240028/149‑1 (MQ=255)                                                                                                                           gCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAg                <  4:30496/149‑1 (MQ=255)                                                                                                                                 tggCAGCCTCGTTGCGCTGTTGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGATTCCGTCCTCTGCGGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTc           <  5:58637/147‑1 (MQ=255)                                                                                                                                 gggCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGAtt            <  4:59260/147‑1 (MQ=255)                                                                                                                                      tccTCGATGCACTACAGCGCATCGAGGCAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAg     >  3:137409/2‑149 (MQ=255)                                                                                                                                        cTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGt   <  6:228510/149‑1 (MQ=255)                                                                                                                                        cTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGt   <  1:217955/149‑1 (MQ=255)                                                                                                                                         tCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGtt  >  4:94440/1‑149 (MQ=255)                                                                                                                                                gcTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         <  4:248503/55‑1 (MQ=255)                                                                                                                                                gcTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  3:248503/1‑55 (MQ=255)                                                                                                                                                    |                                                                                                                                          GGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTT  >  NC_002947/4305461‑4305745

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