Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,563,154	C→T	6.0%	T381M (ACG→ATG)	glpD →	sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,563,154	0	C	T	6.0%	206.4 / 2.4	66	T381M (ACG→ATG)	glpD	sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (2/2);  total (31/35)
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.32e-01

ATTACACCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCGG  >  NC_000913/3563010‑3563301                                                                                                                                                 |                                                                                                                                                    aTTACACCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGaaa                                                                                                                                                 <  8:322373/149‑1 (MQ=255) aTTACACCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGa                                                                                                                                                   >  2:40105/1‑147 (MQ=255)     cacCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGTTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAagaga                                                                                                                                             <  1:363103/149‑1 (MQ=255)     cacCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGATGaaa                                                                                                                                                 >  1:91755/1‑145 (MQ=255)     cacCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGATGAAaga                                                                                                                                               >  1:86620/1‑147 (MQ=255)            gATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGATGAAAGAGAGTGTGc                                                                                                                                       <  2:91755/148‑1 (MQ=255)            gATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGATGAAAGAGAGTGTGc                                                                                                                                       <  2:86620/148‑1 (MQ=255)             atatTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTa                                                                                                                                     >  5:108587/1‑149 (MQ=255)                  cATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGAGCACCTACCGAAAACTGTAGAAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATAGGCCGGGCAAGCACGAAAGCGAGTGTGCTACCggg                                                                                                                                >  8:220913/1‑149 (MQ=255)                    tgatgaAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTg                                                                                                                              >  7:73595/1‑149 (MQ=255)                          aaaTGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCa                                                                                                                        <  8:73595/149‑1 (MQ=255)                           aaTGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCAt                                                                                                                       <  5:368861/149‑1 (MQ=255)                                cAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAg                                                                                                                  <  3:333549/149‑1 (MQ=255)                                cAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAg                                                                                                                  <  6:194261/149‑1 (MQ=255)                                cAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAg                                                                                                                  <  2:10506/149‑1 (MQ=255)                                 aaaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAagaga                                                                                                                                             >  2:293757/1‑121 (MQ=255)                                 aaaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAagaga                                                                                                                                             <  1:293757/121‑1 (MQ=255)                                 aaaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAgg                                                                                                                 >  2:253759/1‑149 (MQ=255)                                 aaaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAgg                                                                                                                 >  6:283687/1‑149 (MQ=255)                                 aaaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAgg                                                                                                                 >  6:288896/1‑149 (MQ=255)                                  aaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGc                                                                                                                <  5:288896/149‑1 (MQ=255)                                  aaGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGc                                                                                                                <  5:283687/149‑1 (MQ=255)                                      aCCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGAcc                                                                                                            >  5:12715/1‑149 (MQ=255)                                            gctgTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGAACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGcgac                                                                                                       >  4:149798/1‑148 (MQ=255)                                               gTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAgt                                                                                                                                           >  1:133572/1‑109 (MQ=255)                                               gTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAgt                                                                                                                                           <  2:133572/109‑1 (MQ=255)                                                tCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGAcc                                                                                                            >  5:107418/1‑139 (MQ=255)                                                tCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGAcc                                                                                                            <  6:107418/139‑1 (MQ=255)                                                tCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGCACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGAt                                                                                                    <  3:149798/147‑1 (MQ=255)                                                  ggTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGAGGATTATg                                                                                                >  2:110646/1‑149 (MQ=255)                                                  ggTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTAt                                                                                                 <  6:245698/148‑1 (MQ=255)                                                   gTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGc                                                                                                                                       >  8:378435/1‑109 (MQ=255)                                                    tATTCGGCGGTAAGCTGACAACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGcc                                                                                              >  6:262903/1‑149 (MQ=255)                                                      ttCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACcg                                                                                                           >  6:163659/1‑134 (MQ=255)                                                      ttCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACcg                                                                                                           <  5:163659/134‑1 (MQ=255)                                                       tCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAgtgt                                                                                                                                         <  8:291986/103‑1 (MQ=255)                                                       tCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAgtgt                                                                                                                                         >  7:291986/1‑103 (MQ=255)                                                         ggcggTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTc                                                                                          <  5:34795/148‑1 (MQ=255)                                                            ggTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCt                                                                                      >  3:330063/1‑149 (MQ=255)                                                            ggTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCt                                                                                      <  2:353328/149‑1 (MQ=255)                                                                  cTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGcgccg                                                                                >  3:148915/1‑149 (MQ=255)                                                                         ccTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTAt                                                                         <  1:253759/149‑1 (MQ=255)                                                                         ccTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTAt                                                                         <  5:49868/149‑1 (MQ=255)                                                                             ccGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGt                                                                     >  4:207942/1‑149 (MQ=255)                                                                               gAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTc                                                                   >  6:230427/1‑149 (MQ=255)                                                                                aaaaCTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTcc                                                                  <  4:148915/149‑1 (MQ=255)                                                                                     tGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCctgact                                                             >  5:170810/1‑149 (MQ=255)                                                                                         ggAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAAt                                                         >  1:303454/1‑149 (MQ=255)                                                                                          gAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATc                                                        <  6:108587/149‑1 (MQ=255)                                                                                                      gAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTgg                                                   >  1:209189/1‑142 (MQ=255)                                                                                                      gAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTgg                                                   <  2:209189/142‑1 (MQ=255)                                                                                                         aaaCTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGc                                                                                        <  6:239113/102‑1 (MQ=255)                                                                                                         aaaCTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGc                                                                                        >  5:239113/1‑102 (MQ=255)                                                                                                               aCGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTc                                    <  4:330063/148‑1 (MQ=255)                                                                                                                   cGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACt                               <  5:262903/149‑1 (MQ=255)                                                                                                                          tCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACgg                          >  5:235596/1‑147 (MQ=255)                                                                                                                            aGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGAtt                                                                                                   <  3:181921/72‑1 (MQ=255)                                                                                                                            aGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGAtt                                                                                                   >  4:181921/1‑72 (MQ=255)                                                                                                                            aGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGgcag                       >  4:296265/1‑148 (MQ=255)                                                                                                                                   tGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGAGGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGc                <  1:110646/148‑1 (MQ=255)                                                                                                                                     gCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGa              <  6:12715/148‑1 (MQ=255)                                                                                                                                               gACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAgctgc         >  1:289436/1‑143 (MQ=255)                                                                                                                                               gACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAgctgc         <  2:289436/143‑1 (MQ=255)                                                                                                                                               gACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCg   <  2:303454/149‑1 (MQ=255)                                                                                                                                               gACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCg   <  8:72348/149‑1 (MQ=255)                                                                                                                                                aCGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCgg  <  6:235596/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    ATTACACCCTTGATATTCATGATGAAAATGGCAAAGCACCGCTGCTGTCGGTATTCGGCGGTAAGCTGACCACCTACCGAAAACTGGCGGAACATGCGCTGGAAAAACTAACGCCGTATTATCAGGGTATTGGCCCGGCATGGACGAAAGAGAGTGTGCTACCGGGTGGCGCCATTGAAGGCGACCGCGACGATTATGCCGCTCGCCTGCGCCGCCGCTATCCGTTCCTGACTGAATCGCTGGCGCGTCATTACGCTCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCGG  >  NC_000913/3563010‑3563301

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