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.

TAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATT  >  NC_000913/3623439‑3623727
                                                                                                                                                   |                                                                                                                                              
tAGGAACTGTGTCATTGTATCTTTGTTTGTTCTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGAt                                                                                                                                              <  1:150449/149‑1 (MQ=255)
   gAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAg                                                                                                                                           >  2:571607/1‑149 (MQ=255)
     acggtgtCATTGTATGTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCt                                                                                                                                         <  2:361131/146‑1 (MQ=255)
       tgtgTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTga                                                                                                                                       <  1:619243/149‑1 (MQ=255)
             aTTGTATCTTTGTTTGTTTTTTTTTATAATTCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAATTATACCGTTAGAAATTTGGTTTAAAATAATGCTATAG‑TATTGAGATTAGCTGAGatata                                                                                                                                 >  2:625622/1‑149 (MQ=255)
                  aTCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTa                                                                                                                            <  1:427177/149‑1 (MQ=255)
                      ttgtttgtTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGaa                                                                                                                        <  2:288809/149‑1 (MQ=255)
                        gtttgtttTTACATACAACACATGGGTGCAGTGTGATGGTAATATTAAAATGATGCGCGAGGATGAGAGGGAAGGAAAATATACCGTTAGAAATTGTGTTAGAAATAAAGCTATAG‑GATTGGGATTAGCGGGGAGATATGTTAAGAAtt                                                                                                                      >  2:467319/1‑148 (MQ=255)
                        gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATc                                                                                                                      <  2:277790/149‑1 (MQ=255)
                        gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATc                                                                                                                      <  1:715233/149‑1 (MQ=255)
                        gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTAAAGAATc                                                                                                                      >  2:256306/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCg                                                                                                                     >  1:142434/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATAAGCTGAGATATATGTTAAGAATCg                                                                                                                     >  1:721189/1‑149 (MQ=255)
                         tttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAAAAATTCATGGTTATGAGAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTAAGATATATGTTAAGAATCg                                                                                                                     >  1:424127/1‑149 (MQ=255)
                          ttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGt                                                                                                                    <  2:648876/149‑1 (MQ=255)
                                  acatacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACa                                                                                                            <  2:188274/149‑1 (MQ=255)
                                  acatacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCAATAG‑AATTGGGAGTCACTGAGAAATAGGTTAAGAATCGTTATGGACa                                                                                                            >  2:304961/1‑149 (MQ=255)
                                   catacaACACATGGGCGCAGTGTAATACTAATATTAAAATAATGCGCAAGTATGAAAGTGAACGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAg                                                                                                           <  1:542552/149‑1 (MQ=255)
                                      acaacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGCTTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                                                                        >  2:517089/1‑149 (MQ=255)
                                      acaacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                                                                        <  1:713192/149‑1 (MQ=255)
                                            caTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCa                                                                                                  <  1:780691/149‑1 (MQ=255)
                                            caTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCa                                                                                                  <  1:780708/149‑1 (MQ=255)
                                              tGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCATTATGGACAGGATGCCGCaga                                                                                                <  2:683560/149‑1 (MQ=255)
                                                         taataataaTATTAAAATAATGCGCAAGTATGAAAGTGAATGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAc                                                                                     <  1:90520/149‑1 (MQ=255)
                                                         taataataaGTTTACAATAAGGCGCAAATATGAAAGTGAAGGTAAATATACTGTTAGAAATTTGGTTAAAAATAAAGCTATAA‑CATGGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGCTGCCGCAGAAGAAGAAAAAc                                                                                     <  1:751142/149‑1 (MQ=255)
                                                             aataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAAATTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATa                                                                                 >  2:514389/1‑149 (MQ=255)
                                                                    ttAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGCAGAAAAAACATACGAAATc                                                                          >  2:750055/1‑149 (MQ=255)
                                                                       aaaTAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGCTAGAAACTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACt                                                                       <  1:544100/149‑1 (MQ=255)
                                                                         ataatGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGa                                                                     <  1:625622/149‑1 (MQ=255)
                                                                            atGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGtt                                                                  >  1:761921/1‑149 (MQ=255)
                                                                             tGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTa                                                                 <  1:240651/149‑1 (MQ=255)
                                                                               cgcAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGCAATCACTGAGTTaac                                                               >  2:97249/1‑149 (MQ=255)
                                                                                  aaGTATGAAAGGGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaacaac                                                            <  1:112402/149‑1 (MQ=255)
                                                                                   aGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGGTAacaaca                                                           >  1:537218/1‑149 (MQ=255)
                                                                                         aaaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAACACCATACGAAATCACTGAGTTAACAACAAGCTgg                                                     >  1:655895/1‑149 (MQ=255)
                                                                                         aaaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTgg                                                     <  1:644763/149‑1 (MQ=255)
                                                                                           aGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgt                                                   >  2:630360/1‑149 (MQ=255)
                                                                                           aGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAC‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTATCAACAAGTTGGgt                                                   <  2:486583/149‑1 (MQ=255)
                                                                                                  gTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGaa                                            <  2:421635/149‑1 (MQ=255)
                                                                                                   tAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAg                                           <  1:441368/149‑1 (MQ=255)
                                                                                                   tAAATATACCGTTAGAAATGTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAg                                           <  2:394301/149‑1 (MQ=255)
                                                                                                    aaaTATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGCTGGGTGGTTGAggg                                          >  2:20961/1‑149 (MQ=255)
                                                                                                     aaTATACCGTTAGGAATTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTGGTTGAAGGt                                         >  2:555896/1‑149 (MQ=255)
                                                                                                                  aaaTTTGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGAGAAATGTAAAGAATCGTGATTGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTGGTTGAAGGTACCATTCACTCAg                            >  2:171279/1‑149 (MQ=255)
                                                                                                                       tGGTTAAAAATAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCGGACCaa                       >  2:45794/1‑149 (MQ=255)
                                                                                                                         gTTAAAAATAAAGCTATAGCCATTGG‑ATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAt                     <  2:424127/149‑1 (MQ=255)
                                                                                                                         gTTAAAAAGAATGGTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAACCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAt                     <  2:142434/149‑1 (MQ=255)
                                                                                                                              aaaTAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtg                <  1:325012/149‑1 (MQ=255)
                                                                                                                              aaaTAAAGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGAACAAATAGCtg                >  2:494609/1‑149 (MQ=255)
                                                                                                                                  aaaGCTATAG‑CATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggtgg            <  2:761921/149‑1 (MQ=255)
                                                                                                                                       tataG‑CATCGGGATAAGCGGAGATATAAGATAAGAATCGATATCGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACAAATCACTCAGACCAAAATGCTGGTgcggaat       >  2:394568/1‑146 (MQ=255)
                                                                                                                                             cATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTttatt  <  1:571607/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAG‑CATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATT  >  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: ‑