Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,930:1	+C	6.0%	coding (173/723 nt)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,930	1	.	C	6.0%	269.0 / 11.0	67	coding (173/723 nt)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (34/29);  new base C (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAT  >  NC_002947/3762784‑3763078                                                                                                                                                    |                                                                                                                                                      ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑Cg                                                                                                                                                     >  1:203943/1‑149 (MQ=255)        gCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCgg                                                                                                                                              <  2:272046/149‑1 (MQ=255)         cccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                             <  4:11748/149‑1 (MQ=255)         cccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCgg                                                                                                                                              >  8:324354/1‑148 (MQ=255)                 gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATg                                                                                                                                      <  3:2173/148‑1 (MQ=255)                  tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATc                                                                                                                                   <  1:144877/149‑1 (MQ=255)                      cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGTCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCgg                                                                                                                               >  1:114452/1‑149 (MQ=255)                      cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCgg                                                                                                                               <  3:366755/149‑1 (MQ=255)                       gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                              >  2:104083/1‑149 (MQ=255)                       gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                              >  8:216081/1‑149 (MQ=255)                        tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                              >  7:8823/1‑148 (MQ=255)                          aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCAc                                                                                                                           >  8:317348/1‑149 (MQ=255)                              cAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTc                                                                                                                       <  8:194063/149‑1 (MQ=255)                                 gcgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAT‑CATGTTCCCGGTCACgtcgaac                                                                                                                    <  8:194223/149‑7 (MQ=255)                                  cgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAgc                                                                                                                   <  4:176200/149‑1 (MQ=255)                                    cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAgcgc                                                                                                                 <  5:38572/149‑1 (MQ=255)                                     ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGcc                                                                                                                >  6:69527/1‑149 (MQ=255)                                          gCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCg                                                                                                           <  3:114736/149‑1 (MQ=255)                                            cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                                                                                         <  8:8823/149‑1 (MQ=255)                                               tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGGTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                                       <  5:34426/148‑1 (MQ=255)                                                ggACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCAt                                                                                                     <  7:216081/149‑1 (MQ=255)                                                 gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                                    >  1:185049/1‑149 (MQ=255)                                                             aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                        >  2:299392/1‑149 (MQ=255)                                                             aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                        >  7:280277/1‑149 (MQ=255)                                                                 ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCa                                                                                                                            >  5:312517/1‑110 (MQ=38)                                                                 ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCa                                                                                                                            <  6:312517/110‑1 (MQ=38)                                                                  gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                                    >  5:20558/1‑148 (MQ=255)                                                                   cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                  >  4:197743/1‑133 (MQ=255)                                                                   cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                  <  3:197743/133‑1 (MQ=255)                                                                   cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                                  >  4:292052/1‑149 (MQ=255)                                                                     atatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACCTTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                                 <  1:104083/149‑1 (MQ=255)                                                                      tatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAccc                                                                               >  7:268578/1‑149 (MQ=255)                                                                         cGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                            >  8:217266/1‑149 (MQ=255)                                                                           aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGATCGCCTCGCGGATCATGCCGTGGTGGCACCACTTACCCCGCCt                                                                          >  5:51578/1‑149 (MQ=255)                                                                                   tACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCAc                                                                   <  7:317348/148‑1 (MQ=255)                                                                                     cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                <  8:280277/149‑1 (MQ=255)                                                                                        ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                             <  2:114452/149‑1 (MQ=255)                                                                                        ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                             >  6:331465/1‑149 (MQ=255)                                                                                         cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                            >  3:264809/1‑149 (MQ=255)                                                                                                       tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                              >  7:18358/1‑149 (MQ=255)                                                                                                         aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                             <  6:19554/148‑1 (MQ=255)                                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                     >  1:154198/1‑106 (MQ=255)                                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                     <  2:154198/106‑1 (MQ=255)                                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                       >  8:320744/1‑136 (MQ=255)                                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                       <  7:320744/136‑1 (MQ=255)                                                                                                             cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCtt                                        <  7:324354/149‑1 (MQ=255)                                                                                                               aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                       <  3:350040/148‑1 (MQ=255)                                                                                                               aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAAGGGCCTTg                                       >  4:350040/1‑148 (MQ=255)                                                                                                                 aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                        <  4:342109/113‑1 (MQ=255)                                                                                                                 aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                        >  3:342109/1‑113 (MQ=255)                                                                                                                 aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCCGGAA‑CATGATACCGGTAACGGGCCGGCGCCGGGCGGGTCATGCCGTGCTGGCACCACTTAACCCCCCTGCGGCACATGCCGTGCTGCGGGGCAACCGGCCTTGGt                                     >  6:50160/1‑148 (MQ=255)                                                                                                                   ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                   >  2:35134/1‑148 (MQ=255)                                                                                                                      ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                               >  3:269936/1‑149 (MQ=255)                                                                                                                       gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                              >  2:50827/1‑149 (MQ=255)                                                                                                                       gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                              >  1:37519/1‑149 (MQ=255)                                                                                                                            gATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCCCGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                          >  7:38630/1‑148 (MQ=255)                                                                                                                                aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGTCAACGGGCCTTGGTCACAGCGCCCGGACaa                     <  1:35134/149‑1 (MQ=255)                                                                                                                                  gcgcgcCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                    >  5:32372/1‑148 (MQ=255)                                                                                                                                      gcCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg               <  6:20558/149‑1 (MQ=255)                                                                                                                                      cgccaggCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCAC‑GTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg               >  4:45687/6‑149 (MQ=37)                                                                                                                                       ccTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg              <  8:268578/149‑1 (MQ=255)                                                                                                                                        cTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg             <  6:32372/149‑1 (MQ=255)                                                                                                                                          gggCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGAGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt            >  3:126177/1‑148 (MQ=255)                                                                                                                                              gTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                 >  7:175760/1‑123 (MQ=255)                                                                                                                                              gTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                 <  8:175760/123‑1 (MQ=255)                                                                                                                                                tCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCGGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg     >  5:196107/1‑149 (MQ=255)                                                                                                                                                 cGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg     <  1:299392/148‑1 (MQ=255)                                                                                                                                                   a‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt  <  8:221480/149‑1 (MQ=255)                                                                                                                                                   a‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt  <  4:269936/149‑1 (MQ=255)                                                                                                                                                   a‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt  >  8:3871/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAT  >  NC_002947/3762784‑3763078

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