Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,906,463	A→T	6.1%	L112Q (CTG→CAG)	PP_4315 ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,906,463	0	A	T	6.1%	201.0 / 3.0	65	L112Q (CTG→CAG)	PP_4315	epimerase
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGG  >  NC_002947/4906316‑4906610                                                                                                                                                    |                                                                                                                                                     aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                     >  4:1607/1‑149 (MQ=255) aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                     >  6:228798/1‑149 (MQ=255) aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                     >  7:105616/1‑149 (MQ=255) aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                     >  6:154636/1‑149 (MQ=255) aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                     >  5:31576/1‑149 (MQ=255) aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGGTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                     >  3:162017/1‑149 (MQ=255)  aCGTTTCGCGCAGCATGGCGTCATAGATCACCACGATATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                    <  4:162017/149‑1 (MQ=255)    gTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTCGTTTGGCGGGGAAGTCCATGGCcagcag                                                                                                                                                  >  2:203683/1‑149 (MQ=255)     tttCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTg‑g                                                                                                                                                >  8:150790/1‑149 (MQ=255)       tCGAGCAGCATGGCGTCATCGATCACCTCGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑Gc                                                                                                                                               <  2:213235/148‑1 (MQ=255)         gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCAt                                                                                                                                            <  5:164431/149‑1 (MQ=255)          agcagcATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTG‑GCCAtg                                                                                                                                           <  7:150790/149‑2 (MQ=255)           gcagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCAt                                                                                                                                            >  5:81279/1‑147 (MQ=255)           gcagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCa                                                                                                                                          <  4:274602/149‑1 (MQ=255)           gcagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCa                                                                                                                                          <  3:263727/149‑1 (MQ=255)            cagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAg                                                                                                                                         <  3:1607/149‑1 (MQ=255)              gcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagcag                                                                                                                                                  >  6:201315/1‑139 (MQ=255)              gcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagcag                                                                                                                                                  <  5:201315/139‑1 (MQ=255)               caTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAgcgc                                                                                                                                      >  3:149678/1‑149 (MQ=255)                          cGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGAcg                                                                                                                           >  8:174776/1‑149 (MQ=255)                          cGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGAcg                                                                                                                           >  6:260241/1‑149 (MQ=255)                            aTCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTG‑GCCAtgg                                                                                                                                          <  5:209879/132‑3 (MQ=255)                             tCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGAcgcg                                                                                                                         >  3:72396/1‑148 (MQ=255)                             tCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACgcgc                                                                                                                        <  6:31576/149‑1 (MQ=255)                                ccacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGCGCCATCAGCGCCACGGCTGAcgcg                                                                                                                         <  7:45305/146‑1 (MQ=255)                                ccacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGAcgcg                                                                                                                         >  8:45305/1‑145 (MQ=255)                                 cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACt                                                                                                                    >  3:172564/1‑149 (MQ=255)                                   cGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCg                                                                                                                  >  4:200330/1‑149 (MQ=255)                                             gtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                         <  1:192105/148‑1 (MQ=255)                                             gtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTCGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                         <  1:203683/148‑1 (MQ=255)                                                 cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAtt                                                                                                    <  4:172564/149‑1 (MQ=255)                                                   aaGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                   >  8:118282/1‑148 (MQ=255)                                                    aGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaa                                                                                                 <  6:81279/149‑1 (MQ=255)                                                      gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagc                                                                                               <  8:105616/149‑1 (MQ=255)                                                       tACTCCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAACAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagc                                                                                               <  7:45160/148‑1 (MQ=255)                                                         cacCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaa                                                                                            <  5:154636/149‑1 (MQ=255)                                                         cacCGCCTGCGCCTGGCTCCAGCCCAGCGCCTGCCAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaa                                                                                            <  6:248561/149‑1 (MQ=255)                                                            cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                         >  7:244458/1‑133 (MQ=255)                                                            cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                         <  8:244458/133‑1 (MQ=255)                                                               cTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                           <  7:118766/128‑1 (MQ=255)                                                               cTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                           >  8:118766/1‑128 (MQ=255)                                                                    ccTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGcc                                                                                 >  4:117822/1‑149 (MQ=255)                                                                        gCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTg                                                                               <  4:149678/147‑1 (MQ=255)                                                                           cAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                          >  3:195352/1‑149 (MQ=255)                                                                             aGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAccc                                                                        <  4:72396/149‑1 (MQ=255)                                                                                  aGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTACGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTg                                                                   <  7:164350/149‑1 (MQ=255)                                                                                    cgcCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCt                                                                 >  6:179999/1‑149 (MQ=255)                                                                                     gcCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTc                                                                <  4:195352/149‑1 (MQ=255)                                                                                         gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGAGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGaa                                                             >  8:114680/1‑148 (MQ=255)                                                                                          cAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                            <  3:200330/148‑1 (MQ=255)                                                                                           aaaaGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAg                                                          >  5:255891/1‑149 (MQ=255)                                                                                            aaaGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGc                                                         <  2:107929/149‑1 (MQ=255)                                                                                            aaaGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGc                                                         <  3:117822/149‑1 (MQ=255)                                                                                              aGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCAc                                                       >  4:220549/1‑149 (MQ=255)                                                                                                gCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGt                                                     <  6:255891/149‑1 (MQ=255)                                                                                                  ccGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGt                                                     >  6:291426/1‑147 (MQ=255)                                                                                                    gCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGc                                                 <  5:228798/149‑1 (MQ=255)                                                                                                        gCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                             >  8:7859/1‑149 (MQ=255)                                                                                                             tCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGt                                                     <  3:101751/136‑1 (MQ=255)                                                                                                             tCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGt                                                     >  4:101751/1‑136 (MQ=255)                                                                                                               gACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAgg                                      <  5:260241/149‑1 (MQ=255)                                                                                                               gACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAgg                                      <  7:174776/149‑1 (MQ=255)                                                                                                                aCCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAgg                                      >  5:30680/1‑148 (MQ=255)                                                                                                                aCCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAgg                                      >  1:32202/1‑148 (MQ=255)                                                                                                                  cGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                            >  5:174063/1‑124 (MQ=255)                                                                                                                  cGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                            <  6:174063/124‑1 (MQ=255)                                                                                                                       cGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCAGAGGCCAAGGTAGCATgg                              >  7:139960/1‑149 (MQ=255)                                                                                                                        gggTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGc                             >  7:16584/1‑149 (MQ=255)                                                                                                                                      ggAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACt               <  7:65339/149‑1 (MQ=255)                                                                                                                                           tCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc          >  2:168310/1‑149 (MQ=255)                                                                                                                                           tCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGAATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCCAGGCCAAGGTAGCATTGCCGAACAGGTCGACTTCCAc          >  6:8363/1‑149 (MQ=255)                                                                                                                                            tgaTGGCCTGCTG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCcgcg         >  4:193608/3‑149 (MQ=255)                                                                                                                                              aTGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt       <  7:114680/149‑1 (MQ=255)                                                                                                                                                 gCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg     <  3:220549/148‑1 (MQ=255)                                                                                                                                                  ccAGCGG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTggg    <  3:247484/148‑1 (MQ=255)                                                                                                                                                   cagcag‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTggggg  >  5:142893/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAG‑GCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGG  >  NC_002947/4906316‑4906610

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