Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,426,288 G→A intergenic (+206/‑166) ydaY → / → tmpR pseudogene, Rac prophage,Phage or Prophage Related/Rac prophage, pseudogene, tail protein family,Phage or Prophage Related, putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,2880GA97.9% 147.6 / ‑4.0 48intergenic (+206/‑166)ydaY/tmpRpseudogene, Rac prophage,Phage or Prophage Related/Rac prophage, pseudogene, tail protein family,Phage or Prophage Related, putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/24);  total (23/25)
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.74e-01

TGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACT  >  NC_000913/1426153‑1426417
                                                                                                                                       |                                                                                                                                 
tGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAattg                                                                                                                                >  2:542725/1‑139 (MQ=255)
 ggggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAattgt                                                                                                                               <  1:641781/139‑1 (MQ=255)
   ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                             >  2:80858/1‑139 (MQ=255)
   ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                             >  2:197491/1‑139 (MQ=255)
             cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTaca                                                                                                                   <  1:707206/139‑1 (MQ=255)
               tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                 <  1:80858/139‑1 (MQ=255)
               tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                 <  1:197491/139‑1 (MQ=255)
                 tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACActct                                                                                                               >  2:730848/1‑139 (MQ=255)
                     tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                                                                                           <  1:594898/139‑1 (MQ=255)
                        gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                        <  1:296498/139‑1 (MQ=255)
                        gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                        <  1:181209/139‑1 (MQ=255)
                        gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                        >  2:409616/1‑139 (MQ=255)
                        gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                        <  1:264329/139‑1 (MQ=255)
                             agaTAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                        >  2:754442/1‑134 (MQ=255)
                             agaTAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                        <  1:754442/134‑1 (MQ=255)
                               aTAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGAtgtg                                                                                                 >  1:415044/1‑139 (MQ=255)
                                         atTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCtt                                                                                       >  2:699527/1‑139 (MQ=255)
                                         atTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCtt                                                                                       >  1:786373/1‑139 (MQ=255)
                                             aGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGAtgtg                                                                                                 <  1:585361/125‑1 (MQ=255)
                                             aGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGAtgtg                                                                                                 >  2:585361/1‑125 (MQ=255)
                                             aGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaat                                                                                   >  2:367385/1‑139 (MQ=255)
                                                         ggTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatgat                                                                       <  2:366445/139‑1 (MQ=255)
                                                                    aGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTaca                                                                                                                   >  1:51403/1‑84 (MQ=255)
                                                                    aGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTaca                                                                                                                   <  2:51403/84‑1 (MQ=255)
                                                                     gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaa                                                           >  2:574402/1‑139 (MQ=255)
                                                                     gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaa                                                           >  1:10588/1‑139 (MQ=255)
                                                                     gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaa                                                           >  1:591262/1‑139 (MQ=255)
                                                                     gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGAATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaa                                                           >  1:424757/1‑139 (MQ=255)
                                                                       gcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaaaa                                                         >  2:24362/1‑139 (MQ=255)
                                                                          aTGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                                                                                              >  1:73112/1‑99 (MQ=255)
                                                                          aTGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                                                                                              <  2:73112/99‑1 (MQ=255)
                                                                               tttGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCg                                                 <  1:730848/139‑1 (MQ=255)
                                                                                          gATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTaa                                      <  2:591262/139‑1 (MQ=255)
                                                                                            tCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATc                                    <  2:415044/139‑1 (MQ=255)
                                                                                             cTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCg                                   <  2:583979/139‑1 (MQ=255)
                                                                                               ttGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                 <  1:574402/139‑1 (MQ=255)
                                                                                               ttGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                 <  1:699527/139‑1 (MQ=255)
                                                                                               ttGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                 <  2:786373/139‑1 (MQ=255)
                                                                                               ttGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                 <  2:10588/139‑1 (MQ=255)
                                                                                                   aCTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTatat                             >  1:667033/1‑139 (MQ=255)
                                                                                                    cTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                            <  1:409616/139‑1 (MQ=255)
                                                                                                     tGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATct                           >  2:705963/1‑139 (MQ=255)
                                                                                                       gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTa                                                                                                                     >  2:124749/1‑47 (MQ=255)
                                                                                                       gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTa                                                                                                                     <  1:124749/47‑1 (MQ=255)
                                                                                                               gTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGt                 <  2:667033/139‑1 (MQ=255)
                                                                                                                            cAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                 <  2:331848/110‑1 (MQ=255)
                                                                                                                            cAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                 >  1:331848/1‑110 (MQ=255)
                                                                                                                              gtgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACt  >  1:103474/1‑139 (MQ=255)
                                                                                                                               tgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                            >  2:603585/1‑112 (MQ=255)
                                                                                                                               tgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                            <  1:603585/112‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
TGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACT  >  NC_000913/1426153‑1426417

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