Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,959,628 A→G 5.3% N73N (AAT→AAC yecM ← putative metal‑binding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,959,6280AG 5.3% 235.3 / 3.0 75N73N (AAT→AACyecMputative metal‑binding enzyme
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (2/2);  total (37/38)
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.08e-01

GATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGGGAAAGTT  >  NC_000913/1959483‑1959770
                                                                                                                                                 |                                                                                                                                              
gATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCAtt                                                                                                                                              >  1:423885/1‑148 (MQ=255)
  tCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGa                                                                                                                                            <  5:421830/148‑1 (MQ=255)
      ccGGCAGAACGATTTCAATATGTTCCCAACATTCGTGTGGGTAAAGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCat                                                                                                                                        <  6:265987/148‑1 (MQ=255)
           aGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATAtttt                                                                                                                                   <  1:250857/148‑1 (MQ=255)
             aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCGTTGATCATATTTTc                                                                                                                                  >  2:214358/1‑147 (MQ=255)
             aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCGTTGATCATATTTTc                                                                                                                                  >  2:208292/1‑147 (MQ=255)
             aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCGTTGATCATATTTTc                                                                                                                                  <  1:208292/147‑1 (MQ=255)
             aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTCCCGTTGATCATATTTTc                                                                                                                                  <  1:214358/147‑1 (MQ=255)
                 aTTTCAATATGTTCCCAACCTTCGTGTGGGTCACGTTTTTCCCCGCGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCAAGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTATCGGAc                                                                                                                              >  1:138397/1‑147 (MQ=255)
                    tCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACaaa                                                                                                                           <  5:324347/147‑1 (MQ=255)
                     cAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAg                                                                                                                         >  3:70998/1‑148 (MQ=255)
                       atatGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCt                                                                                                                       >  8:238833/1‑148 (MQ=255)
                            ttCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCca                                                                                                                  >  8:182540/1‑148 (MQ=255)
                                cAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACt                                                                                                              >  5:395449/1‑148 (MQ=255)
                                 aaCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTg                                                                                                             >  1:77755/1‑148 (MQ=255)
                                             ggTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCcgcg                                                                                                 <  2:196451/148‑1 (MQ=255)
                                               tAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTg                                                                                                             >  1:207466/1‑134 (MQ=255)
                                               tAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTg                                                                                                             <  2:207466/134‑1 (MQ=255)
                                                    tttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAg                                                                                          >  1:10148/1‑148 (MQ=255)
                                                     ttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGc                                                                                         >  8:40701/1‑148 (MQ=255)
                                                             gggCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagca                                                                                 <  6:395449/148‑1 (MQ=255)
                                                                  aCGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGc                                                                            >  2:263890/1‑148 (MQ=255)
                                                                  aCGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGc                                                                            >  3:78010/1‑148 (MQ=255)
                                                                       aGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtt                                                                         >  2:370433/1‑146 (MQ=255)
                                                                       aGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtt                                                                         <  1:370433/146‑1 (MQ=255)
                                                                        gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                                      <  4:78010/148‑1 (MQ=255)
                                                                         cTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGa                                                                     <  7:238833/148‑1 (MQ=255)
                                                                            cacaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGAt                                                                    <  7:182540/146‑1 (MQ=255)
                                                                            cacaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATgg                                                                  >  6:259401/1‑148 (MQ=255)
                                                                            cacaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATg                                                                   <  2:347576/147‑1 (MQ=255)
                                                                             acaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGt                                                                          <  3:188571/139‑1 (MQ=255)
                                                                             acaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGt                                                                          >  4:188571/1‑139 (MQ=255)
                                                                             acaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAAAAGCCAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATgg                                                                  >  4:76300/1‑147 (MQ=255)
                                                                               aaTAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                                  >  8:98012/1‑129 (MQ=255)
                                                                               aaTAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                                  <  7:98012/129‑1 (MQ=255)
                                                                                    aaaaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCa                                                           >  3:90150/1‑147 (MQ=255)
                                                                                      aaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCa                                                           >  1:191459/1‑145 (MQ=255)
                                                                                       aCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAg                                                        >  7:437754/1‑147 (MQ=255)
                                                                                        cTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                                      <  5:333555/148‑1 (MQ=255)
                                                                                           ccAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaaa                                                     <  2:10148/146‑1 (MQ=255)
                                                                                               tGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCa                                                                                           <  3:355282/104‑1 (MQ=255)
                                                                                               tGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCa                                                                                           >  4:355282/1‑104 (MQ=255)
                                                                                               tGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGAt                                               >  5:254335/1‑148 (MQ=255)
                                                                                                gcgcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                                  >  1:154005/1‑112 (MQ=255)
                                                                                                gcgcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                                  <  2:154005/112‑1 (MQ=255)
                                                                                                 cgcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCg                                             <  8:437754/148‑1 (MQ=255)
                                                                                                 cgcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCg                                             <  6:284263/148‑1 (MQ=255)
                                                                                                  gcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTATTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCa                                                                >  4:254336/1‑128 (MQ=255)
                                                                                                  gcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTATTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCa                                                                <  3:254336/128‑1 (MQ=255)
                                                                                                   cAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGc                                           <  3:291047/148‑1 (MQ=255)
                                                                                                   cAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGc                                           <  2:423885/148‑1 (MQ=255)
                                                                                                    aaCCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTgg                                                                             >  6:38792/1‑113 (MQ=255)
                                                                                                    aaCCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTgg                                                                             <  5:38792/113‑1 (MQ=255)
                                                                                                         gTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAg                                     >  8:204133/1‑148 (MQ=255)
                                                                                                          tACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCa                                      <  2:77755/146‑1 (MQ=255)
                                                                                                           aCCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCg                                   >  4:381213/1‑148 (MQ=255)
                                                                                                            ccGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATTTTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCgg                                  <  3:356089/148‑1 (MQ=255)
                                                                                                                ttCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGt                               >  7:138590/1‑147 (MQ=255)
                                                                                                                ttCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTg                              <  5:259401/148‑1 (MQ=255)
                                                                                                                  cATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGa                             >  7:59225/1‑147 (MQ=255)
                                                                                                                   aTGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGata                           <  1:263890/148‑1 (MQ=255)
                                                                                                                   aTGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGat                            <  3:76300/147‑1 (MQ=255)
                                                                                                                     gCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCaa                                                               <  3:133391/110‑1 (MQ=255)
                                                                                                                     gCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGGTTTGATGGCaa                                                               >  4:133391/1‑110 (MQ=255)
                                                                                                                       aaTTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCa                       <  8:138590/148‑1 (MQ=255)
                                                                                                                        aTTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCaa                      >  6:240240/1‑148 (MQ=255)
                                                                                                                           tAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAccc                   >  7:307932/1‑148 (MQ=255)
                                                                                                                              aCAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAg                <  7:40701/148‑1 (MQ=255)
                                                                                                                                   cAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGAcggc           >  8:246570/1‑148 (MQ=255)
                                                                                                                                       tCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                                      <  3:235706/101‑1 (MQ=255)
                                                                                                                                       tCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                                      >  4:235706/1‑101 (MQ=255)
                                                                                                                                              gCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGa                             <  5:71850/119‑1 (MQ=255)
                                                                                                                                              gCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGa                             >  6:71850/1‑119 (MQ=255)
                                                                                                                                              gCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTGGATGGCAACGCAAAGAAATGTGAGCGGCAGTCAGCgg                                  <  1:139455/114‑1 (MQ=255)
                                                                                                                                               ccATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGGGAAAGtt  <  8:59225/145‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGGGAAAGTT  >  NC_000913/1959483‑1959770

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: