Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,305,611	T→C	6.2%	A148A (GCA→GCG)	PP_3777 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,305,611	0	T	C	6.2%	241.4 / 4.3	81	A148A (GCA→GCG)	PP_3777	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/45);  new base C (3/2);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTT  >  NC_002947/4305463‑4305745                                                                                                                                                     |                                                                                                                                       cAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGt                                                                                                                                        <  5:177549/149‑1 (MQ=255) cAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGt                                                                                                                                        <  6:254388/149‑1 (MQ=255)    tGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     <  2:122668/149‑1 (MQ=255)    tGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     >  5:158067/1‑149 (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                                                                                                                                     >  1:141429/1‑135 (MQ=255)                  ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     <  2:141429/135‑1 (MQ=255)                      gCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAg                                                                                                                   >  7:181066/1‑149 (MQ=255)                       cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGTTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTTTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                  <  8:227887/149‑1 (MQ=255)                         cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                >  7:114967/1‑149 (MQ=255)                         cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                >  2:226422/1‑149 (MQ=255)                              atCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCt                                                                                                            <  3:248400/148‑1 (MQ=255)                               tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAgg                                                                                                                               <  5:206334/128‑1 (MQ=255)                               tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAgg                                                                                                                               >  6:206334/1‑128 (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                                                                                                 <  6:243160/149‑1 (MQ=255)                                        cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  5:64318/149‑1 (MQ=255)                                        cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  3:108657/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                                                                                         <  1:54119/144‑1 (MQ=255)                                                     gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  2:54119/1‑144 (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                                                                                            <  1:61490/132‑1 (MQ=255)                                                              cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            >  2:61490/1‑132 (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                                                                                     <  7:178548/117‑1 (MQ=255)                                                                                    aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAgg                                                                                     >  8:178548/1‑117 (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                                                                             >  3:106454/1‑112 (MQ=255)                                                                                                 cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATc                                                                             <  4:106454/112‑1 (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)                                                                                                                                                     |                                                                                                                                       CAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTT  >  NC_002947/4305463‑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: ‑