Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,950,734 A→G R576R (CGT→CGC recD ← exodeoxyribonuclease V subunit RecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,950,7340AG98.2% 205.6 / ‑4.5 54R576R (CGT→CGCrecDexodeoxyribonuclease V subunit RecD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (28/25);  total (29/25)
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

CGTAACACTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGACGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGT  >  NC_000913/2950601‑2950865
                                                                                                                                     |                                                                                                                                   
cGTAACACTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGa                                                                                                                               >  2:431593/1‑140 (MQ=255)
 gTAACACTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAc                                                                                                                              <  1:67938/140‑1 (MQ=255)
   aaCACTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAcgc                                                                                                                            >  2:701693/1‑140 (MQ=255)
    acacTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAcgcg                                                                                                                           <  1:837089/140‑1 (MQ=255)
      acTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAcgcgcg                                                                                                                         <  2:42908/140‑1 (MQ=255)
        tCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAcgcgcgcg                                                                                                                       >  2:197995/1‑140 (MQ=255)
         cGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCgg                                                                                                                                  <  2:473854/128‑1 (MQ=255)
         cGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCgg                                                                                                                                  >  1:473854/1‑128 (MQ=255)
          gTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCggg                                                                                                                     >  1:279078/1‑140 (MQ=255)
               tCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAcgc                                                                                                                            >  2:808482/1‑128 (MQ=255)
               tCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGAcgc                                                                                                                            <  1:808482/128‑1 (MQ=255)
                cGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACcg                                                                                                               <  1:301219/140‑1 (MQ=255)
                cGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACcg                                                                                                               <  2:459041/140‑1 (MQ=255)
                cGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACcg                                                                                                               >  2:876593/1‑140 (MQ=255)
                  cATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGt                                                                                                                    >  2:147977/1‑133 (MQ=255)
                  cATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGt                                                                                                                    <  1:147977/133‑1 (MQ=255)
                  cATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACcgcg                                                                                                             >  1:331071/1‑140 (MQ=255)
                  cATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACcgcg                                                                                                             >  1:829129/1‑140 (MQ=255)
                    tCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGt                                                                                                           >  2:781929/1‑140 (MQ=255)
                            ttACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCa                                                                                                   >  2:331558/1‑140 (MQ=255)
                                  ttATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCtctc                                                                                             <  1:781929/140‑1 (MQ=255)
                                  ttATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCtctc                                                                                             <  1:476082/140‑1 (MQ=255)
                                       cccGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGtt                                                                                        >  1:799509/1‑140 (MQ=255)
                                          gTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTtact                                                                                     >  1:917119/1‑140 (MQ=255)
                                            gAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTtactac                                                                                   >  2:295360/1‑140 (MQ=255)
                                            gAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTtactac                                                                                   <  1:876593/140‑1 (MQ=255)
                                                       cgccgcCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGtt                                                                        >  1:607789/1‑140 (MQ=255)
                                                         ccgccAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTgg                                                                      >  2:203976/1‑140 (MQ=255)
                                                          cgccAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGc                                                                     <  2:331071/140‑1 (MQ=255)
                                                           gccAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCt                                                                    <  1:197995/140‑1 (MQ=255)
                                                                    cTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAAt                                                           >  1:50720/1‑140 (MQ=255)
                                                                         ccgcTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAAcg                                                      >  2:604994/1‑140 (MQ=255)
                                                                            cTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAAcgccg                                                   >  1:299586/1‑140 (MQ=255)
                                                                               aGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCAt                                                >  2:158469/1‑140 (MQ=255)
                                                                                   cGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTc                                            <  2:279078/140‑1 (MQ=255)
                                                                                   cGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTc                                            <  2:829129/140‑1 (MQ=255)
                                                                                   cGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTc                                            <  1:431593/140‑1 (MQ=255)
                                                                                   cGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTc                                            <  1:604994/140‑1 (MQ=255)
                                                                                     aGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGa                                          <  2:917119/140‑1 (MQ=255)
                                                                                           gATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCg                                    >  2:59155/1‑140 (MQ=255)
                                                                                           gATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCg                                    >  2:694066/1‑140 (MQ=255)
                                                                                                ccGCACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATccc                               <  2:299586/140‑1 (MQ=255)
                                                                                                   cACTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGc                            <  1:694066/140‑1 (MQ=255)
                                                                                                     cTTAATATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGa                          <  2:50720/140‑1 (MQ=255)
                                                                                                        aaTATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTc                                      >  1:610303/1‑125 (MQ=255)
                                                                                                        aaTATGCGCTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTc                                      <  2:610303/125‑1 (MQ=255)
                                                                                                                cTCATCGGCATACAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCg               >  1:707559/1‑140 (MQ=255)
                                                                                                                cTCATCGGCATACAGCGACAGACGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCg               >  2:361327/1‑140 (MQ=255)
                                                                                                                            cAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACg   <  2:607789/140‑1 (MQ=255)
                                                                                                                            cAGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACg   <  1:227007/140‑1 (MQ=255)
                                                                                                                             aGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTa                  <  2:353082/124‑1 (MQ=255)
                                                                                                                             aGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTa                  >  1:353082/1‑124 (MQ=255)
                                                                                                                             aGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGt  >  1:6709/1‑140 (MQ=255)
                                                                                                                             aGCGACAGGCGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGt  >  1:328895/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CGTAACACTCGTACGTCGCATCCGGCAATTACGTTTATTCCCGTGAACTAAACAACGCCGCCAGACCACTGCGCCGCTCAGTACGAGTGGCGATTGCCGCACTTAATATGCGCTCATCGGCATACAGCGACAGACGGCGACGCGCGCGGGTCACCGCGGTATAAACCAGCTCTCGCGTTACTACCGGCGTGCGTTGGCTCGGCAAAATCAACGCCGCATGGTCGAACTCCGATCCCTGCGATTTATGTACCGTCATCGCCCACGT  >  NC_000913/2950601‑2950865

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