Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,223 A→T L58Q (CTG→CAG)  metJ ← DNA‑binding transcriptional repressor MetJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,2230AT84.5% 138.5 / 22.0 59L58Q (CTG→CAG) metJDNA‑binding transcriptional repressor MetJ
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base T (27/22);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

GCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAT  >  NC_000913/4127980‑4128447
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
gCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTc                                                                                                                                                                                                                           <  1:72334/251‑1 (MQ=255)
          gcGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGAGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGt                                                                                                                                                                                                                  <  1:218601/249‑1 (MQ=255)
               cAAATTCAATAAATGGCAATATTCCGTAGGCTCGAGAATCGTATCGCATCATGCGATTCCACTCCGCGTCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGTGATTTCGGCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAACCGCTTCGCACTGCAGCTCGCTGTGGGTAGCGTg                                                                                                                                                                                                            <  2:102895/251‑1 (MQ=255)
                 aaTTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAACGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGGGAc                                                                                                                                                                                                          >  2:192001/1‑251 (MQ=255)
                            ttGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGAGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCCTGCAAAAACGCTTCGCCCGGCAGCTCGCTGTTGGTAGCGTGAAGCAGGATGTTc                                                                                                                                                                                               >  2:443541/1‑251 (MQ=255)
                                aaTATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGTTTGTTCACCt                                                                                                                                                                                           >  1:207562/1‑251 (MQ=255)
                                            gccttattagcttAGCGCATCTGGCGATTCCCATCCGCGCCGCTCCTTTTTGCTGTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTGTGCCGCTTCCGGGATTTCGGCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgc                                                                                                                                                                               <  2:548195/239‑1 (MQ=255)
                                                catacgCGTAGCTCCTCAGGCGATTCCGCTTCGCGCCGCTCTTGTTTGCTTTATTATCACCACGGCTACGGGTTAATCCCCATCTCACGCATGAGCTCTTTAGCCGCTTCCTGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCAGCCGTCAAAGGTTGCCCGGTAAAGTCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGCCGACgcgtgcgt                                                                                                                                                                             <  2:154064/244‑1 (MQ=255)
                                                      ggTAGCGCATCAGGGGTTTCCACACCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATGTCGTCGCTGCGCTCTTTACGCATATCGTCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCgtg                                                                                                                                                                     <  2:501746/250‑3 (MQ=255)
                                                      cGTAGCGCAGCAGGCGATTCCACGCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGt                                                                                                                                                                     <  2:372142/251‑1 (MQ=255)
                                                            gcATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAt                                                                                                                                                               >  1:185230/1‑251 (MQ=255)
                                                            gcATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAt                                                                                                                                                               >  2:239376/1‑251 (MQ=255)
                                                               tCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                                                                                                                                            >  1:17559/1‑251 (MQ=255)
                                                                         ccACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCCCGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt                                                                                                                                                   <  1:94567/250‑1 (MQ=255)
                                                                           ccTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                                                                                                                                <  1:86890/250‑1 (MQ=255)
                                                                             tCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                                                                                                                                >  1:525542/1‑249 (MQ=255)
                                                                                  gccgcTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCTCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGTAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg                                                                                                                                         <  2:157186/251‑1 (MQ=255)
                                                                                                 ccactacgctcttccgagctgggTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGAGTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttt                                                                                                                          <  2:251768/231‑1 (MQ=255)
                                                                                                  ttAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttt                                                                                                                          <  1:16454/250‑1 (MQ=255)
                                                                                                      tATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGAc                                                                                                                     >  2:167129/1‑251 (MQ=255)
                                                                                                          ggattaaccccGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                                                                 >  1:336422/10‑251 (MQ=255)
                                                                                                            cACGTCTCCGGGTTACTCCCCATCTCCCGAATGATCTCTTTTGACGCATCCGGGATTTCGTCGCTGCGCTCTTTACGCCGATCAGCATCATCCGGCACAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTACGAGGTATGGCAACCGAAATCTTTTTGACTTGttt                                                                                                               >  2:73024/1‑250 (MQ=255)
                                                                                                                   ccGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCTCTTCCGGGATTTCGTCGCTTCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTATCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTTAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttc                                                                                                        <  2:351230/251‑1 (MQ=255)
                                                                                                                     gggTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttt                                                                                                                          >  1:251768/1‑231 (MQ=255)
                                                                                                                                    tcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGt                                                                                        >  2:35837/1‑250 (MQ=255)
                                                                                                                                     cACGCATTGTCTCTTTTTCCGCTTCCGTGATTTCGTCGCTGCTCGCTTTGCGCAGTTCGTCATCATCCGGCAAAGGTTGCCCGTTAAAGGCATGCTGAAACGCTTCGCACTGCATCTCGCTGTTGGTAGCGTGACGCTGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATTGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTAt                                                                                      <  1:361915/251‑1 (MQ=255)
                                                                                                                                        gCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg                                                                                   <  1:130346/251‑1 (MQ=255)
                                                                                                                                           tGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTTACGACGCGTGCGTTCATCGGTGAGGATATTTAACACCTTAAGAGGACTGGAAACGGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTACTCCGCGAATGGGCTg                                                                                >  2:306533/1‑251 (MQ=255)
                                                                                                                                               ctctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTCGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCCTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGACACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGTCTGatat                                                                            >  2:544531/1‑251 (MQ=255)
                                                                                                                                                tctTTTGCCGCTTCCGGGCTTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatat                                                                            >  1:184319/1‑250 (MQ=255)
                                                                                                                                                  ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                                                                         >  1:495716/1‑251 (MQ=255)
                                                                                                                                                       ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGACTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTAAGCGTATGGGCTGATATATTcgccg                                                                    >  1:281822/1‑251 (MQ=255)
                                                                                                                                                         gCTTCAGGGATTTCGACGCTGCGCTCTTTACGCAGATCGGAATCATCCGGCAAAGGTTGCCCGGGAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGAGGGTAGCGTGGCGCGGGTTGTTCACCTGCCGACGCTTGCGTTCATCGGTGCGGTTTTTTAACACCTTAAGAGGAATGGACACAGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGGGCTCAGCGAATGGGCTGATATATTCGACGCt                                                                  >  2:128337/1‑251 (MQ=255)
                                                                                                                                                                ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAg                                                                                           >  1:306379/1‑219 (MQ=255)
                                                                                                                                                                ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc                                                                                             >  1:517675/1‑217 (MQ=255)
                                                                                                                                                                ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc                                                                                             <  2:517675/217‑1 (MQ=255)
                                                                                                                                                                ggATTGCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAg                                                                                           <  2:306379/219‑1 (MQ=255)
                                                                                                                                                                    ttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                        >  1:312617/1‑250 (MQ=255)
                                                                                                                                                                    gtcgttgCTGCGCTCTTTACGCAGATCGGCATCCTCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGTTTCGCTCTGCAGCTCGCTGTTGGGAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACCCCTTAATAGGAATGGAAACCGTAAGCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCGGCGTATGGGCTGATATATTCGCCGCTCCATTCAGca                                                        <  2:358080/249‑2 (MQ=255)
                                                                                                                                                                             tcgcTCTTTACGCAGAGTGGAATCATCCGGCAAAGTTTTCCAGTAAAAGGCATGCAGCAACGCTTCGCAATGTAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCAGCGTTGAGGATTTTTAACGCCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                <  2:207562/248‑1 (MQ=255)
                                                                                                                                                                              cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAACTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                             <  1:177505/251‑1 (MQ=255)
                                                                                                                                                                                     tACGCAGATCGGCATCATCCGGCAAAGGGTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGACTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGAATGGGCTGATATATTAGCCGCTCCATTCCGCCATGAGATCCTTAATCctc                                      >  2:261850/1‑251 (MQ=255)
                                                                                                                                                                                       cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGCTATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCtt                                    >  2:132335/1‑251 (MQ=255)
                                                                                                                                                                                          tgATCGGCATCATCCGGCAACGGCTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGCTGGTAGCGTGACGCAGGTTGTTCACCTGCCGACGAGTGAGTTCATCGGTGAGGCTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGCTATATTCGCCGCTCCATTCAGCCATGCGATACTTAATCCTCTTCGt                                 >  2:32723/2‑251 (MQ=255)
                                                                                                                                                                                           gATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                <  1:337049/251‑1 (MQ=255)
                                                                                                                                                                                            aTCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                               <  1:440440/251‑1 (MQ=255)
                                                                                                                                                                                               gacatcatCCGGCAAAGGTTGCCCGGTAAAGGCCTGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGGTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACAGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCTTATGGGCTGATATATTCTCCGCTCCATTCAGCCATGAGCAACTTAATACTCTTCGTCAAt                             >  2:140916/3‑250 (MQ=255)
                                                                                                                                                                                                    catcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                       <  1:239376/251‑1 (MQ=255)
                                                                                                                                                                                                     atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGTCGCGTGCGTTCATCGGTGATGATTTTTAACACCTTAAGAgg                                                                                                                                             <  2:558881/132‑1 (MQ=255)
                                                                                                                                                                                                     atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAgg                                                                                                                                             >  1:558881/1‑132 (MQ=255)
                                                                                                                                                                                                     atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAGTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                       >  2:425150/1‑250 (MQ=255)
                                                                                                                                                                                                     atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCCTGAGATACTTAATCCTCTTCGTCAATAAATTg                       >  1:533526/1‑250 (MQ=255)
                                                                                                                                                                                                           cAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCCTGAGCTACTTAATCCTCTTCGTCAATAAATTGagacc                  >  2:304992/1‑249 (MQ=255)
                                                                                                                                                                                                            aaaGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                  >  1:553731/1‑248 (MQ=255)
                                                                                                                                                                                                               ggTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca            <  1:35837/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTAAAGGCATGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTACGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACATTTg    >  2:255398/1‑250 (MQ=255)
                                                                                                                                                                                                                         tAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt  >  1:220844/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCt                                                                                                                            <  2:438834/121‑1 (MQ=255)
                                                                                                                                                                                                                                 tGCAGAAACGCTTCGCACTGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCt                                                                                                                            >  1:438834/1‑121 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
GCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAT  >  NC_000913/4127980‑4128447

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