Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 732,414 G→A 7.1% R944R (AGG→AGA)  rhsC → Rhs protein with 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_000913732,4140GA 7.1% 153.2 / 2.7 56R944R (AGG→AGA) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base A (2/2);  total (24/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 = 5.84e-01

GGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTA  >  NC_000913/732267‑732557
                                                                                                                                                   |                                                                                                                                               
ggATATCCGCATCCCGTATACCACAGACCCGGCAGGTCACCGCCGGCCCGCCCCGGAGCTGCACCCGGACAGCGCCCTCAGCAGGTGGCCGGATACCCGTATCGCCCGTGACGCGCACTCTCTTTACAGGGATGCCCGTCAcggccggc                                                                                                                                                >  6:355198/1‑149 (MQ=11)
ggATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                                                                                                                                                >  8:6505/1‑149 (MQ=34)
  atatCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                                                                                                                                               >  7:250432/1‑148 (MQ=34)
       cGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACaga                                                                                                                                         <  1:232125/149‑1 (MQ=34)
       cGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACaga                                                                                                                                         >  6:295284/1‑149 (MQ=34)
        gCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACa                                                                                                                                           >  7:204085/1‑146 (MQ=34)
          aTCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaa                                                                                                                                      <  4:75174/149‑1 (MQ=34)
           tCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaa                                                                                                                                      >  4:15383/1‑148 (MQ=34)
               gTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCg                                                                                                                                 <  7:274349/149‑1 (MQ=34)
                tataCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGa                                                                                                                                <  3:202869/149‑1 (MQ=34)
                  taCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGAc                                                                                                                               <  7:6505/148‑1 (MQ=33)
                     cacaGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACa                                                                                                                                           <  5:72703/133‑1 (MQ=32)
                     cacaGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACa                                                                                                                                           >  6:72703/1‑133 (MQ=32)
                      acaGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCAt                                                                                                                          >  4:94907/1‑149 (MQ=32)
                        aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATcc                                                                                                                        <  8:122495/149‑1 (MQ=32)
                         gACGCGGCAGGTAACCGCCTGCCCGACCCGGCGCGGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATACCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAACCCGACCTCATccc                                                                                                                       <  3:7065/149‑1 (MQ=11)
                               gCAGGTAACCGCCTGCCCGACCCGGAGCTG‑ACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAggg                                                                                                                <  6:83771/149‑1 (MQ=17)
                                cacgTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGTTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAggg                                                                                                                <  4:49868/146‑1 (MQ=17)
                                 aGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAgggg                                                                                                               <  3:94907/149‑1 (MQ=32)
                                 aGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAgggg                                                                                                               <  6:35335/149‑1 (MQ=32)
                                       ccGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc                                                                                                         <  3:16150/149‑1 (MQ=32)
                                                cGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgat                                                                                                 >  1:205063/1‑148 (MQ=32)
                                                  aCCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgatgat                                                                                              <  4:2326/149‑1 (MQ=32)
                                                    ccGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAGACCGACCTCATCCCGGAAGGGGTTATCCGCACGGAtgatga                                                                                            <  6:182975/149‑1 (MQ=17)
                                                      ggAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgc                                                                                          >  3:9490/1‑149 (MQ=32)
                                                                ccGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                                                                                      >  8:308216/1‑143 (MQ=34)
                                                                ccGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                                                                                      <  7:308216/143‑1 (MQ=34)
                                                                  ggACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAgaga                                                                                                                                       <  4:265660/92‑1 (MQ=33)
                                                                         gcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTa                                                                       <  4:9490/149‑1 (MQ=32)
                                                                          cccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc                                                                      <  8:151513/149‑1 (MQ=32)
                                                                           ccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAACCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc                                                                      >  8:384437/1‑148 (MQ=17)
                                                                           ccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc                                                                      >  3:76247/1‑148 (MQ=32)
                                                                            cTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa                                                                    <  8:250432/149‑1 (MQ=32)
                                                                               aGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACa                                                                  <  2:205063/148‑1 (MQ=32)
                                                                                            aTAACCTTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGACTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACa                                                                  >  7:166768/1‑135 (MQ=17)
                                                                                            aTAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGACTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACa                                                                  <  8:166768/135‑1 (MQ=17)
                                                                                            aTAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGACTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACa                                                                  >  1:348466/1‑135 (MQ=17)
                                                                                            aTAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGACTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACa                                                                  <  2:348466/135‑1 (MQ=17)
                                                                                                        cccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCAc                                                          <  7:360749/131‑1 (MQ=32)
                                                                                                        cccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCAc                                                          >  8:360749/1‑131 (MQ=32)
                                                                                                        cccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcg                                        >  1:258501/1‑149 (MQ=32)
                                                                                                         ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcgc                                       <  5:295284/149‑1 (MQ=32)
                                                                                                                    aCTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                            <  2:258501/149‑1 (MQ=32)
                                                                                                                     cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCa                                                           <  2:94381/117‑1 (MQ=32)
                                                                                                                      tATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATg                           >  4:209911/1‑148 (MQ=32)
                                                                                                                           ttACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGCCACAATATGCAGAGc                     >  5:265647/1‑149 (MQ=17)
                                                                                                                            tACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTATGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGc                          <  4:315978/143‑1 (MQ=17)
                                                                                                                            tACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAAGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGc                          >  3:315978/1‑143 (MQ=17)
                                                                                                                            tACCGGTAGGACCGTCACGGCAGGCTGACAGAGTACACCGAGCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCaga                       <  5:355198/146‑1 (MQ=255)
                                                                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc                                                                      >  6:15922/1‑97 (MQ=33)
                                                                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc                                                                      <  5:15922/97‑1 (MQ=33)
                                                                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGc                  >  8:107517/1‑149 (MQ=32)
                                                                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGc                  <  8:204085/149‑1 (MQ=32)
                                                                                                                                 gTATGCCCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACCCAATATGCAGAGCCGCTgg               >  5:226911/1‑149 (MQ=17)
                                                                                                                                    tGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTc             <  3:15383/148‑1 (MQ=32)
                                                                                                                                     gACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGGCg            >  4:324048/1‑148 (MQ=17)
                                                                                                                                              ggCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTa  >  2:173175/1‑149 (MQ=34)
                                                                                                                                                   |                                                                                                                                               
GGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTA  >  NC_000913/732267‑732557

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