Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,511 Δ1 bp intergenic (+20/+39) 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,5110T.98.5% 298.1 / ‑2.1 68intergenic (+20/+39)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (42/25);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAT  >  NC_000913/4478440‑4478582
                                                                       |                                                                       
cttcttCGCCAGTATTTCACGGGGAGGGACTCCTACTTTACCATTTCTGTACGAGATAAAAAAGGCCGgagcatgc                                                                   <  1:2821555/76‑1 (MQ=255)
cttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatgc                                                                   <  1:2839571/76‑1 (MQ=255)
 ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCt                                                                  >  1:439378/1‑76 (MQ=255)
 ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCt                                                                  >  2:535842/1‑76 (MQ=255)
  tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTc                                                                 >  2:4821635/1‑76 (MQ=255)
    ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                               <  2:4646920/76‑1 (MQ=255)
      aGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGc                                                             >  1:5121425/1‑76 (MQ=255)
         cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                          <  1:2565569/76‑1 (MQ=255)
         cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                          >  2:3546508/1‑76 (MQ=255)
         cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                          <  1:3496131/76‑1 (MQ=255)
          aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                         <  1:5280702/76‑2 (MQ=255)
          aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                         >  2:5280701/1‑75 (MQ=255)
          aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                         >  2:137594/1‑75 (MQ=255)
          aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGTGCATGCTCCGGCCttc                                                         <  1:1613704/76‑2 (MQ=255)
           gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                        >  1:2613792/1‑74 (MQ=255)
           gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                        >  1:4238040/1‑74 (MQ=255)
           gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                        >  2:4199246/1‑74 (MQ=255)
           gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                        >  1:3323015/1‑74 (MQ=255)
            tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgt                                                       <  1:25015/76‑4 (MQ=255)
               ttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGtttt                                                   <  1:5637783/76‑1 (MQ=255)
                tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtc                                                  <  2:5424553/76‑1 (MQ=255)
                 cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                 >  1:3126665/1‑76 (MQ=255)
                 cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                 >  1:173447/1‑76 (MQ=255)
                  aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                <  2:2930857/76‑1 (MQ=255)
                   cGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                               <  2:1214992/76‑1 (MQ=255)
                   cGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                               <  2:2799004/76‑1 (MQ=255)
                    ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                              >  2:110134/1‑76 (MQ=255)
                    ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                              >  2:445725/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                                     >  2:1919643/1‑76 (MQ=255)
                             tGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                     >  1:5716301/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                            <  1:1436625/76‑1 (MQ=255)
                                      tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                            <  2:23312/76‑1 (MQ=255)
                                      tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                            <  2:4874901/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                >  1:2043116/1‑76 (MQ=255)
                                                  aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                >  2:3093254/1‑76 (MQ=255)
                                                  aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                >  2:1243197/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          >  2:2548683/2‑76 (MQ=255)
                                                         taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc          <  1:4821636/75‑1 (MQ=255)
                                                          aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct         >  1:3385609/1‑76 (MQ=255)
                                                          aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct         >  2:1553988/1‑76 (MQ=255)
                                                           aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg        <  1:5155476/76‑1 (MQ=255)
                                                           aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg        <  2:5716299/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)
                                                                       |                                                                       
CTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAT  >  NC_000913/4478440‑4478582

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: