Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,926	T→C	5.9%	F57L (TTC→CTC)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,926	0	T	C	5.9%	196.3 / 2.4	68	F57L (TTC→CTC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

CTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGT  >  NC_002947/3762780‑3763068                                                                                                                                                    |                                                                                                                                                 cTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                                                                                <  7:11235/149‑1 (MQ=255) cTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                                                                                >  2:144877/1‑149 (MQ=255) cTTCCCTGTTGGACGTGATGGTCAGCGTCTACGCAGGCTCCCACGGCCTGTGGACGGAGGAGCGACGTGGCATTATCGAGGCCGTGTACGGGCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                                                                                <  3:337744/149‑1 (MQ=255)     ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑Cg                                                                                                                                           >  1:203943/1‑149 (MQ=255)     ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACCGCCACGGCGCCGAGGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGtt                                                                                                                                                 <  1:268419/145‑1 (MQ=255)     ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACAGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGtt                                                                                                                                                 >  2:268419/1‑144 (MQ=255)            gCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCgg                                                                                                                                    <  2:272046/149‑1 (MQ=255)             cccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                   <  4:11748/149‑1 (MQ=255)             cccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCgg                                                                                                                                    >  8:324354/1‑148 (MQ=255)                     gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATg                                                                                                                            <  3:2173/148‑1 (MQ=255)                      tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATc                                                                                                                         <  1:144877/149‑1 (MQ=255)                          cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGTCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCgg                                                                                                                     >  1:114452/1‑149 (MQ=255)                          cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCgg                                                                                                                     <  3:366755/149‑1 (MQ=255)                           gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                    >  2:104083/1‑149 (MQ=255)                           gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                    >  8:216081/1‑149 (MQ=255)                            tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                    >  7:8823/1‑148 (MQ=255)                              aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCAc                                                                                                                 >  8:317348/1‑149 (MQ=255)                                  cAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTc                                                                                                             <  8:194063/149‑1 (MQ=255)                                     gcgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAT‑CATGTTCCCGGTCACgtcgaac                                                                                                          <  8:194223/149‑7 (MQ=255)                                      cgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAgc                                                                                                         <  4:176200/149‑1 (MQ=255)                                        cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAgcgc                                                                                                       <  5:38572/149‑1 (MQ=255)                                         ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGcc                                                                                                      >  6:69527/1‑149 (MQ=255)                                              gCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCg                                                                                                 <  3:114736/149‑1 (MQ=255)                                                cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                                                                               <  8:8823/149‑1 (MQ=255)                                                   tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGGTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                             <  5:34426/148‑1 (MQ=255)                                                    ggACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCAt                                                                                           <  7:216081/149‑1 (MQ=255)                                                     gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                          >  1:185049/1‑149 (MQ=255)                                                                 aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                              >  7:280277/1‑149 (MQ=255)                                                                 aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                              >  2:299392/1‑149 (MQ=255)                                                                     ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCa                                                                                                                  <  6:312517/110‑1 (MQ=38)                                                                     ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCa                                                                                                                  >  5:312517/1‑110 (MQ=38)                                                                      gCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                          >  5:20558/1‑148 (MQ=255)                                                                       cAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                        <  3:197743/133‑1 (MQ=255)                                                                       cAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                        >  4:197743/1‑133 (MQ=255)                                                                       cAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                        >  4:292052/1‑149 (MQ=255)                                                                         atatCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACCTTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                       <  1:104083/149‑1 (MQ=255)                                                                          tatCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAccc                                                                     >  7:268578/1‑149 (MQ=255)                                                                             cGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                  >  8:217266/1‑149 (MQ=255)                                                                               aGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGATCGCCTCGCGGATCATGCCGTGGTGGCACCACTTACCCCGCCt                                                                >  5:51578/1‑149 (MQ=255)                                                                                       tACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCAc                                                         <  7:317348/148‑1 (MQ=255)                                                                                         cGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                      <  8:280277/149‑1 (MQ=255)                                                                                            ccACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                   >  6:331465/1‑149 (MQ=255)                                                                                            ccACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                   <  2:114452/149‑1 (MQ=255)                                                                                             cACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑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                                             <  7:320744/136‑1 (MQ=255)                                                                                                                gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                             >  8:320744/1‑136 (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                                                              >  3:342109/1‑113 (MQ=255)                                                                                                                      aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                              <  4:342109/113‑1 (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                    >  1:37519/1‑149 (MQ=255)                                                                                                                            gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                    >  2:50827/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)                                                                                                                                                    |                                                                                                                                                 CTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGA‑GGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGT  >  NC_002947/3762780‑3763068

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