Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,289,657 T→C L282L (TTG→CTG)  yraJ → putative fimbrial usher protein YraJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,6570TC96.3% 175.5 / NA 55L282L (TTG→CTG) yraJputative fimbrial usher protein YraJ
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (20/32);  minor base G (2/0);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTAGTTATTTACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACG  >  NC_000913/3289529‑3289789
                                                                                                                                |                                                                                                                                    
gtagtTATTTACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAAcagc                                                                                                                           <  1:512558/140‑1 (MQ=255)
gtagtTATTTACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAAcagc                                                                                                                           >  2:20482/1‑140 (MQ=255)
    ttATTTACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTCCTGGCGTGCAACACGCAGCCGACGAGAATAACCTGCCAAAcagccagc                                                                                                                       >  2:309897/1‑140 (MQ=255)
     tATTTACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAgcg                                                                                                                      <  1:408237/140‑1 (MQ=255)
       tttACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATCCTGCCAAACAGCCAGcgcg                                                                                                                    >  2:415855/1‑140 (MQ=255)
           cAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGAtt                                                                                                                >  2:641224/1‑140 (MQ=255)
           cAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGAtt                                                                                                                >  2:404110/1‑140 (MQ=255)
              cGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTTCGAGGATAATATGCTGTCAAAGAGCCAGCGCGGATTTGc                                                                                                             >  2:186623/1‑140 (MQ=255)
              cGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGc                                                                                                             >  2:428328/1‑140 (MQ=255)
              cGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGc                                                                                                             >  2:353130/1‑140 (MQ=255)
              cGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCAACGATAATATGCTGGCAAACAACCAGCGCGGATTTGc                                                                                                             >  1:194367/1‑140 (MQ=255)
               gTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGcc                                                                                                            >  2:37297/1‑140 (MQ=255)
                tGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGccc                                                                                                           <  1:353130/140‑1 (MQ=255)
                     tCAAGGCGTTGAAGTCTCAGTTACTCCTGTGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACg                                                                                                      <  1:177107/140‑1 (MQ=255)
                        aGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCGGACGATAAGATGGTGCGAAACAGCCGGCGCTGATTTGCCCCAACGGTa                                                                                                   >  2:177107/1‑140 (MQ=255)
                        aGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCGGACGATAAAATGGTGCGAAACAGCCCGCGCGGATTTGCCCCAACGGTa                                                                                                   >  2:216669/1‑140 (MQ=255)
                        aGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTa                                                                                                   >  2:152194/1‑140 (MQ=255)
                          gCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTAcg                                                                                                 <  2:632999/140‑1 (MQ=255)
                          gCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTAcg                                                                                                 <  2:73911/140‑1 (MQ=255)
                             ttGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCgg                                                                                              >  2:578064/1‑140 (MQ=255)
                             ttGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCgg                                                                                              >  1:441969/1‑140 (MQ=255)
                               gAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTa                                                                                            <  1:361604/140‑1 (MQ=255)
                                    ctcAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCa                                                                                       <  2:12304/140‑1 (MQ=255)
                                             ttCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCaaa                                                                                     >  2:490059/1‑133 (MQ=255)
                                             ttCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCaaa                                                                                     <  1:490059/133‑1 (MQ=255)
                                             ttCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGCATTTGCCACAAGGGAACGCGGCATCGCAAACagtagt                                                                              >  1:188768/1‑140 (MQ=255)
                                              tCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTg                                                                             >  2:563468/1‑140 (MQ=255)
                                              tCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTg                                                                             <  2:161883/140‑1 (MQ=255)
                                              tCTGGGAGAAAGCGCCACCAGCCGCCGTACTTTTTCCCGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTg                                                                             <  1:190198/140‑1 (MQ=255)
                                                       aaGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCGCAATGGTGCGCGGTATAGCAATCAGTAGTGCAATCGTGa                                                                    >  2:193502/1‑140 (MQ=255)
                                                          cgcCACCAGCGGCAGTATTTTTACCCGCTTCACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTa                                                                 <  1:275116/140‑1 (MQ=255)
                                                            ccaccaGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATc                                                               <  1:404110/140‑1 (MQ=255)
                                                            ccaccaGCGGCAGTATTTTTTCCAGCTACACCTTTACTGCCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATc                                                               <  2:261941/140‑1 (MQ=255)
                                                               ccaGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAgg                                                            <  1:563468/140‑1 (MQ=255)
                                                                 aGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCa                                                          <  2:631718/140‑1 (MQ=255)
                                                                  gcagCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCaa                                                         <  2:119407/137‑1 (MQ=255)
                                                                  gCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTTCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCaa                                                         <  1:186623/140‑1 (MQ=255)
                                                                  gCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCaa                                                         <  1:428328/140‑1 (MQ=255)
                                                                  gCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCaa                                                         <  1:490552/140‑1 (MQ=255)
                                                                  gCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCaa                                                         <  1:37297/140‑1 (MQ=255)
                                                                    ggCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCaaaa                                                       <  1:578064/140‑1 (MQ=255)
                                                                         tATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGtt                                                  <  1:422491/140‑1 (MQ=255)
                                                                                       acacCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCaa                                    <  1:454254/140‑1 (MQ=255)
                                                                                           cTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCa                                <  2:441969/140‑1 (MQ=255)
                                                                                             ttACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAAc                              <  1:103862/140‑1 (MQ=255)
                                                                                              tACTGGCGTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACg                             >  2:385993/1‑140 (MQ=255)
                                                                                                     gTGCAACTCGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGc                      >  2:138177/1‑140 (MQ=255)
                                                                                                        cAACTCGCTTCCGACGACAATATGCTGCCAACCAGCCAGCGCGGATTTACCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCggg                   <  2:194367/140‑1 (MQ=255)
                                                                                                             cGCTTCCGACGATAATATGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCCGGCAAAATGGTTGTGTGATCTATCAAAGCACCGTGTCAGCGGGTGCCt              >  2:269276/1‑140 (MQ=255)
                                                                                                                      cgaTAATATCCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTa     <  1:415855/140‑1 (MQ=255)
                                                                                                                        ataataTGCTGCCAAACCGCCACCGCGGATTCGCCGCAACGGTACGCGTTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAAc   <  2:174513/140‑1 (MQ=255)
                                                                                                                         taataTGTTGCCAAACAGCCAGTGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACg  <  1:303554/140‑1 (MQ=255)
                                                                                                                         taataTGCTGCCAAATAGCCAGCGCAGATGTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACg  <  1:176300/140‑1 (MQ=255)
                                                                                                                         taataTGCTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACg  <  1:20482/140‑1 (MQ=255)
                                                                                                                         taataTGCTGCCAAACAGCCAGCGCGGATATGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACg  <  1:254792/140‑1 (MQ=255)
                                                                                                                                |                                                                                                                                    
GTAGTTATTTACAACGTGATATCAAGGCGTTGAAGTCTCAGTTGCTTCTGGGAGAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACG  >  NC_000913/3289529‑3289789

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