Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,426,575 G→T S41I (AGC→ATC)  ynaA → Rac prophage; putative prophage tail length tape measure domain‑containing protein YnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,5750GT95.9% 166.0 / ‑0.6 49S41I (AGC→ATC) ynaARac prophage; putative prophage tail length tape measure domain‑containing protein YnaA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (24/23);  total (25/24)
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

GGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATAGCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTATTCTGAATTAAAGAAATTTTTCGATA  >  NC_000913/1426441‑1426709
                                                                                                                                      |                                                                                                                                      
ggAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAg                                                                                                                                   <  1:11395/140‑1 (MQ=255)
 gacttttAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAgg                                                                                                                                  <  2:595621/137‑1 (MQ=255)
    ttttAATAAATGGGAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGaa                                                                                                                               <  1:708986/140‑1 (MQ=255)
       taataaATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGt                                                                                                                            >  1:444672/1‑140 (MQ=255)
           aaaTGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACt                                                                                                                        <  2:444672/140‑1 (MQ=255)
                      cAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACt                                                                                                                        <  1:317951/129‑1 (MQ=255)
                      cAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACt                                                                                                                        >  2:317951/1‑129 (MQ=255)
                         aCCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACagaaga                                                                                                          <  1:729975/140‑1 (MQ=255)
                              gcgTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATgaga                                                                                                     >  2:480161/1‑140 (MQ=255)
                                            ttattGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAGTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCa                                                                                       >  2:236194/1‑140 (MQ=255)
                                            ttattGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCa                                                                                       >  2:332822/1‑140 (MQ=255)
                                              attGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCaaa                                                                                     >  2:601213/1‑140 (MQ=255)
                                                    aGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATAGCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTaaa                                                                               >  1:147477/1‑140 (MQ=255)
                                                      cattagAGCGACAAGTAATGCTGACAATCTTACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATAGCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTaaaaa                                                                             >  2:794082/6‑139 (MQ=255)
                                                      cACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAag                                                                             >  1:396847/1‑140 (MQ=255)
                                                      cACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAag                                                                             <  2:285792/140‑1 (MQ=255)
                                                           gAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAag                                                                             >  1:838421/1‑135 (MQ=255)
                                                           gAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAag                                                                             <  2:838421/135‑1 (MQ=255)
                                                            aGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGATTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTa                                                                       <  2:396847/140‑1 (MQ=255)
                                                              cGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTAtt                                                                     <  2:647729/140‑1 (MQ=255)
                                                              cGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTAtt                                                                     >  1:647729/1‑140 (MQ=255)
                                                                 aaaaTAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGAt                                                                  <  1:601213/140‑1 (MQ=255)
                                                                   aataatGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATcc                                                                >  1:494915/1‑140 (MQ=255)
                                                                           tGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTg                                                        <  2:320766/140‑1 (MQ=255)
                                                                             aCAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGcc                                                      <  1:358368/140‑1 (MQ=255)
                                                                                    gACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTg                                               >  2:957543/1‑140 (MQ=255)
                                                                                    gACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTg                                               >  2:583932/1‑140 (MQ=255)
                                                                                     aCCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGa                                              <  1:583932/140‑1 (MQ=255)
                                                                                            cATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtc                                       >  2:539352/1‑140 (MQ=255)
                                                                                               tAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  1:957543/140‑1 (MQ=255)
                                                                                               tAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    >  2:955976/1‑140 (MQ=255)
                                                                                               tAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  1:260876/140‑1 (MQ=255)
                                                                                               tAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATAGCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  2:147477/140‑1 (MQ=255)
                                                                                                 gTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  1:926382/138‑1 (MQ=255)
                                                                                                 gTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    >  2:926382/1‑138 (MQ=255)
                                                                                                  tCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggtcgc                                 >  2:87627/1‑139 (MQ=255)
                                                                                                    aaagTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACg                               <  1:236194/140‑1 (MQ=255)
                                                                                                    aaaaTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACg                               >  2:30139/1‑140 (MQ=255)
                                                                                                    aaaaTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACg                               <  1:332822/140‑1 (MQ=255)
                                                                                                      aaTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  2:396245/133‑1 (MQ=255)
                                                                                                      aaTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    >  1:396245/1‑133 (MQ=255)
                                                                                                       aTGACGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTa                            >  2:929961/1‑140 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  2:690089/128‑1 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    <  2:837263/128‑1 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    >  1:837263/1‑128 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTgtcggt                                    >  1:690089/1‑128 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTATTCt                        >  2:418128/1‑140 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTATTCt                        >  1:917374/1‑140 (MQ=255)
                                                                                                           cGCAGGCTGGGGAAACTGCTGCAAATATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTATTCt                        >  1:689432/1‑140 (MQ=255)
                                                                                                                                 aaaTATCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTATTCTGAATTAAAGAAATTTTTCGATa  <  1:955976/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCAAATAGCGCAGGGAAAGTGACTAAGGCAACAGAAGATGAGAAGAACGCGCTCGCAAAATTAAAAGCAGCTATTGATCCAGTTGGTGCCGCAATTGATACTGTCGGTCGACGCTATTCTGAATTAAAGAAATTTTTCGATA  >  NC_000913/1426441‑1426709

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