Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,816,531	C→G	5.4%	A114G (GCC→GGC)	PP_5092 →	NLP/P60 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,816,531	0	C	G	5.4%	240.3 / 3.4	74	A114G (GCC→GGC)	PP_5092	NLP/P60 protein
Reads supporting (aligned to +/- strand):  ref base C (40/30);  new base G (2/2);  total (42/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 = 6.73e-01

GTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAC  >  NC_002947/5816386‑5816679                                                                                                                                                    |                                                                                                                                                     gTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGa                                                                                                                                                   >  2:145377/1‑149 (MQ=255)  tCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGAc                                                                                                                                                  >  1:165192/1‑149 (MQ=255)     ccgccATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTc                                                                                                                                               >  8:245863/1‑149 (MQ=255)             tCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATACCGCCTTCAACGACGTCGGCGACGTCgttgaa                                                                                                                                         <  1:240529/148‑6 (MQ=255)             tCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGcc                                                                                                                                       >  4:311971/1‑149 (MQ=255)                ccTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCcgcg                                                                                                                                    >  3:285199/1‑149 (MQ=255)                 cTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCcgcg                                                                                                                                    >  4:161786/1‑148 (MQ=255)                          gcagcGATGTGCTGAGCCGTGCGGTAAATGTGCGCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAAt                                                                                                                          <  8:273437/149‑1 (MQ=255)                             gcGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTCCCGCGCACCTCCAAAGc                                                                                                                        >  8:73874/1‑148 (MQ=255)                                        aGCCGTGCGGTAAATGTCCGCGGTAGTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAACGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCgggg                                                                                                            <  4:307210/149‑1 (MQ=255)                                               cGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     >  6:206119/1‑149 (MQ=255)                                               cGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  3:161786/149‑1 (MQ=255)                                               cGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  4:285199/149‑1 (MQ=255)                                                 gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGt                                                                                                   >  5:261217/1‑149 (MQ=255)                                                 gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGt                                                                                                   >  7:68219/1‑149 (MQ=255)                                                        tgCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTg                                                                                            <  2:165192/149‑1 (MQ=255)                                                               tACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAgg                                                                                     <  8:68219/149‑1 (MQ=255)                                                                    cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGGCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGa                                                                                 >  1:273527/1‑149 (MQ=255)                                                                    cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                >  6:306885/1‑149 (MQ=255)                                                                    cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                >  7:210832/1‑149 (MQ=255)                                                                    cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                <  3:311971/149‑1 (MQ=255)                                                                    cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGGGAAGGTGGCCAAGGGTGAt                                                                                >  3:56948/1‑149 (MQ=255)                                                                     ttATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                               >  5:228590/1‑149 (MQ=255)                                                                     ttATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                               >  4:17968/1‑149 (MQ=255)                                                                     ttATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                               >  3:90700/1‑149 (MQ=255)                                                                       aTGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                               <  1:145377/147‑1 (MQ=255)                                                                        tGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTg                                                                             >  4:161790/1‑148 (MQ=255)                                                                          tttGGGCCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTAAACGACGTCGCCGACGTCGACCTGCCGGGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTg                                                                             >  6:137184/1‑146 (MQ=255)                                                                                gcggcAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                <  1:253970/137‑1 (MQ=255)                                                                                gcggcAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                >  2:253970/1‑137 (MQ=255)                                                                                gcggcAGCAACCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                    >  4:45175/1‑149 (MQ=255)                                                                                   gcagcaGCCCTAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 <  3:239675/149‑1 (MQ=255)                                                                                      gcagcCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTg                                                              <  4:90700/149‑1 (MQ=255)                                                                                        agcCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAt                                                            >  5:351417/1‑149 (MQ=255)                                                                                         gcCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtc                                                           <  6:261217/149‑1 (MQ=255)                                                                                         acccGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtc                                                           <  3:45175/148‑1 (MQ=255)                                                                                          cccGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtct                                                          <  5:306885/149‑1 (MQ=255)                                                                                           ccGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtctt                                                         >  1:76456/1‑149 (MQ=255)                                                                                            cGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACg                                                                                                                                                 <  1:129113/60‑1 (MQ=255)                                                                                            cGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACg                                                                                                                                                 >  2:129113/1‑60 (MQ=255)                                                                                            cGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttc                                                        >  6:200251/1‑149 (MQ=255)                                                                                             gAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 >  7:250909/1‑139 (MQ=255)                                                                                             gAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 <  8:250909/139‑1 (MQ=255)                                                                                             gAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttct                                                       >  6:390627/1‑149 (MQ=255)                                                                                             gAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttct                                                       >  7:219590/1‑149 (MQ=255)                                                                                                 aaaGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATTCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCaa                                                   <  1:333/149‑1 (MQ=255)                                                                                                   aGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAATGTGGCCAAGGGTGATCTGAAGCATGGCGACCTTCTCGTCTTCAACa                                                 >  8:333206/1‑149 (MQ=255)                                                                                                         tGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCaa                                             <  5:200251/147‑1 (MQ=255)                                                                                                                cAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTg                                                              >  7:65269/1‑124 (MQ=255)                                                                                                                cAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTg                                                              <  8:65269/124‑1 (MQ=255)                                                                                                                cAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTc                                     >  8:114355/1‑149 (MQ=255)                                                                                                                cAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTc                                     <  2:273527/149‑1 (MQ=255)                                                                                                                 aGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCAGGGGGTGAAGGTGGCCAAGGGTGATCTGATGCCTGGCGACCTGATCTTCTTCAACATCATGAGCCGTCg                                    >  2:338396/1‑149 (MQ=255)                                                                                                                 aGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCCACCTGATCATCTTCACCATCAAGAGCCGTCg                                    >  2:344034/1‑149 (MQ=255)                                                                                                                  gCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCgg                                   <  3:161790/149‑1 (MQ=255)                                                                                                                          gTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCa                           >  2:362287/1‑149 (MQ=255)                                                                                                                           tCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAGCTTCTTCAACATCAAGAGCCGTCGGGTCAATCAc                          >  4:96223/1‑149 (MQ=255)                                                                                                                                   cGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAg                      >  6:116477/1‑146 (MQ=255)                                                                                                                                   cGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAg                      <  5:116477/146‑1 (MQ=255)                                                                                                                                    gCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACa                                                 <  3:337887/118‑1 (MQ=255)                                                                                                                                    gCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACa                                                 >  4:337887/1‑118 (MQ=255)                                                                                                                                       ttCAACGACGTCGCCGACGTCGCCCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTa               >  3:93266/1‑149 (MQ=255)                                                                                                                                        aggtcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc              >  3:139418/5‑149 (MQ=255)                                                                                                                                         ggtcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAcc             >  8:92203/4‑149 (MQ=255)                                                                                                                                         cAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc              <  3:298539/148‑1 (MQ=255)                                                                                                                                          aaCGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCt            >  7:56811/1‑149 (MQ=255)                                                                                                                                           acgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCACGAGCCGTCGGGTCAATCACGTAGCCATCTACCTg           <  8:322833/149‑1 (MQ=255)                                                                                                                                           acgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTg           >  2:249081/1‑149 (MQ=255)                                                                                                                                            cgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          <  8:210832/149‑1 (MQ=255)                                                                                                                                             gacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          <  4:139418/148‑1 (MQ=255)                                                                                                                                             gacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTggg         <  5:206119/149‑1 (MQ=255)                                                                                                                                             gacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGT‑AATCACGTAGCCATCTACCTGGGt        <  2:76456/149‑1 (MQ=255)                                                                                                                                               cgTCGCCGACGTCGACCTGCCTCTCACCTCCAATGCCATGGCCCGGGGCCACGGTGTGAATGTGGCAAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTa       <  3:295367/149‑1 (MQ=255)                                                                                                                                                 tCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAAc     <  7:114355/149‑1 (MQ=255)                                                                                                                                                    ccGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAc  <  8:56811/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACT‑GCAGCGGGCTGGTCAAATA‑CGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAC  >  NC_002947/5816386‑5816679

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑