Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,816	T→A	7.3%	Q177L (CAG→CTG)	PP_4911 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,582,816	0	T	A	7.3%	199.1 / 5.4	70	Q177L (CAG→CTG)	PP_4911	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (3/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGC  >  NC_002947/5582679‑5582964                                                                                                                                          |                                                                                                                                                     cAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCGGGTTGGCCAGCGGCTTCCCCGCCGTGCTCCGCGGTTCGGCCAGCGTCTGCCACAGTTGGCCCTGCGCCTGCCGCTGCGCCGACTCACGGCCCTCGGGCTGGTGGTCGCgg                                                                                                                                           >  8:117935/1‑149 (MQ=255)    cGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTgg                                                                                                                                        <  6:298043/149‑1 (MQ=255)    cGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTgg                                                                                                                                        <  3:149433/149‑1 (MQ=255)        cTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGt                                                                                                                                          >  4:283355/1‑143 (MQ=255)        cTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGt                                                                                                                                          <  3:283355/143‑1 (MQ=255)         tttttCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGgatga                                                                                                                                    >  6:138249/1‑148 (MQ=255)         tttttCGCTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGAAACAGTTCGACCTGAGCCTGCACCTGCGACTACTCACTGCCCTCGTGCTGCTCCACGAGGTGgatga                                                                                                                                    >  4:173884/1‑148 (MQ=255)                   ccTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCa                                                                                                                         >  8:160975/1‑149 (MQ=255)                   ccTCTTCGGCAATCGCCGGCTGGTTGACCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCa                                                                                                                         >  1:9969/1‑149 (MQ=255)                         cGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTACACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGc                                                                                                                      <  6:52949/146‑1 (MQ=255)                           gCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTGTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATg                                                                                                                 <  3:278296/149‑1 (MQ=255)                           gCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATg                                                                                                                 <  2:11650/149‑1 (MQ=255)                           gCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATg                                                                                                                 >  2:88004/1‑149 (MQ=255)                            cAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGCGTAACTCCCCTCGTGCGGCTGCTCGAGGTGGTTGATACGGTGTTCCCGCACATgg                                                                                                                >  3:48943/1‑149 (MQ=255)                                 gCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcg                                                                                                           >  6:36949/1‑149 (MQ=255)                                      cTGGTTGGCCAGCGGCTTCACCTCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGt                                                                                                      <  6:290331/149‑1 (MQ=255)                                         gTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCaa                                                                                                   <  4:122657/149‑1 (MQ=255)                                         gTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCCGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTc                                                                                                     >  1:57410/1‑147 (MQ=255)                                          ttGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAg                                                                                                  <  4:93032/149‑1 (MQ=255)                                            ggCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGc                                                                                                 <  5:36949/148‑1 (MQ=255)                                              ccAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGa                                                                                              >  7:18220/1‑149 (MQ=255)                                               cAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGAc                                                                                             <  7:117935/149‑1 (MQ=255)                                                    gCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCAc                                                                                        <  5:138249/149‑1 (MQ=255)                                                     cTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGgat                                                                                                                                      <  7:162755/102‑1 (MQ=255)                                                     cTTCACCGCCGTGCTCAGCTGTTCCGACAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCCATGTCgat                                                                                                                                      >  8:162755/1‑102 (MQ=255)                                                      ttCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTg                                                                                      >  6:163798/1‑149 (MQ=255)                                                        cACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGGCTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTggg                                                                                    >  4:312803/1‑149 (MQ=255)                                                        cACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTgg                                                                                     <  2:9969/148‑1 (MQ=255)                                                          ccgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACCGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCa                                                                                  >  1:214545/1‑149 (MQ=255)                                                           cgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCAt                                                                                 <  8:55858/149‑1 (MQ=255)                                                           cgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCAt                                                                                 >  7:55858/1‑149 (MQ=255)                                                                   tCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATg                                                                         >  3:26992/1‑149 (MQ=255)                                                                        tGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCgg                                                                    >  7:255562/1‑149 (MQ=255)                                                                             ggCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTgg                                                                <  7:142732/148‑1 (MQ=255)                                                                             ggCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTgg                                                                <  7:160975/148‑1 (MQ=255)                                                                               ccAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTg                                                                                      >  4:223465/1‑124 (MQ=255)                                                                               ccAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTg                                                                                      <  3:223465/124‑1 (MQ=255)                                                                                   cTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCTCTGGGCATCCAGGATGGCCGGCTGGTTCTTGc                                                         <  8:18220/149‑1 (MQ=255)                                                                                         cAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAAt                                                   >  4:86034/1‑149 (MQ=255)                                                                                         cAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAAt                                                   >  6:137396/1‑149 (MQ=255)                                                                                            cAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCa                                                <  8:316707/149‑1 (MQ=255)                                                                                               ttGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGc                                             <  5:163798/149‑1 (MQ=255)                                                                                                 ggCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgt                                           >  6:327515/1‑149 (MQ=255)                                                                                                  gCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtg                                          >  5:218253/1‑149 (MQ=255)                                                                                                   cccTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtgt                                         <  1:88004/149‑1 (MQ=255)                                                                                                    ccTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGt                                                               >  5:316680/1‑126 (MQ=255)                                                                                                    ccTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGt                                                               <  6:316680/126‑1 (MQ=255)                                                                                                        aGCCTGCACCTGCGACGGCTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATc                                     <  3:312803/148‑1 (MQ=255)                                                                                                        aGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGc                                             >  1:189199/1‑140 (MQ=255)                                                                                                        aGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGc                                             <  2:189199/140‑1 (MQ=255)                                                                                                         gcgtGCACCTGCGAAGACTCACTGCCCTGGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGa                                   <  3:173884/146‑1 (MQ=255)                                                                                                         gCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGa                                   <  4:26992/149‑1 (MQ=255)                                                                                                          ccTGCACCTGCGTCGACTAACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGac                                  <  6:201409/149‑1 (MQ=255)                                                                                                           cTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGAATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtg                                          >  5:194674/1‑140 (MQ=255)                                                                                                               aCCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAg                             >  5:1117/1‑149 (MQ=255)                                                                                                                        gacTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCAc                    <  6:218253/149‑1 (MQ=255)                                                                                                                             aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTc                 <  5:327515/147‑1 (MQ=255)                                                                                                                             aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCg               >  8:167464/1‑149 (MQ=255)                                                                                                                             aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCg               >  7:6778/1‑149 (MQ=255)                                                                                                                              tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGTTGGTTCTTGCCAGAATCCAGGCGtgtg                                        <  6:343523/119‑1 (MQ=255)                                                                                                                              tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGTTGGTTCTTGCCAGAATCCAGGCGtgtg                                        >  5:343523/5‑123 (MQ=255)                                                                                                                               cgCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCggg             <  2:214545/148‑1 (MQ=255)                                                                                                                               ccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGTGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCtt                                                           <  3:318269/100‑1 (MQ=255)                                                                                                                               ccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGTGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCtt                                                           >  4:318269/4‑103 (MQ=255)                                                                                                                                caccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGc            >  3:105436/3‑149 (MQ=255)                                                                                                                                 cccTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGc            <  6:1117/148‑1 (MQ=255)                                                                                                                                   cTCGAGCTGCTGCTCGAGGTGGATGATACGCTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGCGCATCCAGGATGGACGTCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGacaa                                >  5:228976/1‑126 (MQ=255)                                                                                                                                    tCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGTCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTg          >  5:41566/1‑147 (MQ=255)                                                                                                                                    tCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTg          >  6:107594/1‑147 (MQ=255)                                                                                                                                      gTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGc      <  5:137396/149‑1 (MQ=255)                                                                                                                                        tctgctgctCGAGGTGGATGATACGGTTCTCACGCACATGGACACGCACGTCAAGCAGACGCAACTGGACATCCATGATGGCCGGCTGGTTCTTGCAAGAATACAGGCGTGGCATCGCCAACAGCACTTCCACGTACGGGCTGAAGCCt    >  7:145292/2‑149 (MQ=255)                                                                                                                                          tgctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATc                                     >  2:210379/1‑114 (MQ=255)                                                                                                                                          tgctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATc                                     >  2:205314/1‑114 (MQ=255)                                                                                                                                          tgctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATc                                     <  1:210379/114‑1 (MQ=255)                                                                                                                                          tgctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATc                                     <  1:205314/114‑1 (MQ=255)                                                                                                                                          tgctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGc  >  2:224325/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                     CAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGC  >  NC_002947/5582679‑5582964

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