Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,054,925	A→G	6.7%	A118A (GCT→GCC)	glnG ←	fused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,054,925	0	A	G	6.7%	182.2 / 3.6	60	A118A (GCT→GCC)	glnG	fused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

ATCAGCACGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTT  >  NC_000913/4054782‑4055070                                                                                                                                                |                                                                                                                                                  aTCAGCACGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAgcg                                                                                                                                                <  8:192908/147‑1 (MQ=255)      cACGCTAATAGAAGAACGCGAAAGCCGACCGATAATCCGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCa                                                                                                                                          <  1:174260/148‑1 (MQ=255)         gCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAAc                                                                                                                                        >  3:283138/1‑147 (MQ=255)          cTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAAc                                                                                                                                        >  2:261770/1‑146 (MQ=255)             aTAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAgcg                                                                                                                                   >  3:21391/1‑148 (MQ=255)                gaagaaCGCGACAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGcc                                                                                                                                 <  7:327134/147‑1 (MQ=255)                     aCGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCt                                                                                                                           >  2:210987/1‑148 (MQ=255)                          aaaGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTtcgtc                                                                                                                      >  7:233066/1‑148 (MQ=255)                                 accgaTAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATc                                                                                                                >  1:75587/1‑147 (MQ=255)                                   cgaTAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCaaa                                                                                                             >  2:13013/1‑148 (MQ=255)                                      taataCGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACgg                                                                                                          >  6:35023/1‑148 (MQ=255)                                        ataCGGAACTCGTCCTGCATGGATGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATATCGCGTCTGCTGCTGTTCCTGGTAATGACTGATAGCGAGCTCAACCAGCTAAACTGATTCGGCGATATCAAACGGt                                                                                                         <  4:319625/147‑1 (MQ=255)                                         taCGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGttt                                                                                                       >  3:321372/1‑148 (MQ=255)                                               acacGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAg                                                                                                 <  6:312609/148‑1 (MQ=255)                                                    tCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAAt                                                                                            <  4:150764/148‑1 (MQ=255)                                                       tGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCa                                                                                          >  2:8860/1‑147 (MQ=255)                                                        gCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCaca                                                                                        >  8:201176/1‑146 (MQ=255)                                                         cATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCaa                                                                                         <  4:315750/146‑1 (MQ=255)                                                         cATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAACGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAAc                                                                                       >  4:228563/1‑148 (MQ=255)                                                          aTGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCAGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGTTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACg                                                                                      <  2:363012/148‑1 (MQ=255)                                                                    cTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCttgttg                                                                            <  2:142239/148‑1 (MQ=255)                                                                     ttCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGa                                                                           >  7:303249/1‑148 (MQ=255)                                                                       cGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGAt                                                                          <  1:8860/147‑1 (MQ=255)                                                                          cGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATGGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGgc                                                                      >  2:21382/1‑148 (MQ=255)                                                                          cGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATGGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGgc                                                                      >  2:15208/1‑148 (MQ=255)                                                                               atatCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTg                                                                  <  4:283138/147‑1 (MQ=255)                                                                               atatCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGa                                                                 >  7:41232/1‑148 (MQ=255)                                                                                  tCGGTCGTTGGGTCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACg                                                               <  5:143341/147‑1 (MQ=255)                                                                                    ggTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATAAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACgg                                                              >  1:305127/1‑146 (MQ=255)                                                                                       cGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGgcagca                                                         >  2:179392/1‑148 (MQ=255)                                                                                           gggCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATcc                                                      <  8:233066/147‑1 (MQ=255)                                                                                              ccGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGAt                                                  <  2:374778/148‑1 (MQ=255)                                                                                                  tAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAACCGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATc                                                       >  8:9906/1‑139 (MQ=255)                                                                                                   aaGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTACTGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGa                                              <  2:75658/147‑1 (MQ=255)                                                                                                   aaGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGa                                              <  8:303249/147‑1 (MQ=255)                                                                                                      cTGAACTTTACGCGGCTGCTGCTGTTCCTGGTTATGACTGATAGCGCGCTCAACCAGCGCCACTTCTTCGTCGATTTCATACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGGAGCATCCAGATCGGAAtgt                                          <  3:98029/148‑1 (MQ=255)                                                                                                      cTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAAtgt                                          <  5:35023/148‑1 (MQ=255)                                                                                                         aaCATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCg                                       >  6:49110/1‑148 (MQ=255)                                                                                                         aaCATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGGTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCg                                       <  7:178885/148‑1 (MQ=255)                                                                                                          aCATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCgg                                      <  1:179392/148‑1 (MQ=255)                                                                                                              tACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGa                                              >  2:73434/1‑136 (MQ=255)                                                                                                              tACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGa                                              <  1:73434/136‑1 (MQ=255)                                                                                                                cgcgGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATaa                                >  2:209284/1‑148 (MQ=255)                                                                                                                 gcgGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAat                               <  5:49110/148‑1 (MQ=255)                                                                                                                 gcgGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAat                               <  1:13013/148‑1 (MQ=255)                                                                                                                   ggCTGCTGCTGTTCCTGGTAATGACTGATGGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAatga                             <  1:15208/148‑1 (MQ=255)                                                                                                                   ggCTGCTGCTGTTCCTGGTAATGACTGATGGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAatga                             <  1:21382/148‑1 (MQ=255)                                                                                                                             gTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAg                   >  3:360018/1‑148 (MQ=255)                                                                                                                               tCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCa                 >  8:328055/1‑148 (MQ=255)                                                                                                                                   ggTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTg              >  5:82235/1‑147 (MQ=255)                                                                                                                                    gTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTgg             <  7:9906/147‑1 (MQ=255)                                                                                                                                     tAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGAt           >  3:116461/1‑148 (MQ=255)                                                                                                                                     tAACGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGAt           <  3:228563/148‑1 (MQ=255)                                                                                                                                         gACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGgcg       >  8:296491/1‑148 (MQ=255)                                                                                                                                          aCTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGgcgc      >  3:192955/1‑148 (MQ=255)                                                                                                                                             gATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATaa                                <  4:360921/119‑1 (MQ=255)                                                                                                                                             gATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATaa                                >  3:360921/1‑119 (MQ=255)                                                                                                                                              aTAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCACACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGATGAAGGAAGCATCCAGATCGGAATGTGCGGTCATACTGATGACCGGAAGCATTGGATGGCGCTg    <  5:280316/146‑1 (MQ=255)                                                                                                                                              aTAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGtt  >  4:127179/1‑148 (MQ=255)                                                                                                                                              aTAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGt   <  1:261770/147‑1 (MQ=255)                                                                                                                                                aGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAatg                              <  4:67318/118‑1 (MQ=255)                                                                                                                                                aGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAatg                              >  3:67318/1‑118 (MQ=255)                                                                                                                                                |                                                                                                                                                  ATCAGCACGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTT  >  NC_000913/4054782‑4055070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑