Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,527 Δ1 bp intergenic (+36/+23) yjgN → / ← valS hypothetical protein/valine tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,5270G.98.6% 314.2 / ‑3.6 71intergenic (+36/+23)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (45/25);  total (46/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCG  >  NC_000913/4478455‑4478596
                                                                        |                                                                     
ttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGtttt                                                                 <  1:5637783/76‑1 (MQ=255)
 tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtc                                                                <  2:5424553/76‑1 (MQ=255)
  cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                               >  1:173447/1‑76 (MQ=255)
  cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                               >  1:3126665/1‑76 (MQ=255)
   aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                              <  2:2930857/76‑1 (MQ=255)
    cGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                             <  2:2799004/76‑1 (MQ=255)
    cGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                             <  2:1214992/76‑1 (MQ=255)
     ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                                            >  2:445725/1‑76 (MQ=255)
     ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                                            >  2:110134/1‑76 (MQ=255)
         aGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAACGGCCGGAGCATGCTCCGGCCTTCGTTGTCATCACtgt                                                        >  1:5565557/1‑76 (MQ=255)
         aGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                        >  1:5536370/1‑76 (MQ=255)
          gTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                       >  2:186423/1‑76 (MQ=255)
           ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTgtgt                                                      >  2:2030042/1‑76 (MQ=255)
            tATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                     >  1:2352476/1‑76 (MQ=255)
              tGCCTTCTTTACCATTTTTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                   >  1:5651911/1‑76 (MQ=255)
              tGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                   >  1:5716301/1‑76 (MQ=255)
              tGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                   >  2:1919643/1‑76 (MQ=255)
              tGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                   <  1:140896/76‑1 (MQ=255)
                     tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                                            >  1:3581796/1‑76 (MQ=255)
                      ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                           <  2:988408/76‑1 (MQ=255)
                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                          <  2:23312/76‑1 (MQ=255)
                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                          <  2:4874901/76‑1 (MQ=255)
                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                          <  1:1436625/76‑1 (MQ=255)
                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                          >  2:1980049/1‑76 (MQ=255)
                        aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcg                                         >  1:1974792/1‑76 (MQ=255)
                        aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcg                                         <  2:4360065/76‑1 (MQ=255)
                         ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                        <  2:3543240/76‑1 (MQ=255)
                            tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                     >  1:3050028/1‑76 (MQ=255)
                             ctCTGTAAGAGATAAAAAAGGCCGGAGCAGGCTCCGGCCTTCGTTTTCATGACTGTTTTTTTATTACAGCgcggcg                                    >  1:2313164/2‑76 (MQ=255)
                              tCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                   <  2:4967062/76‑1 (MQ=255)
                               cTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGAt                                  >  1:4810318/1‑76 (MQ=255)
                                tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                 >  1:441915/1‑76 (MQ=255)
                                tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGGTTTGATTACAGCGCGGCGATa                                 >  2:2919129/1‑76 (MQ=255)
                                   aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                              >  2:3093254/1‑76 (MQ=255)
                                   aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                              >  2:1243197/1‑76 (MQ=255)
                                   aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                              >  1:2043116/1‑76 (MQ=255)
                                     gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTATAGCGCGGAGATAACAGc                            >  1:3377699/1‑76 (MQ=255)
                                     gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGc                            >  1:2207485/1‑76 (MQ=255)
                                        gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                          >  1:2974092/4‑76 (MQ=255)
                                          taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                        <  1:4821636/75‑1 (MQ=255)
                                          taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                        >  2:2548683/2‑76 (MQ=255)
                                           aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                       >  1:3385609/1‑76 (MQ=255)
                                           aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                       >  2:1553988/1‑76 (MQ=255)
                                            aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg                      <  2:5716299/76‑1 (MQ=255)
                                            aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg                      <  1:5155476/76‑1 (MQ=255)
                                              aaaGGCCGGAGCATGCTCGGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt                    <  1:266011/76‑1 (MQ=255)
                                              aaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt                    >  2:5207248/1‑76 (MQ=255)
                                               agggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                   <  1:487678/74‑1 (MQ=255)
                                               aaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                   >  1:728293/1‑76 (MQ=255)
                                                aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa                  <  2:5121424/76‑1 (MQ=255)
                                                 ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCaa                 <  1:799780/76‑1 (MQ=255)
                                                  gCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                >  2:3729522/1‑76 (MQ=255)
                                                      ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg             >  2:2223179/4‑76 (MQ=255)
                                                        agcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt           >  1:5125926/2‑76 (MQ=255)
                                                         gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt          >  2:3146269/1‑76 (MQ=255)
                                                         gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt          >  2:3669546/1‑76 (MQ=255)
                                                         gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt          >  1:1250280/1‑76 (MQ=255)
                                                         gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt          >  2:2694859/1‑76 (MQ=255)
                                                         gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt          <  2:4657852/76‑1 (MQ=255)
                                                         gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt          >  1:5638229/1‑76 (MQ=255)
                                                          catgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCCGtttc         >  2:85101/1‑76 (MQ=255)
                                                          catgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc         >  2:3477436/1‑76 (MQ=255)
                                                            tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc       >  2:2434130/1‑76 (MQ=255)
                                                            tgcTCCGGCCTGCGTTTTCATCACTGTGTTCTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc       >  2:2402950/1‑76 (MQ=255)
                                                             gcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgct      >  1:3313475/1‑76 (MQ=255)
                                                              cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt     <  1:5701769/76‑1 (MQ=255)
                                                              cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt     <  1:737922/76‑1 (MQ=255)
                                                               tCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt    >  2:3533811/1‑76 (MQ=255)
                                                                ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc   <  2:5651909/76‑1 (MQ=255)
                                                                ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc   >  1:291973/1‑76 (MQ=255)
                                                                 cGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcg  <  1:4394668/76‑1 (MQ=255)
                                                                        |                                                                     
TTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCG  >  NC_000913/4478455‑4478596

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