Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,960,844 T→C T529A (ACG→GCG)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,8440TC97.4% 137.1 / ‑3.9 39T529A (ACG→GCG) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (22/16);  total (23/16)
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.88e-01

CCCGCGACGCCAGATGTTTAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCG  >  NC_000913/2960698‑2960985
                                                                                                                                                  |                                                                                                                                             
cccGCGACGCCAGATGTTTAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAgg                                                                                                                                           >  2:169002/1‑151 (MQ=255)
cccGCGACGCCAGATGTTTAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAgg                                                                                                                                           >  1:469597/1‑151 (MQ=255)
         ccAGATGTTTAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGc                                                                                                                                  >  2:444078/1‑151 (MQ=255)
           aGATGTTTAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCa                                                                                                                                <  2:389718/151‑1 (MQ=255)
                tttAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGgtgt                                                                                                                           >  1:396079/1‑151 (MQ=255)
                 ttAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGgtgt                                                                                                                           >  1:249518/1‑150 (MQ=255)
                  tAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGgtgtg                                                                                                                          >  2:577297/1‑150 (MQ=255)
                             gCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGc                                                                                                              >  2:287008/1‑151 (MQ=255)
                               aGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCgg                                                                                                            <  1:61949/151‑1 (MQ=255)
                                  gAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCTCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCgggg                                                                                                          >  1:598128/1‑150 (MQ=255)
                                  gAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCgggg                                                                                                          >  1:470347/1‑150 (MQ=255)
                                   aaGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCgggg                                                                                                          >  1:287570/1‑149 (MQ=255)
                                        aGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTcc                                                                                                    >  2:465232/1‑150 (MQ=255)
                                         gATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTcc                                                                                                    >  2:268620/1‑149 (MQ=255)
                                             ccccACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTcg                                                                                              >  1:502545/1‑151 (MQ=255)
                                             ccccACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTcg                                                                                              >  2:423304/1‑151 (MQ=255)
                                                  ccaGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAAc                                                                                         <  2:438200/151‑1 (MQ=255)
                                                    aGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGt                                                                                       <  2:181545/151‑1 (MQ=255)
                                                                         gATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAAc                                                                  >  1:548146/1‑151 (MQ=255)
                                                                            tcatcaTAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGa                                                                >  1:183469/1‑150 (MQ=255)
                                                                             catcatAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAAt                                                              >  1:225658/1‑151 (MQ=255)
                                                                                 atAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCg                                                          <  1:515099/151‑1 (MQ=255)
                                                                                      gTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTc                                                     <  1:558469/151‑1 (MQ=255)
                                                                                        aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTcct                                                   >  2:474420/1‑149 (MQ=255)
                                                                                        aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGt                                                   >  2:118162/1‑151 (MQ=255)
                                                                                         gAACCGATTGTCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                  >  1:583266/1‑151 (MQ=255)
                                                                                              gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAccc                                             <  2:360256/151‑1 (MQ=255)
                                                                                              gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAccc                                             >  1:84540/1‑151 (MQ=255)
                                                                                                  gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                          <  2:225653/150‑1 (MQ=255)
                                                                                                  gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                          <  1:380563/150‑1 (MQ=255)
                                                                                                  gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTg                                         <  1:287015/151‑1 (MQ=255)
                                                                                                    cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGc                                       >  1:510549/1‑151 (MQ=255)
                                                                                                         gCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                  <  1:56305/151‑1 (MQ=255)
                                                                                                              gcgcgcTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAAc                              <  2:141586/150‑1 (MQ=255)
                                                                                                                 cgcTCTCCATCGCGTAGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAg                          <  2:488549/151‑1 (MQ=255)
                                                                                                                                aGCCCAACAACATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGAt           <  2:412648/151‑1 (MQ=255)
                                                                                                                                    caacaaCATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc        >  1:40371/1‑150 (MQ=255)
                                                                                                                                     aacaacATACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCaa      <  2:510796/151‑1 (MQ=255)
                                                                                                                                         acaTACGCGCCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcg  <  2:542828/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCCGCGACGCCAGATGTTTAGCTGCATTAGCAGTGAAGCCAGATGCCCCACCAGTTCTGCAATTAAGCCGCTCGATTCATCATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCG  >  NC_000913/2960698‑2960985

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