Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,976,139	T→C	6.2%	Q211R (CAG→CGG)	flgJ ←	peptidoglycan hydrolase FlgJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,976,139	0	T	C	6.2%	181.1 / 2.4	64	Q211R (CAG→CGG)	flgJ	peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base C (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TTGCCCCAACCGGTTTCCAGCGCGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTT  >  NC_002947/4975992‑4976286                                                                                                                                                      |                                                                                                                                                    ttGCCCCAACCGGTTTCCAGCGCGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTg                                                                                                                                                    >  7:180538/1‑149 (MQ=255)   gcaccACCCGGTTTCCAGCGCGGCCGGGTCGACCAGGTTGCGCGGGTCTACACCAATTCGCCTGGGGGGCTGCTCT‑‑GCGATTGGCAGCATTGTGGCGACGAAGGCGTCGCTGTCGGTGTAAGCCTTGTTCGGCGCCAGCGGTGGCt                                                                                                                                                     <  8:94038/143‑1 (MQ=255)   gCCCCAACCGGTTTCCAGCGCGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTgg                                                                                                                                                   <  6:111933/148‑1 (MQ=255)     cccAACCGGTTTCCAGCGCGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGcc                                                                                                                                                <  4:45116/149‑1 (MQ=255)        aaCCGGTTTCCAGCGCGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCcacg                                                                                                                                             <  4:154965/149‑1 (MQ=255)                    gcgcGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCGGAGGCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGa                                                                                                                                   >  1:93238/1‑149 (MQ=255)                      gcgGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCcacgg                                                                                                                                            <  7:118899/136‑1 (MQ=255)                      gcgGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCcacgg                                                                                                                                            >  8:118899/1‑136 (MQ=255)                          ccTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATgcgc                                                                                                                           >  8:127914/1‑149 (MQ=255)                            tGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAg                                                                                                                         >  5:95988/1‑149 (MQ=255)                              ggCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAgac                                                                                                                       <  1:213932/149‑3 (MQ=255)                                ccaccaGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCg                                                                                                                     >  1:316826/1‑149 (MQ=255)                                ccaccaGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCg                                                                                                                     >  1:316828/1‑149 (MQ=255)                                ccaccaGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCg                                                                                                                     <  7:121268/149‑1 (MQ=255)                                  accaGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGa                                                                                                                                   <  7:124888/133‑1 (MQ=255)                                  accaGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGa                                                                                                                                   >  8:124888/1‑133 (MQ=255)                                           cgcgGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCa                                                                                                           <  2:54547/148‑1 (MQ=255)                                             cgGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACgg                                                                                                        <  5:221044/149‑1 (MQ=255)                                                   gtcacCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAgg                                                                                                  >  3:69788/3‑149 (MQ=255)                                                   gACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAgg                                                                                                  >  5:52615/1‑149 (MQ=255)                                                   gACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAg                                                                                                   <  6:67508/148‑1 (MQ=255)                                                    acacCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGc                                                                                                 >  4:123296/1‑149 (MQ=255)                                                      acCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCCGTGccgt                                                                                                                                              <  5:269093/102‑3 (MQ=255)                                                      acCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCCGTGccgt                                                                                                                                              <  5:264675/102‑3 (MQ=255)                                                      acCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCCGTGccgt                                                                                                                                              >  6:269093/1‑100 (MQ=255)                                                      acCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCCGTGccgt                                                                                                                                              >  6:264675/1‑100 (MQ=255)                                                       ccAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCtt                                                                                              <  4:69788/149‑1 (MQ=255)                                                        cAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAGGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTg                                                                                             >  7:80861/1‑149 (MQ=255)                                                         aaTGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTgg                                                                                            >  3:78052/1‑149 (MQ=255)                                                         aaTGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTgg                                                                                            >  3:78895/1‑149 (MQ=255)                                                          aTGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGc                                                                                           <  7:127914/149‑1 (MQ=255)                                                             cgcCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTa                                                                                         <  8:180538/148‑1 (MQ=255)                                                              gcCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACt                                                                                       <  6:113599/149‑1 (MQ=255)                                                               ccTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAGGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTg                                                                                      <  8:80861/149‑1 (MQ=255)                                                                       ccTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCg                                                                              <  2:316826/149‑1 (MQ=255)                                                                       ccTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCg                                                                              <  2:316828/149‑1 (MQ=255)                                                                        cgcgatg‑‑gCCATTGGCAGCATGGTGGCGACGAAGGCGGCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGt                                                                             <  4:268692/143‑1 (MQ=255)                                                                         tGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGt                                                                             >  4:54298/1‑148 (MQ=255)                                                                                   cATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTg                                                                                             <  5:180516/124‑1 (MQ=255)                                                                                   cATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTg                                                                                             >  6:180516/1‑124 (MQ=255)                                                                                   cATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggcag                                                                    >  8:92001/1‑149 (MQ=255)                                                                                        gcagcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGATGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcg                                                               >  1:266449/1‑149 (MQ=255)                                                                                        gcagcaTGGTGGCGACGAAGGCGTCGCTGGCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcaggcagt                                                               <  5:75749/149‑3 (MQ=255)                                                                                         cagcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcaggcagta                                                              <  4:78895/149‑4 (MQ=255)                                                                                         cagcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcaggcagta                                                              <  4:78052/149‑4 (MQ=255)                                                                                         cagcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcgg                                                              >  3:88121/1‑149 (MQ=255)                                                                                           gcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcggcg                                                            >  4:318707/1‑149 (MQ=255)                                                                                           gcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcggcg                                                            <  8:80251/149‑1 (MQ=255)                                                                                                   gCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGAACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTAtt                                                    <  2:93238/149‑1 (MQ=255)                                                                                                      acgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTgg                                                  >  3:216970/1‑148 (MQ=255)                                                                                                       cgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCa                                                <  3:54298/149‑1 (MQ=255)                                                                                                         aaGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGt                                              <  3:123296/149‑1 (MQ=255)                                                                                                            gCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTc                                           >  7:108333/1‑149 (MQ=255)                                                                                                                       ggTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCa                                >  5:125862/1‑149 (MQ=255)                                                                                                                          gAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGACCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGa                             >  3:98018/1‑149 (MQ=255)                                                                                                                           aaaGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGa                             >  3:186996/1‑148 (MQ=255)                                                                                                                             aGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGTCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGATCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATAc                          <  1:84214/149‑1 (MQ=255)                                                                                                                                cTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCgg                       >  3:247507/1‑149 (MQ=255)                                                                                                                                cTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCgg                       >  2:27889/1‑149 (MQ=255)                                                                                                                                     tCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGAc                                    <  7:57099/131‑1 (MQ=255)                                                                                                                                     tCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGAc                                    >  8:57099/1‑131 (MQ=255)                                                                                                                                       ggCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAg                >  8:122215/1‑149 (MQ=255)                                                                                                                                            cAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGTAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggt           <  2:266449/149‑1 (MQ=255)                                                                                                                                                 gTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCGGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTg      <  6:95988/149‑1 (MQ=255)                                                                                                                                                     cTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGtt  <  4:216970/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TTGCCCCAACCGGTTTCCAGCGCGGCCTGGGCCACCAGGTAGCGCGGGTCGACACCAATGCGCCTGGCGGCCTGCTCG‑‑GCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTT  >  NC_002947/4975992‑4976286

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