Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,002,224 (G)7→8 coding (418/1071 nt) elfG → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,002,2171.G93.5% 100.6 / 0.4 31coding (411/1071 nt)elfGputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (16/13);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TGGGCATTCGCCGCGGATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGGTCGTC  >  NC_000913/1002080‑1002340
                                                                                                                                          |                                                                                                                           
tGGGCATTCGCCGCGGATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAg                                                                                                                            >  2:1306768/1‑139 (MQ=255)
                                     gCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGa                                                                                                              <  1:1003709/117‑1 (MQ=255)
                                     gCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGa                                                                                                              >  2:1003709/1‑117 (MQ=255)
                                       gACAACTGGTACAAAATTAATGACTATCTTGCCGCAAAAACTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGCCACAGGGGGGGGATACCCGATGGCATATTTGCGACCATGGCGGt                                                                                      >  2:1113717/1‑139 (MQ=255)
                                          aaCTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCaa                                                                                   >  1:495566/1‑139 (MQ=255)
                                               gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCa                                                                                                         <  1:590309/112‑1 (MQ=255)
                                               gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCa                                                                                                         >  2:590309/1‑112 (MQ=255)
                                                 aCAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAg                                                                            <  2:615262/139‑1 (MQ=255)
                                                   aaaaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGt                                                                          >  1:499752/1‑139 (MQ=255)
                                                          aTGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAggg                                                                  >  2:1057604/1‑139 (MQ=255)
                                                                  cTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCg                                                          <  1:1057604/139‑1 (MQ=255)
                                                                         ccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGa                                                                                                              <  1:715454/81‑1 (MQ=255)
                                                                         ccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGa                                                                                                              >  2:715454/1‑81 (MQ=255)
                                                                           aaaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAAt                                                     >  2:683682/1‑136 (MQ=255)
                                                                           aaaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAAt                                                     <  1:683682/136‑1 (MQ=255)
                                                                                         ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGt                                                                          >  2:653579/1‑101 (MQ=255)
                                                                                         ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGt                                                                          <  1:653579/101‑1 (MQ=255)
                                                                                                     ttCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATAcc                                                                                                                 >  2:960408/1‑50 (MQ=255)
                                                                                                     ttCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATAcc                                                                                                                 <  1:960408/50‑1 (MQ=255)
                                                                                                               gCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGa                                                                                                              >  2:193440/1‑43 (MQ=255)
                                                                                                               gCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGa                                                                                                              <  1:193440/43‑1 (MQ=255)
                                                                                                                  gTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTccc                                     >  1:930338/1‑113 (MQ=255)
                                                                                                                  gTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTccc                                     <  2:930338/113‑1 (MQ=255)
                                                                                                                  gTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTTGTTGGc           >  1:581663/1‑139 (MQ=255)
                                                                                                                  gTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGc           >  2:130612/1‑139 (MQ=255)
                                                                                                                  gTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGc           >  2:280524/1‑139 (MQ=255)
                                                                                                                     ccTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCg        <  2:499752/139‑1 (MQ=255)
                                                                                                                           gAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGgtcgtc  <  2:1091947/139‑1 (MQ=255)
                                                                                                                                   gggCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCttt                                         <  1:683816/92‑1 (MQ=255)
                                                                                                                                   gggCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCttt                                         >  2:683816/1‑92 (MQ=255)
                                                                                                                                    ggCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAg                                                                            >  1:697042/1‑56 (MQ=255)
                                                                                                                                    ggCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAg                                                                            <  2:697042/56‑1 (MQ=255)
                                                                                                                                          |                                                                                                                           
TGGGCATTCGCCGCGGATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGGTCGTC  >  NC_000913/1002080‑1002340

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: