Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,816,730	G→C	6.3%	L180F (TTG→TTC)	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,730	0	G	C	6.3%	193.1 / 2.5	63	L180F (TTG→TTC)	PP_5092	NLP/P60 protein
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base C (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCG  >  NC_002947/5816586‑5816874                                                                                                                                                   |                                                                                                                                                 ggTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTCAGc                                                                                                                                               >  1:25016/1‑148 (MQ=255) ggTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCAACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTCAGc                                                                                                                                               >  1:25013/1‑148 (MQ=255)     gCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGc                                                                                                                                           <  8:7975/148‑1 (MQ=255)       cAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCt                                                                                                                                         >  7:129085/1‑148 (MQ=255)        aaGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGCGCAAGCCTTa                                                                                                                                       >  4:283385/1‑149 (MQ=255)         agtctGATCTGAAGCCTGGCGACCTGATCATCTTCAACATCAAGAGCCGTCTGGTCAATCACGTAGACATCTACCTGGGT‑AACGACCGCTTCATCAACGCTCAACCCC‑GAGGCAAGCGGGTGAGAATCGACACCTTGAGCAAGCCTTa                                                                                                                                       >  8:110232/5‑148 (MQ=255)              gATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTTGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGc                                                                                                                                  <  6:323692/148‑1 (MQ=255)                 cTGAAGCCTGGCGACCTGATCTTCTACAAAATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAg                                                                                                                                >  1:114952/1‑147 (MQ=255)                        cTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTAc                                                                                                                       >  7:40511/1‑149 (MQ=255)                          ggCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCa                                                                                                                           <  1:360211/143‑1 (MQ=255)                          ggCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCa                                                                                                                           >  2:360211/1‑143 (MQ=255)                           gCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAATCACTACgtg                                                                                                                    >  2:357665/1‑149 (MQ=255)                           gCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACgtg                                                                                                                    <  8:129085/149‑1 (MQ=255)                              aCCTGATCTTATTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACgtggtt                                                                                                                 <  7:110232/149‑2 (MQ=255)                                       tcttcAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCgg                                                                                                        >  5:120790/1‑149 (MQ=255)                                       tcttcAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCgg                                                                                                        >  8:42256/1‑149 (MQ=255)                                              cATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGtgctgc                                                                                                 >  7:95940/1‑149 (MQ=255)                                               aTCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAAACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGtgctgc                                                                                                 >  1:78060/1‑149 (MQ=255)                                                 cAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCg                                                                                               >  1:379711/1‑148 (MQ=255)                                                 cAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCATTACTGGAAGAAGCACTACGTGGTGGCCAACAGGGTGCTGACg                                                                                               >  8:347506/1‑148 (MQ=255)                                                  aaGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGa                                                                                              <  3:346784/148‑1 (MQ=255)                                                   agagCCGTCGGGTCAATCACGTAGCCTTCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGaaa                                                                                            <  4:274979/149‑1 (MQ=255)                                                   agagCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGaaa                                                                                            <  7:347506/149‑1 (MQ=255)                                                   agagCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGa                                                                                              >  3:62351/1‑147 (MQ=255)                                                     agCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAaga                                                                                          <  2:78060/149‑1 (MQ=255)                                                     agCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTCACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAaga                                                                                          <  8:80086/149‑1 (MQ=255)                                                        cGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGagca                                                                                       <  2:34342/149‑1 (MQ=255)                                                        cGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGagca                                                                                       <  2:379711/149‑1 (MQ=255)                                                        cGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGagca                                                                                       <  4:42132/149‑1 (MQ=255)                                                         gTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGagcag                                                                                      >  5:256547/1‑149 (MQ=255)                                                          tCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGagcagc                                                                                     <  2:107013/149‑1 (MQ=255)                                                                 aaTCACGTAGCCATCTACCTGTGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGTAGCAGCTGa                                                                               <  2:114952/148‑1 (MQ=255)                                                                 aaTCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCCGGGGGCAGGCGGGTGAGTATCGCCACCTTGAGCAAGCGTTAGTGGCAGAAGCACTACGTGGTGGCCATCCGGGTCCTGCCGAAAGAGAGGCAGCTgc                                                                               >  1:99349/1‑148 (MQ=255)                                                                 aaTCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGACCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAg                                                                              >  3:107657/1‑149 (MQ=255)                                                                                   cTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGc                                                            <  8:95940/149‑1 (MQ=255)                                                                                             gACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGgcg                                                   >  1:346738/1‑149 (MQ=255)                                                                                             gACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGgcg                                                   >  4:342686/1‑149 (MQ=255)                                                                                             gACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGATAAGCTGATCGgcg                                                   >  1:319129/1‑149 (MQ=255)                                                                                              aCCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGgcgc                                                  >  5:256158/1‑149 (MQ=255)                                                                                                cGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTa                                                >  2:40171/1‑149 (MQ=255)                                                                                                 gCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTAt                                               >  5:251425/1‑149 (MQ=255)                                                                                                  cTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATc                                              >  7:163656/1‑149 (MQ=255)                                                                                                      aTCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCgg                                           >  2:161095/1‑148 (MQ=255)                                                                                                          acgctccacgcC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCg                                      <  3:342686/149‑1 (MQ=255)                                                                                                            gctccacgcC‑GTGGCAAGCGGGTGAGTATCGACACCTTCAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGc                                    <  2:25013/149‑1 (MQ=255)                                                                                                            gctccacgcC‑GTGGCAAGCGCGTGAGTATCGACACCTTCAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGc                                    <  2:25016/149‑1 (MQ=255)                                                                                                              tccacgcC‑GTGGCAAGCGGGTGAGTATCGACACCATGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGccc                                  >  6:84082/1‑149 (MQ=255)                                                                                                                 acgcC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCtt                               >  7:256234/1‑149 (MQ=255)                                                                                                                 acgcC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCt                                >  6:336348/1‑148 (MQ=255)                                                                                                                       gTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg                          <  2:346738/149‑1 (MQ=255)                                                                                                                       gTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg                          <  7:42256/149‑1 (MQ=255)                                                                                                                            aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGc                                                                                                          <  7:339634/64‑1 (MQ=255)                                                                                                                            aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGc                                                                                                          >  8:339634/1‑64 (MQ=255)                                                                                                                              gCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCC                   >  2:249401‑M1/1‑149 (MQ=255)                                                                                                                                     aGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCC             <  2:99349‑M1/148‑1 (MQ=255)                                                                                                                                     aGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCC             <  6:256158‑M1/148‑1 (MQ=255)                                                                                                                                      gTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCA           <  5:352491‑M1/149‑1 (MQ=255)                                                                                                                                      gTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCAATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCA           >  3:315583‑M1/1‑149 (MQ=255)                                                                                                                                         tCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAG        <  8:40511‑M1/149‑1 (MQ=255)                                                                                                                                         tCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAG        >  5:148974‑M1/1‑149 (MQ=255)                                                                                                                                              acCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGTGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAgc   <  1:357665/149‑1 (MQ=255)                                                                                                                                              acCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAgc   <  6:120790/149‑1 (MQ=255)                                                                                                                                                cTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAtgt  >  2:228699‑M1/1‑145 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGT‑AACGACCGCTTCATCCACGCTCCACGCC‑GTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCG  >  NC_002947/5816586‑5816874

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