Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,620,414 A→C 6.1% K408T (AAG→ACG) ‡ rhsB → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,620,4140AC 6.1% 223.9 / 4.2 63K408T (AAG→ACG) ‡rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base C (2/2);  total (31/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 = 9.99e-01

TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCA  >  NC_000913/3620272‑3620551
                                                                                                                                              |                                                                                                                                         
tACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     >  2:1486987‑M2/1‑149 (MQ=255)
 aCGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    >  2:1350118‑M2/1‑149 (MQ=255)
  cGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   >  2:828302‑M2/1‑149 (MQ=255)
   gACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  >  1:718922‑M2/1‑149 (MQ=255)
    aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 <  1:28019‑M2/149‑1 (MQ=255)
    aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 <  1:1535035‑M2/149‑1 (MQ=255)
    aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 >  1:1820414‑M2/1‑149 (MQ=255)
    aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCAGTGAAGTACTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGGGGTG                                                                                                                                 <  4:1595313‑M2/149‑1 (MQ=255)
     cAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGA                                                                                                                                <  4:1425553‑M2/149‑1 (MQ=255)
        cGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    <  1:1093183‑M2/142‑1 (MQ=255)
        cGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    >  2:1115574‑M2/1‑142 (MQ=255)
        cGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    >  2:1093183‑M2/1‑142 (MQ=255)
        cGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    <  1:1115574‑M2/142‑1 (MQ=255)
              gcgggTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACGCGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  >  2:276802‑M2/1‑138 (MQ=255)
              gcgggTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  <  1:276802‑M2/138‑1 (MQ=255)
                 ggTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 <  2:24368‑M2/136‑1 (MQ=255)
                 ggTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 >  1:24368‑M2/1‑136 (MQ=255)
                  gTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCGCACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    >  2:1692961‑M2/1‑132 (MQ=255)
                  gTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCGCACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    <  1:1692961‑M2/132‑1 (MQ=255)
                   tGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     >  3:171802‑M2/1‑130 (MQ=255)
                   tGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     <  4:171802‑M2/130‑1 (MQ=255)
                              cTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    <  4:308021‑M2/120‑1 (MQ=255)
                              cTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    >  3:308021‑M2/1‑120 (MQ=255)
                               tAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   >  2:1217964‑M2/1‑120 (MQ=255)
                               tAAACCCGGCAGGCTTAAGCTACACGTATCAGGATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGAGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   <  1:1217964‑M2/120‑1 (MQ=255)
                                 aaCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  >  4:801236‑M2/1‑119 (MQ=255)
                                 aaCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  <  3:801236‑M2/119‑1 (MQ=255)
                                  aCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   >  3:1516043‑M2/1‑117 (MQ=255)
                                  aCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   <  4:1516043‑M2/117‑1 (MQ=255)
                                    ccGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     >  4:371436‑M2/1‑113 (MQ=255)
                                    ccGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     <  3:371436‑M2/113‑1 (MQ=255)
                                         aGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   >  2:1067755‑M2/1‑110 (MQ=255)
                                         aGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   <  1:1067755‑M2/110‑1 (MQ=255)
                                             ttAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    <  2:97320‑M2/105‑1 (MQ=255)
                                             ttAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                                                                                    >  1:97320‑M2/1‑105 (MQ=255)
                                                   tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     <  3:1551723‑M2/98‑1 (MQ=255)
                                                   tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     >  4:1551723‑M2/1‑98 (MQ=255)
                                                   tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     <  1:1031647‑M2/98‑1 (MQ=255)
                                                   tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     >  2:1031647‑M2/1‑98 (MQ=255)
                                                    acacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 >  3:1669288‑M2/1‑101 (MQ=255)
                                                    acacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                                                                                                                                 <  4:1669288‑M2/101‑1 (MQ=255)
                                                      acGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   <  1:233920‑M2/97‑1 (MQ=255)
                                                      acGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                                                                                   >  2:233920‑M2/1‑97 (MQ=255)
                                                        gTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  >  4:801970‑M2/1‑96 (MQ=255)
                                                        gTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  <  3:801970‑M2/96‑1 (MQ=255)
                                                                        gACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgccc                                                             <  3:1303230‑M2/149‑68 (MQ=255)
                                                                              atcaccatcaccGACAGCCTGGACCGCGGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgct                                                       <  2:358993‑M2/149‑74 (MQ=255)
                                                                              atcaccatcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgct                                                       <  1:1434590‑M2/149‑74 (MQ=255)
                                                                                caccatcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctctctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgctgt                                                     <  2:710622‑M2/149‑76 (MQ=255)
                                                                                  ccatcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     >  1:1575894‑M2/1‑67 (MQ=255)
                                                                                  ccatcaccGACAGCCTGGACCGCCGGGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                                                                                                                                     <  2:1575894‑M2/67‑1 (MQ=255)
                                                                                   catcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgctgtcgt                                                  <  1:653212‑M2/149‑79 (MQ=255)
                                                                                   catcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgt                                                                           >  4:862898‑M2/1‑71 (MQ=255)
                                                                                   catcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgt                                                                           <  3:862898‑M2/124‑54 (MQ=255)
                                                                                     tcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgctgtcgtaa                                                >  2:653212‑M2/1‑69 (MQ=255)
                                                                                            cAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgctgtcgtaacggtaac                                         >  4:1303230‑M2/1‑62 (MQ=255)
                                                                                                cTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGTCGGTGATGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagct                                                                             <  3:489648‑M2/109‑52 (MQ=255)
                                                                                                cTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGTCGGTGATGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagct                                                                             >  4:489648‑M2/1‑58 (MQ=255)
                                                                                                cTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  <  2:1059373‑M2/56‑1 (MQ=255)
                                                                                                cTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                                                                                  >  1:1059373‑M2/1‑56 (MQ=255)
                                                                                                  ggACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGTCGGTGATGGTGAtgcggtctttctca                                                                                                                  <  4:1018885‑M2/70‑15 (MQ=255)
                                                                                                  ggACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGTCGGTGATGGTGAtgcggtctttctca                                                                                                                  >  3:1018885‑M2/1‑56 (MQ=255)
                                                                                                                                   aGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgctgtcgtaacggtaacggatttccggccggcccgtgtgacggtgcgccaccatcc  >  1:710622‑M2/1‑23 (MQ=255)
                                                                                                                                              |                                                                                                                                         
TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCA  >  NC_000913/3620272‑3620551

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