Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,477,041 Δ1 bp intergenic (‑151/‑269) b4256 ← / → yjgN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,477,0410G.98.4% 553.5 / ‑1.9 126intergenic (‑151/‑269)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (34/90);  total (34/92)
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

Only 100 of 128 total aligned reads displayed.

CGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCC  >  NC_000913/4476971‑4477110
                                                                      |                                                                      
ggccGGTGTGGCCTTGATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGAtt                                                                  <  2:206190/73‑1 (MQ=255)
cGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATACGGCGGATACCTGACGCATGACGCGGGCGATAGCAGGAtt                                                                  >  2:1463704/1‑76 (MQ=255)
  gcGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTAt                                                                <  1:1090617/75‑1 (MQ=255)
  aCGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTAt                                                                <  1:5653419/76‑1 (MQ=255)
    ggTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATcc                                                              >  2:1700146/1‑76 (MQ=255)
    ggTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATcc                                                              >  1:4883828/1‑76 (MQ=255)
     gtgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCt                                                             <  1:4209125/76‑1 (MQ=255)
     gtgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCt                                                             <  2:2352158/76‑1 (MQ=255)
      tgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTg                                                            <  2:5151990/76‑1 (MQ=255)
      tgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTg                                                            <  1:2015057/76‑1 (MQ=255)
       gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCCGGATTATCCTGc                                                           >  2:31259/1‑76 (MQ=255)
       gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGc                                                           >  1:3046510/1‑76 (MQ=255)
       gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGc                                                           >  1:2094399/1‑76 (MQ=255)
       ggaGCCTTTATCAGCGGGAAGACCGTATTCGGCGGATACCTGGCGAATGACGCGGGCGATAGCAGGATTATCCTGc                                                           <  2:111314/74‑1 (MQ=255)
        tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCa                                                          <  1:1530096/76‑1 (MQ=255)
        tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCa                                                          >  1:3522800/1‑76 (MQ=255)
        tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCa                                                          <  1:346212/76‑1 (MQ=255)
         aGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAg                                                         <  1:2544058/76‑1 (MQ=255)
          gCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCCGGATTATCCTGCCGc                                                        >  2:172755/1‑76 (MQ=255)
          gCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGc                                                        <  2:3183671/76‑1 (MQ=255)
          gCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGc                                                        >  2:1640773/1‑76 (MQ=255)
            cTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGt                                                      <  1:4950934/76‑1 (MQ=255)
            cTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGt                                                      <  2:188528/76‑1 (MQ=255)
             tttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTc                                                     <  2:4680390/76‑1 (MQ=255)
             tttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTc                                                     >  2:4304749/1‑76 (MQ=255)
             tttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTc                                                     >  1:1397867/1‑76 (MQ=255)
              ttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCa                                                    <  2:1885762/76‑1 (MQ=255)
               tATCAGCGTTAAGGCCGTATTCGGCGGATCCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAg                                                   <  2:181139/76‑1 (MQ=255)
               tATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAg                                                   <  2:5310218/76‑1 (MQ=255)
               tATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAg                                                   <  2:3143974/76‑1 (MQ=255)
               tATCAGCGGTAAGACCGGATTCGGCGGATACCTGACGACTGACGCGGGCGATAGCAGGGTTATCCTGCAGCGTCAg                                                   >  2:706195/1‑76 (MQ=255)
                aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAgg                                                  <  2:4858112/76‑1 (MQ=255)
                aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAgg                                                  <  2:2759579/76‑1 (MQ=255)
                aTCACCGATAAGACCTTATACGGCGGATACCTGACGAATGCCGCGTGCTCTAGCAGGACTATCCTGCATCGTCAgg                                                  <  1:5044632/76‑1 (MQ=255)
                 tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                 <  1:3148516/76‑1 (MQ=255)
                  cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcg                                                <  1:2656795/76‑1 (MQ=255)
                  cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcg                                                <  2:1213795/76‑1 (MQ=255)
                  cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcg                                                <  2:80514/76‑1 (MQ=255)
                   aGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcga                                               >  1:3102997/1‑76 (MQ=255)
                     cGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacg                                             >  1:281892/1‑76 (MQ=255)
                      ggTAAGACCGTATTCGGCGGATGCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgc                                            <  1:2807350/76‑2 (MQ=255)
                      ggTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgc                                            >  1:788257/1‑75 (MQ=255)
                      ggTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgc                                            <  1:1022861/76‑2 (MQ=255)
                           gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTAc                                       >  2:3951259/1‑76 (MQ=255)
                             ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                     <  2:2541291/76‑1 (MQ=255)
                             ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                     >  2:2737224/1‑76 (MQ=255)
                             ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                     <  2:3102997/76‑1 (MQ=255)
                               gTATTCGGCGGATACCTGACGAATGACGCGGGTGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                   <  1:3850921/76‑1 (MQ=255)
                               gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                   >  2:2472888/1‑76 (MQ=255)
                               gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTAAAGGt                                   >  1:311865/1‑76 (MQ=255)
                                tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTg                                  <  1:5498250/76‑1 (MQ=255)
                                  ttCGGCGGATACCTGCCGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACATGTGAt                                <  1:1096466/76‑1 (MQ=255)
                                  ttCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                <  1:5746444/76‑1 (MQ=255)
                                   gcggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                               <  2:335454/75‑1 (MQ=255)
                                    cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTg                              <  2:4170384/76‑1 (MQ=255)
                                     ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                             >  1:3839462/1‑76 (MQ=255)
                                     ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                             >  2:3260897/1‑76 (MQ=255)
                                     ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGCGACTgc                             >  2:1740792/1‑76 (MQ=255)
                                            aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGgcgcca                      <  2:1864855/76‑3 (MQ=255)
                                            aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                      <  2:4883827/76‑1 (MQ=255)
                                            aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                      <  1:3171742/76‑1 (MQ=255)
                                             ccTGACTAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                     <  1:2220056/76‑1 (MQ=255)
                                             ccTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                     <  1:19286/76‑1 (MQ=255)
                                             ccTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                     <  1:361005/76‑1 (MQ=255)
                                             ccTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                     <  1:262312/76‑1 (MQ=255)
                                              cTGACGAGTGACGCGGGCGATAGGAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                    <  2:3052689/76‑1 (MQ=255)
                                              cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                >  2:2335051/1‑64 (MQ=255)
                                              cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                <  1:2335051/64‑1 (MQ=255)
                                              cTGACGAATAACGCGGGCGGTAGCAGGATGATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                    <  2:3021201/76‑1 (MQ=255)
                                                gACGAATGACTCGTGCGTTAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGtt                  <  2:3651079/76‑1 (MQ=255)
                                                gACGAATGACGCGGGCGATAGCATGATTATCCTGCGGCGTCAGGCGACGCATTGCAGGTGATTGCGGCGCAATGtt                  <  1:4607719/76‑1 (MQ=255)
                                                gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                             >  2:5801250/1‑65 (MQ=255)
                                                gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGtt                  <  1:5138595/76‑1 (MQ=255)
                                                 aCGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                 <  1:1637162/76‑1 (MQ=255)
                                                  ggAAGGACGGGGGCGATAGGAGGATTATCCTACAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAt                <  2:5129799/75‑1 (MQ=255)
                                                   gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGCCGCATTACAGGTGATTGCGGCGCAATGTTAtt               <  2:1076214/76‑1 (MQ=255)
                                                   gAATGACGCGGGCGATAGCAGGAGTATCCTGCAGCGTCAGGCGCCGCATTGCTGGTGATTGTGGCGCGCTGTTAtt               <  2:1122822/76‑1 (MQ=255)
                                                    aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAttt              >  1:1117532/1‑75 (MQ=255)
                                                    aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc              <  2:281892/76‑1 (MQ=255)
                                                    aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc              <  2:3852652/76‑1 (MQ=255)
                                                    aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc              <  1:1900909/76‑1 (MQ=255)
                                                     aTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa             >  1:1452537/1‑76 (MQ=255)
                                                      tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt            <  1:5590518/76‑1 (MQ=255)
                                                      tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt            <  1:664469/76‑1 (MQ=255)
                                                       gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa           <  1:3286893/76‑1 (MQ=255)
                                                       gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa           <  2:2101086/76‑1 (MQ=255)
                                                        aCGCGGGCGATAGCAGGATTATCCTGCATCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAg          <  2:3649993/76‑1 (MQ=255)
                                                         cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt         <  1:3951260/76‑1 (MQ=255)
                                                         cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt         <  1:5806381/76‑1 (MQ=255)
                                                         cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt         <  2:1414661/76‑1 (MQ=255)
                                                         cgcgGGCGATAGCAGGATGATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt         <  1:2640525/76‑1 (MQ=255)
                                                          gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt        <  1:4553067/76‑1 (MQ=255)
                                                           cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCCGCGCAATGTTATTCATAGttt       >  2:856091/1‑76 (MQ=255)
                                                             ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa     <  2:4423599/76‑1 (MQ=255)
                                                             ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa     <  1:3629538/76‑1 (MQ=255)
                                                             ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa     <  1:327221/76‑1 (MQ=255)
                                                              gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAc    <  1:250940/76‑1 (MQ=255)
                                                               cGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc   <  2:2100745/76‑1 (MQ=255)
                                                               cGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc   <  1:189223/76‑1 (MQ=255)
                                                                gATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc  <  2:1051811/76‑1 (MQ=255)
                                                                      |                                                                      
CGACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCC  >  NC_000913/4476971‑4477110

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