Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,050	+G	coding (277/666 nt)	yieH →	phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,046	1	.	G	98.7%	265.9 / ‑3.9	74	coding (273/666 nt)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (32/41);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑GGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896908‑3897177                                                                                                                                             |                                                                                                                                    tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑aa                                                                                                                                      <  1:978576/139‑1 (MQ=255)      aaGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑ggggg                                                                                                                                <  1:426148/139‑2 (MQ=255)            tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑ggggg                                                                                                                                <  1:378301/133‑2 (MQ=255)            tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑ggggg                                                                                                                                >  2:378301/1‑132 (MQ=255)            tAAAACTGCACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAg                                                                                                                           <  1:1005945/139‑1 (MQ=255)               aaCTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAgcgc                                                                                                                        <  2:1451046/139‑1 (MQ=255)                   gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCt                                                                                                                       >  1:452206/1‑136 (MQ=255)                   gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCt                                                                                                                       <  2:452206/136‑1 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAgc                                                                                                                          >  2:1491615/1‑132 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAgc                                                                                                                          <  1:1491615/132‑1 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  1:957003/139‑1 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  2:75646/139‑1 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  2:957003/1‑139 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  2:1175448/1‑139 (MQ=255)                          aTTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                                             <  2:1377551/139‑1 (MQ=255)                          aTTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                                             <  1:98933/139‑1 (MQ=255)                                     tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                                  >  2:98842/1‑139 (MQ=255)                                       tttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCaa                                                                                                >  2:1222351/1‑139 (MQ=255)                                        ttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               <  1:1222351/139‑1 (MQ=255)                                         tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAgcgc                                                                                                                        >  1:1454851/1‑113 (MQ=255)                                         tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAgcgc                                                                                                                        <  2:1454851/113‑1 (MQ=255)                                            cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    <  1:1445957/130‑1 (MQ=255)                                            cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    >  2:1445957/1‑130 (MQ=255)                                             ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  1:98842/139‑1 (MQ=255)                                             ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  1:179017/139‑1 (MQ=255)                                              tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgt                                                                                         >  1:886577/1‑139 (MQ=255)                                               gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    <  1:115647/127‑1 (MQ=255)                                               gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    >  2:115647/1‑127 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGTAG‑GGGGCTGGAGCGCTCCTg                                                                                                                   >  1:773329/1‑106 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑TAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  2:773329/106‑1 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGa                                                                                                                            >  1:1013989/1‑97 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGa                                                                                                                            <  2:1013989/97‑1 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                              >  2:1494886/1‑111 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                              <  1:1494886/111‑1 (MQ=255)                                                          tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                                             >  2:650839/1‑139 (MQ=255)                                                               tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 >  1:379694/1‑114 (MQ=255)                                                               tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  2:379694/114‑1 (MQ=255)                                                               tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaa                                                                        >  2:235956/1‑139 (MQ=255)                                                                aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCt                                                                                                                       <  2:1256814/91‑1 (MQ=255)                                                                aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCt                                                                                                                       >  1:1256814/1‑91 (MQ=255)                                                                aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTcc                                                                                                                     >  2:1571739/1‑93 (MQ=255)                                                                aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTcc                                                                                                                     <  1:1571739/93‑1 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  1:1366396/119‑1 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  2:1366396/1‑119 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTa                                                                                 >  1:153145/1‑128 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTa                                                                                 <  2:153145/128‑1 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACg                                                                              <  1:1419594/131‑1 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACg                                                                              >  2:1419594/1‑131 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                      >  1:129268/1‑139 (MQ=255)                                                                  cGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                            >  2:839274/1‑132 (MQ=255)                                                                  cGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                            <  1:839274/132‑1 (MQ=255)                                                                     aaaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                  <  1:650839/139‑1 (MQ=255)                                                                       aCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                <  2:1060443/139‑1 (MQ=255)                                                                                             gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                                      <  1:250864/95‑1 (MQ=255)                                                                                             gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                                      >  2:250864/1‑95 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                              <  1:743195/68‑1 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                              >  2:1016639/1‑68 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                              >  2:743195/1‑68 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                              <  1:1016639/68‑1 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              <  2:1653237/132‑1 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              >  1:1653237/1‑132 (MQ=255)                                                                                                cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                       <  1:317496/139‑1 (MQ=255)                                                                                                   aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                    <  1:235956/139‑1 (MQ=255)                                                                                                     gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑aa                                                                                                                                      >  1:745239/1‑39 (MQ=255)                                                                                                     gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑aa                                                                                                                                      <  2:745239/39‑1 (MQ=255)                                                                                                     gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTg                                                                                                                   <  1:879728/58‑1 (MQ=255)                                                                                                     gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTg                                                                                                                   >  2:879728/1‑58 (MQ=255)                                                                                                        gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCTATCACTGCGCCAATGTGTGTGGTa                                                                                     <  2:1146452/85‑1 (MQ=255)                                                                                                        gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTa                                                                                     >  1:1146452/1‑85 (MQ=255)                                                                                                        gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑ggg                                                                                                                                  <  2:170079/39‑1 (MQ=255)                                                                                                        gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑ggg                                                                                                                                  >  1:170079/1‑39 (MQ=255)                                                                                                             gCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:411837/75‑1 (MQ=255)                                                                                                             gCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:411837/1‑75 (MQ=255)                                                                                                             gCTGTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCacta                          >  1:1367447/1‑139 (MQ=255)                                                                                                                gTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt                       <  2:1367447/139‑1 (MQ=255)                                                                                                                gTTCGATTCAGAACTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt                       <  2:886577/139‑1 (MQ=255)                                                                                                                           aaCTGGAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAAc            >  2:629057/1‑139 (MQ=255)                                                                                                                                gAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                >  1:815785/1‑82 (MQ=255)                                                                                                                                gAAGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                <  2:815785/82‑1 (MQ=255)                                                                                                                                 aaGCCATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCa      <  1:1175448/139‑1 (MQ=255)                                                                                                                                     cATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  <  1:629057/139‑1 (MQ=255)                                                                                                                                     cATCG‑AAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCACCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  <  1:890869/139‑1 (MQ=255)                                                                                                                                             |                                                                                                                                    TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCG‑AA‑GGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896908‑3897177

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑