Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,477,068 +C intergenic (‑178/‑242) 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,0681.C98.9% 230.4 / ‑4.3 91intergenic (‑178/‑242)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (33/57);  total (33/58)
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

CGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTG  >  NC_000913/4476992‑4477137
                                                                              |                                                                     
cGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacg                                                                        >  1:3437306/1‑76 (MQ=255)
 ggTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgc                                                                       >  1:3004857/1‑75 (MQ=255)
      gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTAc                                                                  >  1:1564502/1‑76 (MQ=255)
       aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                 <  1:2257034/76‑1 (MQ=255)
        ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                                                <  1:3481326/76‑1 (MQ=255)
        ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                                                <  2:1761506/76‑1 (MQ=255)
         cggaTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAgg                                                               <  2:3704070/73‑1 (MQ=255)
         cGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGTGACGCATTACAgg                                                               <  2:4023304/76‑1 (MQ=255)
          gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAggg                                                              >  1:634606/1‑75 (MQ=255)
          gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                                              >  2:2972790/1‑76 (MQ=255)
            aTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                                            >  2:3960893/1‑76 (MQ=255)
              tCGGCGGATACCTGGCGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                                                          <  2:446383/76‑1 (MQ=255)
              tCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                                                          <  1:3837435/76‑1 (MQ=255)
               cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTg                                                         <  2:3024876/76‑1 (MQ=255)
                ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                                                        >  2:2343336/1‑76 (MQ=255)
                ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                                                        <  2:1478775/76‑1 (MQ=255)
                ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                                                        <  1:3298575/76‑1 (MQ=255)
                 gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcg                                                       >  1:1747763/1‑76 (MQ=255)
                  cggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                                                      >  1:875595/1‑76 (MQ=255)
                  cggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                                                      <  1:1159088/76‑1 (MQ=255)
                   ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggc                                                     <  2:4809603/76‑1 (MQ=255)
                   ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggc                                                     <  2:3605509/76‑1 (MQ=255)
                   ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTccggc                                                     >  2:176809/1‑76 (MQ=255)
                     aTACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGgcgc                                                   <  1:2685434/76‑1 (MQ=255)
                       aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                 >  1:2972853/1‑76 (MQ=255)
                       aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                 >  1:393854/1‑76 (MQ=255)
                        ccTGACGAGTCAGGCGGGCGAGAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGATGATTGCGGCGCAAt                                                <  1:3320466/76‑1 (MQ=255)
                         cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                                               >  2:2731245/1‑76 (MQ=255)
                          tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGTGACGCATTACAGGTGATTGCGGCGCAATGt                                              >  2:3424175/1‑76 (MQ=255)
                          tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGt                                              <  2:2725327/76‑1 (MQ=255)
                              gactGACGCGGGCGATAGGTGGATTATCCTGCAGCGTCAGGCGACGCATTCCAGGAGATTGCGGCGCAATGTTCtt                                          <  2:3000782/73‑1 (MQ=255)
                              gAATGGCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                          <  2:4700800/76‑1 (MQ=255)
                              gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                          <  2:1636352/76‑1 (MQ=255)
                              gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                          >  1:1819241/1‑76 (MQ=255)
                              gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                          >  1:4534887/1‑76 (MQ=255)
                              gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                          <  1:2963121/76‑1 (MQ=255)
                               aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                                         >  2:2520471/1‑76 (MQ=255)
                               aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                                         <  1:2576690/76‑1 (MQ=255)
                                aTGACGCGGGCGACAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                        >  2:1396884/1‑76 (MQ=255)
                                 tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt                                       >  2:829264/1‑76 (MQ=255)
                                   aCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAg                                     <  2:2572941/76‑1 (MQ=255)
                                    cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                                    <  2:3088664/76‑1 (MQ=255)
                                    cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                                    <  1:3656851/76‑1 (MQ=255)
                                     gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                                   >  1:1711655/1‑76 (MQ=255)
                                     gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                                   >  1:3697519/1‑76 (MQ=255)
                                     gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                                   >  2:3755647/1‑76 (MQ=255)
                                      cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt                                  >  2:706424/1‑76 (MQ=255)
                                      cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt                                  <  1:3367764/76‑1 (MQ=255)
                                      catgCGATAGCAGGATTATCCGGCAGCGTCAGGCGACGCATGACAGGTGATTGCGGCGCAAGGTTATTCATAGttt                                  <  1:3392416/73‑1 (MQ=255)
                                       gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                 >  1:1845650/1‑76 (MQ=255)
                                       gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                 >  1:2323139/1‑76 (MQ=255)
                                       gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                 <  2:3517267/76‑1 (MQ=255)
                                        ggCGATAGCAGGGTTATCCTGCAGCGTCAGGCGACGGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                                <  2:4580832/76‑1 (MQ=255)
                                        ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                                >  2:248973/1‑76 (MQ=255)
                                         gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAc                               <  1:1564600/76‑1 (MQ=255)
                                          cGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTGCAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc                              <  1:3843536/76‑1 (MQ=255)
                                           gATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc                             <  1:537851/76‑1 (MQ=255)
                                                gcaGGATTATCCTGCAGCGTCAGTCGACGCGTTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTACCGCaata                         <  1:3718377/74‑1 (MQ=255)
                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata                         <  2:1821172/74‑1 (MQ=255)
                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata                         >  2:4736114/3‑76 (MQ=255)
                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata                         <  1:871569/74‑1 (MQ=255)
                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata                         <  2:1956799/74‑1 (MQ=255)
                                                 caGGATTATCCTGCAGCTTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa                        <  2:3165854/75‑1 (MQ=255)
                                                  aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaa                       >  1:2240957/1‑76 (MQ=255)
                                                    gATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaa                     <  2:3033430/76‑1 (MQ=255)
                                                     aTTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaat                    >  1:1324261/1‑76 (MQ=255)
                                                     aTTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaat                    <  1:4689305/76‑1 (MQ=255)
                                                        aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaataaa                 <  1:166691/76‑1 (MQ=255)
                                                        aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaataaa                 <  2:1678163/76‑1 (MQ=255)
                                                           cTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa              <  2:3727400/76‑1 (MQ=255)
                                                           cTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa              <  2:4711035/76‑1 (MQ=255)
                                                            tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTaa             <  1:972485/76‑1 (MQ=255)
                                                            tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTaa             <  2:759321/76‑1 (MQ=255)
                                                            tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTaa             <  2:3105074/76‑1 (MQ=255)
                                                             gCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAt            >  1:1225863/1‑76 (MQ=255)
                                                              cAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAtt           >  2:4521844/1‑76 (MQ=255)
                                                               agggTCAGGCGACGCATTTCAGGTGATTTCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta          <  2:1647837/73‑1 (MQ=255)
                                                               aGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta          <  2:3877639/76‑1 (MQ=255)
                                                                gCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtat         <  2:3697519/76‑1 (MQ=255)
                                                                gCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtat         <  2:2816120/76‑1 (MQ=255)
                                                                 cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata        <  2:3193306/76‑1 (MQ=255)
                                                                 cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata        >  1:1735627/1‑76 (MQ=255)
                                                                 cGTCACGCGACGCATTAAAGGTGATTGCTACGCAATTTTATTCATAGTTTACCCCAATAAAAATAACGTAATtatc        <  1:2110833/76‑2 (MQ=255)
                                                                  gTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATAc       <  1:4285004/76‑1 (MQ=255)
                                                                   tCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa      <  2:4278099/76‑1 (MQ=255)
                                                                   tCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa      >  2:3038774/1‑76 (MQ=255)
                                                                    cAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAt     <  1:2872748/76‑1 (MQ=255)
                                                                     aGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAtt    <  2:1608839/76‑1 (MQ=255)
                                                                     aGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAtt    <  1:3027433/76‑1 (MQ=255)
                                                                     aGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATGCATAGGTTAACCCAATAAAAATAAAGTAAGGATACAtt    <  1:3784596/76‑1 (MQ=255)
                                                                     aGGCGACGCAGTACAGGTGATTGCGGCGCAATGTTATTCATAGATTAACCCAACAAAAATAAAGCACTCATACAtt    >  2:222794/1‑76 (MQ=255)
                                                                       gCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg  <  1:724444/76‑1 (MQ=255)
                                                                       gCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg  >  2:800703/1‑76 (MQ=255)
                                                                              |                                                                     
CGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTG  >  NC_000913/4476992‑4477137

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