Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,698,281 A→G 67.3% V300V (GTT→GTC)  malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,698,2810AG67.3% 86.7 / 56.5 70V300V (GTT→GTC) malItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base A (12/11);  new base G (24/23);  total (36/34)
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

ATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTT  >  NC_000913/1698183‑1698378
                                                                                                  |                                                                                                 
ataaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTaa                                                                                                  >  2:25721/1‑100 (MQ=255)
   aGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAg                                                                                               <  1:1562153/100‑1 (MQ=255)
    gATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAgg                                                                                              >  2:306118/1‑100 (MQ=255)
    gATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAgg                                                                                              <  2:1682420/100‑1 (MQ=255)
     aTTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGt                                                                                             >  2:1667076/1‑100 (MQ=255)
         cGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTatat                                                                                         >  1:99071/1‑100 (MQ=255)
          gTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATatc                                                                                        <  1:623788/100‑1 (MQ=255)
           tGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATatca                                                                                       >  2:186184/1‑100 (MQ=255)
           gatctGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATatca                                                                                       <  2:1158180/96‑1 (MQ=255)
            gAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATatcat                                                                                      >  1:1833450/1‑100 (MQ=255)
               tGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATatcatc                                                                                     <  2:88806/98‑1 (MQ=255)
               tGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCaa                                                                                   >  2:1958633/1‑100 (MQ=255)
                gCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATatcaagatc                                                                                  >  1:1158190/1‑95 (MQ=255)
                 cGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAgt                                                                                 <  2:2649361/100‑1 (MQ=255)
                 cGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAgt                                                                                 >  2:1802897/1‑100 (MQ=255)
                     tctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGtgg                                                                             >  1:1733418/1‑100 (MQ=255)
                     tctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGtgg                                                                             >  1:1730125/1‑100 (MQ=255)
                      ctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggt                                                                            >  2:1762616/1‑100 (MQ=255)
                      ctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGtggt                                                                            <  2:765887/100‑1 (MQ=255)
                        tCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggttg                                                                          <  1:2594029/100‑1 (MQ=255)
                            gggTGATTTTTTGCATCATGCGATCCGCACGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCTGGTGACAGGTATATCATCAAGTGTGGTTGgtgt                                                                      <  2:1588323/100‑1 (MQ=255)
                            gggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGgtgt                                                                      <  2:321970/100‑1 (MQ=255)
                               tGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCCGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGc                                                                   >  1:697742/1‑100 (MQ=255)
                                 aTTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCAt                                                                 <  2:628387/100‑1 (MQ=255)
                                  ttttttGCATTTTTTTATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATc                                                                >  2:2413748/1‑100 (MQ=255)
                                      ttGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCAACAAGTGTGGTTGGTGTCGCATCGGTa                                                            <  1:1887104/100‑1 (MQ=255)
                                       tGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTaa                                                           >  2:600293/1‑100 (MQ=255)
                                          atcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATg                                                        <  1:172006/100‑1 (MQ=255)
                                          atcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATg                                                        >  1:110533/1‑100 (MQ=255)
                                             atGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccg                                                     <  1:1588975/100‑1 (MQ=255)
                                             atGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccg                                                     <  2:415382/100‑1 (MQ=255)
                                              tGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgc                                                    <  1:364396/100‑1 (MQ=255)
                                                 gATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAg                                                 <  2:694371/100‑1 (MQ=255)
                                                   tCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCg                                               <  2:1988111/100‑1 (MQ=255)
                                                     cGCAAGTGTGATCCCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGaa                                             >  1:1738609/1‑100 (MQ=255)
                                                     cGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGaa                                             >  1:389750/1‑100 (MQ=255)
                                                        aaGTGTGATACCAAGTTCCCGCGCTGGTGTGCTGGCTCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGaaa                                            <  1:513958/98‑1 (MQ=255)
                                                        aaGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAAcc                                          <  1:865139/100‑1 (MQ=255)
                                                         aGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCt                                         <  2:1002957/100‑1 (MQ=255)
                                                           tgtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGt                                       >  1:120166/1‑100 (MQ=255)
                                                           tgtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGt                                       >  1:2014043/1‑100 (MQ=255)
                                                            gtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGtt                                      >  1:346497/1‑100 (MQ=255)
                                                             tgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTg                                     <  2:1956529/100‑1 (MQ=255)
                                                              gATACCAAGTTCCCGCGCTTGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGc                                    >  1:1500381/1‑100 (MQ=255)
                                                               aTACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCt                                   <  1:186186/100‑1 (MQ=255)
                                                                 aCCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCa                                 >  1:2469251/1‑100 (MQ=255)
                                                                 aCCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTATATGCCGCCAGCGAAACCTGTTGCTCa                                 <  1:1615542/100‑1 (MQ=255)
                                                                 aCCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCa                                 >  1:640924/1‑100 (MQ=255)
                                                                     aGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGt                             <  2:80980/100‑1 (MQ=255)
                                                                     aGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGt                             <  2:1895585/100‑1 (MQ=255)
                                                                       ttCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTaa                           >  1:2534808/1‑100 (MQ=255)
                                                                       ttCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTaa                           >  1:220800/1‑100 (MQ=255)
                                                                       ttCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCACAACGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTaa                           >  1:683572/1‑100 (MQ=255)
                                                                        tCCCGCGCTGGCGTGCTGGCCTAGGTAACATTTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAAc                          <  1:838465/100‑1 (MQ=255)
                                                                        tCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAAc                          <  1:321108/100‑1 (MQ=255)
                                                                         cccGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACg                         >  2:621019/1‑100 (MQ=255)
                                                                            gcgcTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATc                      <  1:423568/100‑1 (MQ=255)
                                                                               cTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGAc                   <  2:2487406/100‑1 (MQ=255)
                                                                                 ggCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGGCGCATCGGTAAATTCCGCCAACGAAACCAGATGCTAAAAGAAACGATCGACTc                 >  1:1069150/1‑100 (MQ=255)
                                                                                  gCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTcc                <  2:992800/100‑1 (MQ=255)
                                                                                   cGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCg               >  2:88704/1‑100 (MQ=255)
                                                                                    gTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGc              <  2:2535345/100‑1 (MQ=255)
                                                                                        tGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTc          <  2:1909701/100‑1 (MQ=255)
                                                                                        tGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTc          >  2:1259076/1‑100 (MQ=255)
                                                                                         ggCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcc         >  2:2423429/1‑100 (MQ=255)
                                                                                         ggCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcc         >  2:2092343/1‑100 (MQ=255)
                                                                                             cAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGc     >  2:2381927/1‑100 (MQ=255)
                                                                                               ggTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCtt   >  1:388013/1‑100 (MQ=255)
                                                                                               ggTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCtt   <  2:647141/100‑1 (MQ=255)
                                                                                                gTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCttt  >  2:2033812/1‑100 (MQ=255)
                                                                                                gTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCttt  >  2:2031362/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
ATAAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTT  >  NC_000913/1698183‑1698378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑