Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,149,779	G→C	9.3%	Q77H (CAG→CAC)	PP_3654 →	leucine‑responsive regulatory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,149,779	0	G	C	9.3%	214.9 / 12.2	75	Q77H (CAG→CAC)	PP_3654	leucine‑responsive regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (5/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

GGCTCAGCCTGTCACCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGACGAGCACCTGAC  >  NC_002947/4149631‑4149922                                                                                                                                                     |                                                                                                                                                    ggCTCAGCCTGTCACCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAg                                                                                                                                                     >  3:307770/1‑149 (MQ=255)       gCCGGTCACCCTGCCATCGGCGCTTTAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGTGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTg                                                                                                                                               <  4:301308/149‑1 (MQ=255)        ccTGTCACCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTgg                                                                                                                                              >  8:86186/1‑149 (MQ=255)          tGTCACCCTGACATCGGCGCTTGAAGGCACTGGAAGAGTCCGCAGCCATTCTGGACTATCACCCACGGCTCGCACCGAGCGCCCTTGCCCTCAACTTTGCCGCGCTGGTAATCGTGACGATACGCGAGGTAACCCGCCAGCCGGtctcg                                                                                                                                            >  2:96297/1‑145 (MQ=255)           gTCACCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCgg                                                                                                                                           >  7:350185/1‑149 (MQ=255)               cccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACtt                                                                                                                                       >  1:38530/1‑149 (MQ=255)               cccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACtt                                                                                                                                       >  2:69237/1‑149 (MQ=255)                ccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACttt                                                                                                                                      <  6:380109/149‑1 (MQ=255)                ccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACttt                                                                                                                                      <  6:192902/149‑1 (MQ=255)                ccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACttt                                                                                                                                      >  5:15844/1‑149 (MQ=255)                ccTGACATCGGAGCTTGAAGGCACTGGAAGAGTCCGGAGCCATACTGGGCTATCACGCAAGGCTGGCACCGAGCGCCCTAGGGCTGAACTTTGACGCGCTGCTATTCGTGACGCTACGCGAGGTAACCCGCCAGCAGGTGGCGGACttt                                                                                                                                      >  5:338069/1‑149 (MQ=255)                  tGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTg                                                                                                                                     >  2:107431/1‑148 (MQ=255)                      aTCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACg                                                                                                                               >  8:238518/1‑149 (MQ=255)                      aTCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACg                                                                                                                               >  3:52049/1‑149 (MQ=255)                                    gCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGgcg                                                                                                                      >  1:353739/1‑144 (MQ=255)                                      aCTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATcc                                                                                                              >  2:86503/1‑149 (MQ=255)                                      aCTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATcc                                                                                                              >  7:220954/1‑149 (MQ=255)                                           aaGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGAC‑TTGAA‑ACGGCACTGGCGCA‑GATCCCg                                                                                                            >  7:273420/1‑145 (MQ=255)                                            agagTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGa                                                                                                        <  5:3976/149‑1 (MQ=255)                                             gagTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGAt                                                                                                       >  5:126362/1‑149 (MQ=255)                                                tCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGt                                                                                                    <  2:38530/149‑1 (MQ=255)                                                   ggAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCg                                                                                                  >  6:242846/1‑148 (MQ=255)                                                   ggAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGCGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCg                                                                                                  >  8:80811/1‑148 (MQ=255)                                                   ggAGCCATGCTGGGCTTTCTCGCACGGGTGGCACCGAGCGCCCTTGGGATGAACTTTGCCGCGCTGTTATTCTTGACGCTACGCGAGGTAACCCTCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGa                                                                                                 <  1:96297/149‑1 (MQ=255)                                                    gAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCg                                                                                                  >  3:259159/1‑147 (MQ=255)                                                     aGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GAt                                                                                                                <  4:118580/132‑1 (MQ=255)                                                     aGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GAt                                                                                                                >  3:118580/1‑132 (MQ=255)                                                          ttCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCTTGACGCTACGCGAGGTAACGCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAg                                                                                          <  7:238518/149‑1 (MQ=255)                                                              gggCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGc                                                                                      >  4:85277/1‑149 (MQ=255)                                                              gggCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGc                                                                                      >  8:83159/1‑149 (MQ=255)                                                              gggCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGc                                                                                      >  6:266896/1‑149 (MQ=255)                                                                    tCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGTGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTc                                                                                 <  4:307770/148‑1 (MQ=255)                                                                    tCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAACACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTc                                                                                 <  1:86503/149‑1 (MQ=255)                                                                     cacgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCgg                                                                               >  8:24307/1‑149 (MQ=255)                                                                      acgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCgg                                                                               <  8:39535/148‑1 (MQ=255)                                                                       cgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGc                                                                                                                                                    >  8:72757/1‑80 (MQ=255)                                                                       cgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGc                                                                                                                                                    <  7:72757/80‑1 (MQ=255)                                                                       cgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCTCTGGCGCA‑GATCCCGCTGATTGTCGAAGCGCAGCGGCTGTTCGGCg                                                                             <  7:83159/149‑1 (MQ=255)                                                                       cgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCCGCGGCTGTTCGGc                                                                              >  2:269551/1‑148 (MQ=255)                                                                        gcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑a                                                                           <  4:259159/149‑1 (MQ=255)                                                                          acgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑Acc                                                                         >  8:183174/1‑149 (MQ=255)                                                                              cTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGa                                                                     >  1:149633/1‑149 (MQ=255)                                                                                    ccGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACct                                                               <  4:52049/149‑1 (MQ=255)                                                                                         cgcCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCtgctgc                                                          <  1:69237/149‑1 (MQ=255)                                                                                                gggCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCAtgtg                                                     <  8:350185/147‑1 (MQ=255)                                                                                                gggCTGAACTTTGCCGCGCTGGTATTCGCGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTgg                                                    <  7:80811/148‑1 (MQ=255)                                                                                                         tttGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGa                                          >  4:269048/1‑149 (MQ=255)                                                                                                         tttGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGa                                          >  7:127335/1‑149 (MQ=255)                                                                                                           tGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGa                                          <  2:11864/147‑1 (MQ=255)                                                                                                           tGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTAACTGCTGCATGTGGTAGCCAAGGAc                                         <  5:242846/148‑1 (MQ=255)                                                                                                             ccGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCt                                       <  7:183174/148‑1 (MQ=255)                                                                                                              cgcgCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGc                                     <  5:10433/149‑1 (MQ=255)                                                                                                                   tGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCg                               >  7:309437/1‑149 (MQ=255)                                                                                                                     gTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCg                               <  3:85277/147‑1 (MQ=255)                                                                                                                      tATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCACGGCgct                             <  6:126362/149‑3 (MQ=255)                                                                                                                        ttCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCa                          <  5:177083/149‑1 (MQ=255)                                                                                                                           gTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAg                         >  4:357251/1‑147 (MQ=255)                                                                                                                            tGAAGCTAGGCGATGTAACCGGCCAGCAGGTGGCGGACTTTGAA‑AGGGAAGTGGCGCA‑GATCCGGCAGTTTGTCGAAGCGCAGCGGCTGTTCGTCG‑ACCCGGATTAGCTGCTGCATGTGGTTCCCAAGGACCTGC‑CGGCGTTTCAGAAg                      <  6:338069/149‑1 (MQ=255)                                                                                                                              aCGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCt                    <  3:269048/149‑1 (MQ=255)                                                                                                                               cGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGACGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTg                   <  6:15844/149‑1 (MQ=255)                                                                                                                                   aCGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCAGGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTAt                >  1:279756/1‑149 (MQ=255)                                                                                                                                    cgcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTAt                <  1:107431/148‑1 (MQ=255)                                                                                                                                     gcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGAc             <  7:24307/149‑1 (MQ=255)                                                                                                                                     gcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGAc             >  4:137548/1‑149 (MQ=255)                                                                                                                                       gAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTa                                                  <  8:131372/111‑1 (MQ=255)                                                                                                                                       gAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTa                                                  >  7:131372/1‑111 (MQ=255)                                                                                                                                          aaagtccGCCACCGGCTGGCGGAGTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGACGAGCa        >  5:336044/6‑149 (MQ=255)                                                                                                                                            agtccGCCACCGGCTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GCTCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑CCCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGGTTCAGAAGCTGTAGGACGCGCAcc      >  7:5747/4‑149 (MQ=255)                                                                                                                                            agtccGCCACCGGCTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑CCCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTGTCAGAAGCTGTATGACGAGCAcc      >  7:11106/4‑149 (MQ=255)                                                                                                                                            aaCCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGACGAGCAcc      >  1:329667/1‑149 (MQ=255)                                                                                                                                             gtccGCCACCGGCTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGa              >  7:304507/3‑140 (MQ=255)                                                                                                                                             gtccGCCACCGGCTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGa              <  8:304507/138‑1 (MQ=255)                                                                                                                                             gtccGCCACCGGCTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGa              <  8:299180/138‑1 (MQ=255)                                                                                                                                             gtccGCCACCGGCTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGa              >  7:299180/3‑140 (MQ=255)                                                                                                                                             aCCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGACGAGCACCt     >  3:121312/1‑149 (MQ=255)                                                                                                                                                cGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGACGAGCACCTGAc  <  7:86186/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGCTCAGCCTGTCACCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAA‑ACGGCACTGGCGCA‑GATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCG‑ACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGC‑CGGCGTTTCAGAAGCTGTATGACGAGCACCTGAC  >  NC_002947/4149631‑4149922

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