Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,337,024 G→A 100% S3L (TCG→TTG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,337,0240GA99.0% 273.3 / ‑5.9 79S3L (TCG→TTG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/45);  minor base C (0/1);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336927‑4337120
                                                                                                 |                                                                                                
cAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCaaa                                                                                                <  2:1759491/100‑1 (MQ=255)
 aGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAg                                                                                               >  1:65097/1‑100 (MQ=255)
   gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             >  1:1167467/1‑100 (MQ=255)
   gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  1:900514/100‑1 (MQ=255)
   gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  1:604927/100‑1 (MQ=255)
   gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  1:2202686/100‑1 (MQ=255)
     ttGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACAt                                                                                           <  2:147754/100‑1 (MQ=255)
     ttGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACAt                                                                                           <  2:861005/100‑1 (MQ=255)
          ggTAACAGAAAAGCACCTGACGCCATAATATTCCCCGACACCAGCAGGGTGACGCGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      <  2:1979047/100‑1 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      >  1:1739675/1‑100 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      >  1:597858/1‑100 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      >  2:1628168/1‑100 (MQ=255)
           gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaa                                                                                     <  2:1112104/100‑1 (MQ=255)
           gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaa                                                                                     <  1:424250/100‑1 (MQ=255)
            tAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaaa                                                                                    >  1:2197558/1‑100 (MQ=255)
            tAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaaa                                                                                    >  1:579927/1‑100 (MQ=255)
             aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                   <  2:2060411/100‑1 (MQ=255)
              aCAGAAAAACACCTGACCCCATAATATTCCCCGACGCCAGCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACt                                                                                  <  2:1830156/100‑1 (MQ=255)
                aGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTcc                                                                                >  1:302589/1‑100 (MQ=255)
                 gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCt                                                                               >  1:1718012/1‑100 (MQ=255)
                  aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTg                                                                              <  1:429926/100‑1 (MQ=255)
                    aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                            <  2:1041793/100‑1 (MQ=255)
                       cacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAg                                                                         <  1:800787/100‑1 (MQ=255)
                        acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAgg                                                                        <  1:1464134/100‑1 (MQ=255)
                         ccTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGc                                                                       <  2:1936647/100‑1 (MQ=255)
                           tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGa                                                                     <  1:2208899/100‑1 (MQ=255)
                              ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                  >  2:459362/1‑100 (MQ=255)
                               cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCACAATTTGTGAGCATCAGCATCCAAAGAAATAATTAAACTACTGCGAAGGCGAGCtt                                                                 >  2:1869442/1‑100 (MQ=255)
                                  ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                              <  1:1840680/100‑1 (MQ=255)
                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:1485113/100‑1 (MQ=255)
                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             >  2:1291552/1‑100 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTCAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                           <  2:199174/100‑1 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTTATCTCCTGCGAAGGCGAGCTTCg                                                               <  2:2220402/96‑1 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                           <  1:1548961/100‑1 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                           <  2:1591803/100‑1 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                           <  2:1818271/100‑1 (MQ=255)
                                      tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                          <  1:1033283/100‑1 (MQ=255)
                                      tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                          <  1:45963/100‑1 (MQ=255)
                                       atTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACaa                                                         <  1:2273830/100‑1 (MQ=255)
                                          ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:64740/1‑100 (MQ=255)
                                           cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaaa                                                     >  1:1694934/1‑100 (MQ=255)
                                            ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAt                                                    >  2:149188/1‑100 (MQ=255)
                                              gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                  <  2:2326875/100‑1 (MQ=255)
                                               acacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                 <  2:2333205/100‑1 (MQ=255)
                                                cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                >  2:1474114/1‑100 (MQ=255)
                                                 acCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                               <  2:2197539/100‑1 (MQ=255)
                                                  ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              >  1:4352/1‑100 (MQ=255)
                                                  ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              >  2:2101833/1‑100 (MQ=255)
                                                   cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                             <  2:71519/100‑1 (MQ=255)
                                                   cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                             >  2:1511040/1‑100 (MQ=255)
                                                   cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                             >  2:880127/1‑100 (MQ=255)
                                                       aGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatc                                         <  2:752096/100‑1 (MQ=255)
                                                       aGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACAtcatatc                                         >  2:490193/1‑100 (MQ=255)
                                                        gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                        >  2:982208/1‑100 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGTATCCCAAGACATACTTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  <  1:106806/100‑1 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  2:599183/1‑100 (MQ=255)
                                                                gggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                <  2:1786295/100‑1 (MQ=255)
                                                                gggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                >  1:599646/1‑100 (MQ=255)
                                                                    ttATGCGCACTTTGTATGCATCAGCATCCAAATACATAATTAAACCCCTGCGAAGGCGAGCTTTGTAACAATAAATTACGTCAGATCTTAGTCCTgt                               <  2:914836/97‑2 (MQ=255)
                                                                     tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                           <  1:492988/100‑1 (MQ=255)
                                                                       aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  1:907707/1‑100 (MQ=255)
                                                                       aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  1:843541/1‑100 (MQ=255)
                                                                       aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGACCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  2:1506070/1‑100 (MQ=255)
                                                                         cccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATa                       <  1:1511671/100‑1 (MQ=255)
                                                                          ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                      <  2:261186/100‑1 (MQ=255)
                                                                          ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                      >  1:452539/1‑100 (MQ=255)
                                                                            aCTTTGTGAGCATCAGCATCCAAAGACATCATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGtt                    <  1:1238287/100‑1 (MQ=255)
                                                                             cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  1:887580/100‑1 (MQ=255)
                                                                             cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  2:429786/100‑1 (MQ=255)
                                                                             cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  1:297698/100‑1 (MQ=255)
                                                                                tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  1:588388/100‑1 (MQ=255)
                                                                                tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  1:1982606/100‑1 (MQ=255)
                                                                                    agcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt            <  2:1227662/100‑1 (MQ=255)
                                                                                     gcatcagcatcCAAAGACATAATTAAACTCATGAGAAGGAGAGAATCGTGACAATACATTACCACCTATAATACGCCTGCATTCAATAGTTGAAtcatta           >  2:437323/1‑100 (MQ=255)
                                                                                        tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        <  2:1167462/100‑1 (MQ=255)
                                                                                          agcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc      <  2:1101577/100‑1 (MQ=255)
                                                                                          agcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc      >  1:1312742/1‑100 (MQ=255)
                                                                                           gcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATAACATACGCCTGCATTCGATAGTTTAATTATTATTGAc      >  2:2189799/1‑99 (MQ=255)
                                                                                            catcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg    >  1:328145/1‑100 (MQ=255)
                                                                                              tcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg  >  2:2130327/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336927‑4337120

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: ‑