Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,281,193 C→T 6.9% R188C (CGT→TGT)  radD → putative DNA repair helicase RadD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,281,1930CT 6.9% 177.5 / 2.6 58R188C (CGT→TGT) radDputative DNA repair helicase RadD
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base T (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAA  >  NC_000913/2281050‑2281336
                                                                                                                                               |                                                                                                                                               
tGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTtata                                                                                                                                           <  1:1764982/150‑1 (MQ=255)
  aaTCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATg                                                                                                                                         <  2:1040412/150‑1 (MQ=255)
    tCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGAt                                                                                                                                       <  1:38860/150‑1 (MQ=255)
     ccccACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGAtt                                                                                                                                      >  2:1956340/1‑150 (MQ=255)
      cccACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTa                                                                                                                                     <  1:950354/150‑1 (MQ=255)
          cTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAAAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAaca                                                                                                                                 >  2:1105582/1‑150 (MQ=255)
              cGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGc                                                                                                                             <  1:698264/150‑1 (MQ=255)
                  tgctgGGGCTGACTGACACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAAAACgcataaa                                                                                                                         >  2:88093/1‑144 (MQ=255)
                   gctgGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCt                                                                                                                        >  1:1614568/1‑150 (MQ=255)
                   gctgGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCt                                                                                                                        <  1:1067365/150‑1 (MQ=255)
                    ctgGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTg                                                                                                                       >  1:1812883/1‑150 (MQ=255)
                        ggCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcgcc                                                                                                                  <  2:681339/150‑1 (MQ=255)
                          cggACTGCCACGCCGTTTCGATTGGGCAAAGGCGGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcgccg                                                                                                                 <  2:1294725/148‑1 (MQ=255)
                          cTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcgccg                                                                                                                 <  2:269883/150‑1 (MQ=255)
                           tGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcgccgc                                                                                                                <  2:196682/150‑1 (MQ=255)
                            gACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACgccgcc                                                                                                               >  2:1557599/1‑150 (MQ=255)
                            gACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACgccgcc                                                                                                               <  2:1091138/150‑1 (MQ=255)
                            gACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACgccgcc                                                                                                               <  2:1075018/150‑1 (MQ=255)
                                gCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGaa                                                                                                           >  1:329475/1‑150 (MQ=255)
                                 ccACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAAc                                                                                                          <  1:1067285/150‑1 (MQ=255)
                                  cACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACg                                                                                                         <  2:1614568/150‑1 (MQ=255)
                                       cTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTgg                                                                                                    >  1:554877/1‑150 (MQ=255)
                                       cTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTgg                                                                                                    <  2:1812883/150‑1 (MQ=255)
                                            cGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcg                                                                                                                    >  2:127756/1‑129 (MQ=255)
                                            cGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcg                                                                                                                    <  1:127756/129‑1 (MQ=255)
                                              aTTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGcc                                                                                             <  1:481946/150‑1 (MQ=255)
                                              aTTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGcc                                                                                             <  1:1956340/150‑1 (MQ=255)
                                                  ggCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCagta                                                                                         >  2:327426/1‑150 (MQ=255)
                                                  ggCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCagta                                                                                         >  2:1114910/1‑150 (MQ=255)
                                                      aaGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGc                                                                                     <  2:297794/150‑1 (MQ=255)
                                                           tGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATAc                                                                                <  1:1855187/150‑1 (MQ=255)
                                                            ggATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACg                                                                               >  1:1944830/1‑150 (MQ=255)
                                                                cTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGAttt                                                                           >  1:869725/1‑150 (MQ=255)
                                                                  aCCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATTTTTTTAAAAACGGCTATCTTACGAGTCCATAACGACTGTATATTCTAGTAGTTCAATACGATTTCa                                                                         <  2:509914/150‑1 (MQ=255)
                                                                        ttCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGTGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGc                                                                                     >  1:718391/1‑132 (MQ=255)
                                                                        ttCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGTGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGc                                                                                     <  2:711320/132‑1 (MQ=255)
                                                                        ttCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGTGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGc                                                                                     >  1:711320/1‑132 (MQ=255)
                                                                        ttCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGTGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGc                                                                                     <  2:718391/132‑1 (MQ=255)
                                                                              aTCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAgg                                                             >  2:1227571/1‑150 (MQ=255)
                                                                                cACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGca                                                           >  1:1023248/1‑150 (MQ=255)
                                                                                cACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGca                                                           >  1:82783/1‑150 (MQ=255)
                                                                                   ggCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACag                                                        <  2:486745/150‑1 (MQ=255)
                                                                                              cggcgATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCt                                             >  1:1409746/1‑150 (MQ=255)
                                                                                                   aTGAAAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCa                                        <  1:1105582/150‑1 (MQ=255)
                                                                                                    tGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAg                                       >  1:1865082/1‑150 (MQ=255)
                                                                                                      agaAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTg                                            >  1:602166/1‑143 (MQ=255)
                                                                                                      agaAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTg                                            <  2:602166/143‑1 (MQ=255)
                                                                                                      agaAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCg                                     >  1:1095528/1‑150 (MQ=255)
                                                                                                          aGCCCCTTTCCGTTCATTCATTTATTCTCTGCCGGTTCGTTATATTATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGAGaata                                 <  1:88093/150‑3 (MQ=255)
                                                                                                            cccTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCg                               <  2:795486/150‑1 (MQ=255)
                                                                                                               ttttCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGAtc                            <  2:1643304/150‑1 (MQ=255)
                                                                                                                        aCTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTg                   >  1:1847856/1‑150 (MQ=255)
                                                                                                                                    aGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaac       >  1:233781/1‑150 (MQ=255)
                                                                                                                                     gCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGTCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaaca      <  2:215832/150‑1 (MQ=255)
                                                                                                                                     gCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaaca      >  1:440281/1‑150 (MQ=255)
                                                                                                                                     gCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCACTAGTGCAATACGATTTCAGCCGCTTGCATGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaaca      >  2:402419/1‑150 (MQ=255)
                                                                                                                                        gCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaacaaca   <  2:1409746/150‑1 (MQ=255)
                                                                                                                                         ccGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAaacaacaa  <  1:1557599/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGATATGCCAGTAGTGCAATACGATTTCAGCCGCTTGCAGGCACAGAGTAACGGGCTGTTCAGCGAAGCCGATCTCAACCGTGAGCTGAAAAAACAACAA  >  NC_000913/2281050‑2281336

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