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.7% 344.6 / ‑3.8 77intergenic (+20/+39)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (41/35);  minor base C (1/0);  total (42/35)
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

GCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAT  >  NC_000913/4478439‑4478582
                                                                        |                                                                       
gCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatg                                                                    >  2:1501711/1‑76 (MQ=255)
 cttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatgc                                                                   >  1:1686183/1‑76 (MQ=255)
  gtcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCt                                                                  <  2:1490776/75‑1 (MQ=255)
   tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTc                                                                 <  1:3603145/76‑1 (MQ=255)
    cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTcc                                                                <  1:4085748/76‑1 (MQ=255)
    cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTcc                                                                <  2:444812/76‑1 (MQ=255)
     ttAGGCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                               <  2:4094229/76‑1 (MQ=255)
        gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGcc                                                            >  1:125426/1‑76 (MQ=255)
         ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCt                                                           <  1:3161219/76‑1 (MQ=255)
          cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                          >  1:1402123/1‑76 (MQ=255)
           aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                         >  1:4238140/1‑75 (MQ=255)
           aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                         <  1:4726168/76‑2 (MQ=255)
            gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                        >  1:367421/1‑74 (MQ=255)
             tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgt                                                       <  2:4238140/76‑4 (MQ=255)
               tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGttt                                                    >  1:707463/1‑76 (MQ=255)
                ttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGtttt                                                   >  1:2343422/1‑76 (MQ=255)
                 tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtc                                                  >  2:2033553/1‑76 (MQ=255)
                  cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                 <  2:3636799/76‑1 (MQ=255)
                   aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAACGCCGGCGCATGCTCCGGCCTTCGTTTtcac                                                >  2:1519463/1‑75 (MQ=255)
                    cGGGGAGTGATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                               <  1:4388812/76‑1 (MQ=255)
                     ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatcc                                              >  2:56451/1‑75 (MQ=255)
                     ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                              <  2:985894/76‑1 (MQ=255)
                      gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCAc                                             <  1:1948184/76‑1 (MQ=255)
                       ggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACt                                            <  1:548008/76‑1 (MQ=255)
                        gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                           <  1:4321047/76‑1 (MQ=255)
                        gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                           <  1:1256798/76‑1 (MQ=255)
                        gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                           >  2:4741804/1‑76 (MQ=255)
                         aGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                          >  1:842519/1‑76 (MQ=255)
                          gTTAGGCCTTCTTTACCATTTCTGTAAGAGATAAAGACGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                         <  1:3970872/76‑1 (MQ=255)
                           ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTgtgt                                        >  2:2329213/1‑76 (MQ=255)
                            tATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                       <  2:2912015/76‑1 (MQ=255)
                            tATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                       >  2:422308/1‑76 (MQ=255)
                             aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGttt                                      <  1:2631366/76‑1 (MQ=255)
                               gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTg                                    <  2:4405138/76‑1 (MQ=255)
                               gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTg                                    >  2:3120087/1‑76 (MQ=255)
                                ccTTCTTTCCCATTTCTGTAAGCGCTAAAAACGGCCGGAGCATGCTCCGGCATTCGTTTTCATCACTGTGtttttt                                   >  2:1315322/1‑74 (MQ=255)
                                ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGa                                   >  2:4342939/1‑76 (MQ=255)
                                 cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAt                                  <  2:4163031/76‑1 (MQ=255)
                                 cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAt                                  >  2:2234281/1‑76 (MQ=255)
                                  ttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                 <  2:3576600/76‑1 (MQ=255)
                                   tcttTACCATTTCTGTAAGAGATAAAATAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAttc                                >  2:1330722/1‑75 (MQ=255)
                                   tcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTa                                >  2:340851/1‑76 (MQ=255)
                                   tcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCAGGCTCCGGCCTTCGTTTTCAGCACTGCGTTTTGCttc                                >  2:303361/1‑75 (MQ=255)
                                     tttACCATTTCTGTAAGAGATNAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                              >  1:3276041/1‑76 (MQ=255)
                                     tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                              >  2:4433068/1‑76 (MQ=255)
                                      ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                             <  2:2229050/76‑1 (MQ=255)
                                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                            <  2:367421/76‑1 (MQ=255)
                                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                            >  2:2246381/1‑76 (MQ=255)
                                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                            <  1:733602/76‑1 (MQ=255)
                                       tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                            >  2:4703852/1‑76 (MQ=255)
                                         ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                          <  2:234122/76‑1 (MQ=255)
                                         ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                          <  1:2380398/76‑1 (MQ=255)
                                          cATTTCTGTAAGAGATAACAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGctcg                         >  1:1218092/1‑76 (MQ=255)
                                          cATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGcgcg                         >  1:1205146/1‑76 (MQ=255)
                                           aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcgg                        >  1:4253217/1‑76 (MQ=255)
                                            tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                       >  2:2894803/1‑76 (MQ=255)
                                            tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                       >  2:2302779/1‑76 (MQ=255)
                                              tCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                     <  1:419864/76‑1 (MQ=255)
                                                 gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                  >  1:831207/1‑76 (MQ=255)
                                                 gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                  <  2:4432840/76‑1 (MQ=255)
                                                 gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                  >  1:1225971/1‑76 (MQ=255)
                                                 gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                  <  2:603731/76‑1 (MQ=255)
                                                  tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                 >  2:2358754/1‑76 (MQ=255)
                                                   aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                <  1:1866192/76‑1 (MQ=255)
                                                    agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACac               >  2:2236830/1‑75 (MQ=255)
                                                        gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct            <  2:3902864/73‑1 (MQ=255)
                                                        gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct            <  1:594893/73‑1 (MQ=255)
                                                         ataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctg           >  1:1802882/3‑76 (MQ=255)
                                                          taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc          >  1:1912756/2‑76 (MQ=255)
                                                          taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc          >  1:2013989/2‑76 (MQ=255)
                                                             aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt       >  1:1284515/1‑76 (MQ=255)
                                                              aaaGGCGGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt      <  2:222011/76‑1 (MQ=255)
                                                              aaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATCACAGCctgctgct      >  2:383760/1‑74 (MQ=255)
                                                               aaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc     <  1:2059/76‑1 (MQ=255)
                                                               aaGGCCGGAGCATGCGCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc     <  2:1843363/76‑1 (MQ=255)
                                                                aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa    >  2:1447437/1‑76 (MQ=255)
                                                                 ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCaa   >  2:4337619/1‑76 (MQ=255)
                                                                  gCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCCAt  >  2:2707458/1‑76 (MQ=255)
                                                                  gCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt  >  2:2675674/1‑76 (MQ=255)
                                                                        |                                                                       
GCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAT  >  NC_000913/4478439‑4478582

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