Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,623,585 A→G 100% E63G (GAA→GGA)  yhhH → putative NTF2 fold immunity protein for polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,5850AG98.1% 168.7 / ‑3.7 52coding (188/384 nt)yhhHputative NTF2 fold immunity protein for polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/29);  minor base . (0/1);  total (22/30)
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.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATT  >  NC_000913/3623439‑3623727
                                                                                                                                                   |                                                                                                                                              
tAGGAACTGTGTCATTGTATCTTTGTTTGTTCTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGAt                                                                                                                                              <  1:150449/149‑1 (MQ=255)
   gAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAg                                                                                                                                           >  2:571607/1‑149 (MQ=255)
     acggtgtCATTGTATGTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCt                                                                                                                                         <  2:361131/146‑1 (MQ=255)
       tgtgTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTga                                                                                                                                       <  1:619243/149‑1 (MQ=255)
             aTTGTATCTTTGTTTGTTTTTTTTTATAATTCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAATTATACCGTTAGAAATTTGGTTTAAAATAATGCTATAGTATTGAGATTAGCTGAGatata                                                                                                                                 >  2:625622/1‑149 (MQ=255)
                  aTCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTa                                                                                                                            <  1:427177/149‑1 (MQ=255)
                      ttgtttgtTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGaa                                                                                                                        <  2:288809/149‑1 (MQ=255)
                        gtttgtttTTACATACAACACATGGGTGCAGTGTGATGGTAATATTAAAATGATGCGCGAGGATGAGAGGGAAGGAAAATATACCGTTAGAAATTGTGTTAGAAATAAAGCTATAGGATTGGGATTAGCGGGGAGATATGTTAAGAAtt                                                                                                                      >  2:467319/1‑148 (MQ=255)
                        gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATc                                                                                                                      <  2:277790/149‑1 (MQ=255)
                        gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATc                                                                                                                      <  1:715233/149‑1 (MQ=255)
                        gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTAAAGAATc                                                                                                                      >  2:256306/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCg                                                                                                                     >  1:142434/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATAAGCTGAGATATATGTTAAGAATCg                                                                                                                     >  1:721189/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAAAAATTCATGGTTATGAGAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTAAGATATATGTTAAGAATCg                                                                                                                     >  1:424127/1‑149 (MQ=255)
                          ttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGt                                                                                                                    <  2:648876/149‑1 (MQ=255)
                                  acatacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACa                                                                                                            <  2:188274/149‑1 (MQ=255)
                                  acatacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCAATAGAATTGGGAGTCACTGAGAAATAGGTTAAGAATCGTTATGGACa                                                                                                            >  2:304961/1‑149 (MQ=255)
                                   catacaACACATGGGCGCAGTGTAATACTAATATTAAAATAATGCGCAAGTATGAAAGTGAACGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAg                                                                                                           <  1:542552/149‑1 (MQ=255)
                                      acaacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGCTTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                                                                        >  2:517089/1‑149 (MQ=255)
                                      acaacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                                                                        <  1:713192/149‑1 (MQ=255)
                                            caTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCa                                                                                                  <  1:780691/149‑1 (MQ=255)
                                            caTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCa                                                                                                  <  1:780708/149‑1 (MQ=255)
                                              tGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCATTATGGACAGGATGCCGCaga                                                                                                <  2:683560/149‑1 (MQ=255)
                                                         taataataaTATTAAAATAATGCGCAAGTATGAAAGTGAATGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAc                                                                                     <  1:90520/149‑1 (MQ=255)
                                                         taataataaGTTTACAATAAGGCGCAAATATGAAAGTGAAGGTAAATATACTGTTAGAAATTTGGTTAAAAATAAAGCTATAACATGGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGCTGCCGCAGAAGAAGAAAAAc                                                                                     <  1:751142/149‑1 (MQ=255)
                                                             aataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAAATTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATa                                                                                 >  2:514389/1‑149 (MQ=255)
                                                                    ttAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGCAGAAAAAACATACGAAATc                                                                          >  2:750055/1‑149 (MQ=255)
                                                                       aaaTAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGCTAGAAACTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACt                                                                       <  1:544100/149‑1 (MQ=255)
                                                                         ataatGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGa                                                                     <  1:625622/149‑1 (MQ=255)
                                                                            atGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGtt                                                                  >  1:761921/1‑149 (MQ=255)
                                                                             tGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTa                                                                 <  1:240651/149‑1 (MQ=255)
                                                                               cgcAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGCAATCACTGAGTTaac                                                               >  2:97249/1‑149 (MQ=255)
                                                                                  aaGTATGAAAGGGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaacaac                                                            <  1:112402/149‑1 (MQ=255)
                                                                                   aGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGGTAacaaca                                                           >  1:537218/1‑149 (MQ=255)
                                                                                         aaaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAACACCATACGAAATCACTGAGTTAACAACAAGCTgg                                                     >  1:655895/1‑149 (MQ=255)
                                                                                         aaaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTgg                                                     <  1:644763/149‑1 (MQ=255)
                                                                                           aGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgt                                                   >  2:630360/1‑149 (MQ=255)
                                                                                           aGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATACCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTATCAACAAGTTGGgt                                                   <  2:486583/149‑1 (MQ=255)
                                                                                                  gTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGaa                                            <  2:421635/149‑1 (MQ=255)
                                                                                                   tAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAg                                           <  1:441368/149‑1 (MQ=255)
                                                                                                   tAAATATACCGTTAGAAATGTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAg                                           <  2:394301/149‑1 (MQ=255)
                                                                                                    aaaTATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGCTGGGTGGTTGAggg                                          >  2:20961/1‑149 (MQ=255)
                                                                                                     aaTATACCGTTAGGAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTGGTTGAAGGt                                         >  2:555896/1‑149 (MQ=255)
                                                                                                                  aaaTTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGAGAAATGTAAAGAATCGTGATTGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTGGTTGAAGGTACCATTCACTCAg                            >  2:171279/1‑149 (MQ=255)
                                                                                                                       tGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCGGACCaa                       >  2:45794/1‑149 (MQ=255)
                                                                                                                         gTTAAAAATAAAGCTATAGCCATTGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAt                     <  2:424127/149‑1 (MQ=255)
                                                                                                                         gTTAAAAAGAATGGTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAACCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAt                     <  2:142434/149‑1 (MQ=255)
                                                                                                                              aaaTAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtg                <  1:325012/149‑1 (MQ=255)
                                                                                                                              aaaTAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGAACAAATAGCtg                >  2:494609/1‑149 (MQ=255)
                                                                                                                                  aaaGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggtgg            <  2:761921/149‑1 (MQ=255)
                                                                                                                                       tataGCATCGGGATAAGCGGAGATATAAGATAAGAATCGATATCGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACAAATCACTCAGACCAAAATGCTGGTgcggaat       >  2:394568/1‑146 (MQ=255)
                                                                                                                                             cATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTttatt  <  1:571607/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATT  >  NC_000913/3623439‑3623727

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