Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,527,044	T→C	6.1%	L186P (CTG→CCG)	PP_4860 →	mannose‑6‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,527,044	0	T	C	6.1%	193.7 / 2.7	65	L186P (CTG→CCG)	PP_4860	mannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (32/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 = 8.17e-01

TACTGGGGCAAAGTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCGGCCCATGGCCTGATGCTGGAAAA  >  NC_002947/5526902‑5527191                                                                                                                                               |                                                                                                                                                     tactGGGGCAAAGTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGa                                                                                                                                                 >  6:213729/1‑148 (MQ=255)  actGGGGCAAAGTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATAATCGACCAGCAGTTCGCCCGTGACCATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTTCAGAAt                                                                                                                                               >  1:247959/1‑149 (MQ=255)     ggggCAAAGTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCg                                                                                                                                            >  2:291023/1‑149 (MQ=255)       ggCAAAGTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCg                                                                                                                                            >  5:103018/1‑147 (MQ=255)             gTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCa                                                                                                                                   <  1:291023/149‑1 (MQ=255)                cgcgAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCAcc                                                                                                                                 <  5:213729/148‑1 (MQ=255)                   gAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGc                                                                                                                             <  5:168125/149‑1 (MQ=255)                      aaCCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGa                                                                                                                          <  2:343657/149‑1 (MQ=255)                       aCCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGaa                                                                                                                         <  4:372752/149‑1 (MQ=255)                           ggtggAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGt                                                                                                                     >  2:41459/1‑149 (MQ=255)                            gtggAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGt                                                                                                                     <  6:103018/148‑1 (MQ=255)                                     aCGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTCGACTAAGCGTTCCTGCAgg                                                                                                            >  8:88259/1‑148 (MQ=255)                                            aTGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAg                                                                                                             <  7:119464/140‑1 (MQ=255)                                            aTGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAg                                                                                                             >  8:119464/1‑140 (MQ=255)                                                cgcgCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGc                                                                                                <  1:326373/149‑1 (MQ=255)                                                  cgcTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGt                                                                                              <  1:41459/149‑1 (MQ=255)                                                    cTGGACATCATCGACCAGCAGTGCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGACGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCTTTCCTGCAGGTACTGGCTGTGCGTga                                                                                            <  2:247959/149‑1 (MQ=255)                                                       gACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTgacga                                                                                         >  7:232792/1‑149 (MQ=255)                                                          atcatcGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCCGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTgacgac                                                                                        >  7:34348/1‑147 (MQ=255)                                                          atcatcGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCCGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTgacgac                                                                                        <  8:34348/147‑1 (MQ=255)                                                            catcGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCa                                                                                                                                                   <  2:304670/87‑1 (MQ=255)                                                            catcGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCa                                                                                                                                                   >  1:304670/1‑87 (MQ=255)                                                              tcGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACAt                                                                                  >  8:135983/1‑149 (MQ=255)                                                               cGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTc                                                                                                                   <  6:163308/115‑1 (MQ=255)                                                               cGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTc                                                                                                                   >  5:163308/1‑115 (MQ=255)                                                                 aCCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCCGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTgacga                                                                                         <  8:110955/139‑1 (MQ=255)                                                                 aCCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCCGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTgacga                                                                                         >  7:110955/1‑139 (MQ=255)                                                                  ccAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTcc                                                                              <  6:49696/149‑1 (MQ=255)                                                                   cagcagTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCg                                                                             >  5:110135/1‑149 (MQ=255)                                                                    agcagTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGc                                                                            >  8:141117/1‑149 (MQ=255)                                                                      cagTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTga                                                                                            >  3:368813/1‑131 (MQ=255)                                                                      cagTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTga                                                                                            <  4:368813/131‑1 (MQ=255)                                                                              ccGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCgctgc                                                                  >  1:188821/1‑149 (MQ=255)                                                                              ccGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAC‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCt                                                                                                    >  7:160094/1‑115 (MQ=255)                                                                              ccGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAC‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCt                                                                                                    <  8:160094/115‑1 (MQ=255)                                                                               cGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTgg                                                                                                      <  2:94493/112‑1 (MQ=255)                                                                               cGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTgg                                                                                                      >  1:94493/1‑112 (MQ=255)                                                                                   aCGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTa                                                                                                          <  3:243240/104‑1 (MQ=255)                                                                                   aCGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTa                                                                                                          >  4:243240/1‑104 (MQ=255)                                                                                    cGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAAt                                                            <  7:135983/149‑1 (MQ=255)                                                                                     gATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGATCGCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTga                                                                                            <  1:264332/117‑1 (MQ=255)                                                                                     gATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTga                                                                                            >  2:264332/1‑116 (MQ=255)                                                                                          ccTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCtgt                                                                                                  >  8:322812/1‑105 (MQ=255)                                                                                          ccTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCtgt                                                                                                  <  7:322812/105‑1 (MQ=255)                                                                                            tGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAg                                                                         >  6:152465/1‑128 (MQ=255)                                                                                              tACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGTCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAgg                                                   >  4:332929/1‑148 (MQ=255)                                                                                              tACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGc                                                  >  5:325033/1‑149 (MQ=255)                                                                                                cGAAGCCTGCCTGGGCGAAGACTGTGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCa                                                 <  6:110135/148‑1 (MQ=255)                                                                                                cGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGATGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCAc                                                >  1:93739/1‑149 (MQ=255)                                                                                                cGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCa                                                 <  3:332929/148‑1 (MQ=255)                                                                                                  aaGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTc                                              >  4:134484/1‑149 (MQ=255)                                                                                                     cctgcctgGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAg                                           >  2:371743/1‑149 (MQ=255)                                                                                                      ctgcctgGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAgg                                          <  7:141117/149‑1 (MQ=255)                                                                                                           tgGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCAt                                     >  4:267158/1‑149 (MQ=255)                                                                                                             ggCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCAtt                                    <  1:39726/148‑1 (MQ=255)                                                                                                                gAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCa                                                                          <  8:162054/107‑1 (MQ=255)                                                                                                                 aaGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCAtt                               <  3:134484/149‑1 (MQ=255)                                                                                                                   gACTGGGCCGACCTCGGCAGTGGCCCGCTGAAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGGGTGACGACGAACATGTCCGGCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTg                              >  7:88380/1‑148 (MQ=255)                                                                                                                      tGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAAcc                          >  8:236668/1‑149 (MQ=255)                                                                                                                         gCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCTGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCgg                        <  7:88259/148‑1 (MQ=255)                                                                                                                             aCCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCAttt                                   <  7:133418/133‑1 (MQ=255)                                                                                                                                tCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCGGCCCATGGc                >  1:67085/1‑149 (MQ=255)                                                                                                                                       tGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCGGCCCATGGCCTGATGc         <  6:325033/149‑1 (MQ=255)                                                                                                                                         gCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCGGCCCATGGCCTGATGCTg       <  2:188821/149‑1 (MQ=255)                                                                                                                                              cTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCGGCCCATGGCCTGATGCTGGaaaa  <  1:371743/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     TACTGGGGCAAAGTGCGCGAAAACCTGGTGGAGTCCACGTTGAATGCCGCGCTGGACATCATCGACCAGCAGTTCGCCCGTGACGATGGCCTGTACGAAGCCTGCCTGGGCGAAGACTGGGCCGACCTCGGCAGTGGCCCGCTGCAGAATCCGCAGAT‑GCACCTGGCCGAAGCGTTCCTGCAGGTACTGGCTGTGCGTGACGACGAACATGTCCGCCAGTCGCTGCTGCAATTGTGTCAGGCACTCGAGGCGCATTTCATTGAACCGGCCCATGGCCTGATGCTGGAAAA  >  NC_002947/5526902‑5527191

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