Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 150,947 T→C T3A (ACT→GCT)  yadC ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913150,9470TC97.8% 155.8 / ‑6.8 45T3A (ACT→GCT) yadCputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/19);  minor base G (0/1);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGA  >  NC_000913/150832‑151073
                                                                                                                   |                                                                                                                              
gAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTc                                                                                                         >  1:1233190/1‑139 (MQ=255)
        tCGATAGTTATCAACTCCTGCATTATTTCTCACAATAAAACTTGTCTATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGt                                                                                                 >  2:1001099/1‑139 (MQ=255)
              gTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGt                                                                                           >  2:409286/1‑139 (MQ=255)
                    aaCTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGttt                                                                                                          >  1:1004141/1‑118 (MQ=255)
                    aaCTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGttt                                                                                                          <  2:1004141/118‑1 (MQ=255)
                           gCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGttttt                                                                              <  2:1233190/139‑1 (MQ=255)
                                attTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgtttg                                                                         <  2:1118263/139‑1 (MQ=255)
                                   tCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTa                                                                      >  2:990864/1‑139 (MQ=255)
                                           tAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgtagt                                                              >  2:1267136/1‑139 (MQ=255)
                                                cTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                         <  1:1267136/139‑1 (MQ=255)
                                                  tGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                       <  1:1001099/139‑1 (MQ=255)
                                                   gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATcc                                                      >  1:1094230/1‑139 (MQ=255)
                                                            ctgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGGCTTCATTTTTACGCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                             <  2:1604329/139‑1 (MQ=255)
                                                            ctgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                             >  1:1459136/1‑139 (MQ=255)
                                                             tgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACAt                                            <  2:1094230/139‑1 (MQ=255)
                                                              gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTa                                                                      <  2:528767/112‑1 (MQ=255)
                                                              gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTa                                                                      >  1:528767/1‑112 (MQ=255)
                                                              gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATg                                           >  2:404714/1‑139 (MQ=255)
                                                              gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATg                                           >  2:1389054/1‑139 (MQ=255)
                                                               actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                 <  1:1282631/116‑1 (MQ=255)
                                                               actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                 >  2:1282631/1‑116 (MQ=255)
                                                               actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATgg                                          >  2:684368/1‑139 (MQ=255)
                                                                    aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     >  2:437065/1‑75 (MQ=255)
                                                                    aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     <  1:437065/75‑1 (MQ=255)
                                                                    aTTACAACAAGAATACAGAGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     <  1:1283780/75‑1 (MQ=255)
                                                                    aTTACAACAAGAATACAGAGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     >  2:1283780/1‑75 (MQ=255)
                                                                               aaTACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                >  2:1243232/1‑69 (MQ=255)
                                                                               aaTACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                <  1:1243232/69‑1 (MQ=255)
                                                                                   cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                >  1:5910/1‑65 (MQ=255)
                                                                                   cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                <  2:5910/65‑1 (MQ=255)
                                                                                   cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAg                          <  2:1195738/135‑1 (MQ=255)
                                                                                   cAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAg                          >  1:1195738/1‑135 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTtgtttgtt                                                                       <  2:1211342/87‑1 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTtgtttgtt                                                                       >  1:1211342/1‑87 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                 <  2:1301018/93‑1 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                 >  1:1301018/1‑93 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAaca                                               >  1:739745/1‑111 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAaca                                               <  2:739745/111‑1 (MQ=255)
                                                                                      tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac                   >  2:151166/1‑139 (MQ=255)
                                                                                        cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAacac                 >  2:945452/1‑139 (MQ=255)
                                                                                         tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                <  1:151166/139‑1 (MQ=255)
                                                                                              aaaGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATg           <  1:945452/139‑1 (MQ=255)
                                                                                                       aCCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                 >  2:1045640/1‑76 (MQ=255)
                                                                                                       aCCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                 <  1:1045640/76‑1 (MQ=255)
                                                                                                       aCCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGa  >  1:288756/1‑139 (MQ=255)
                                                                                                                   |                                                                                                                              
GAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGA  >  NC_000913/150832‑151073

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