Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 568,239 A→G pseudogene (205/213 nt) renD → DLP12 prophage; protein RenD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,2390AG98.1% 170.0 / ‑4.4 53pseudogene (205/213 nt)renDDLP12 prophage; protein RenD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/31);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTG  >  NC_000913/568109‑568363
                                                                                                                                  |                                                                                                                            
ggTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAgca                                                                                                                            <  2:473174/133‑1 (MQ=255)
      ttACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAatga                                                                                                                      <  1:316200/133‑1 (MQ=255)
         ccGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAAtgatga                                                                                                                   <  1:369517/133‑1 (MQ=255)
           ggTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGaaa                                                                                                                 <  1:310385/133‑1 (MQ=255)
           ggTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGaaa                                                                                                                 <  1:401134/133‑1 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGGGGGTATATGAAAGAACGTCAATGGGGGCACAATGATGAAAc                                                                                                                >  2:135192/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGGGGGTATATCAAATAACTTCAATGGGTGCACAATGATGAAAc                                                                                                                >  2:175992/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATGAAAAAACAACAATGGGAGCACCTTGATGAAAc                                                                                                                >  2:251409/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                                                                >  2:72593/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                                                                >  2:403231/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                                                                >  1:494600/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                                                                >  1:256854/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                                                                >  1:29398/1‑133 (MQ=255)
            gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                                                                >  2:315446/1‑133 (MQ=255)
              tCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGTGTATATAAGAGAGCATGAATGGGAGGACAAGGATGAAACAg                                                                                                              >  2:132127/1‑133 (MQ=255)
              tCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAg                                                                                                              >  2:349290/1‑133 (MQ=255)
              tCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAg                                                                                                              >  1:53417/1‑133 (MQ=255)
                              aTGGGAAGGAAAGGTGAGCATGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCaaa                                                                                              <  1:409492/133‑1 (MQ=255)
                                gggAAGGAAACGTGATCACGAATCGGATTTTTGAGGAGTGTAGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACt                                                                                            <  2:195634/133‑1 (MQ=255)
                                 ggaaggaaAGGTGAGTACGAATCTGATTTTCAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTg                                                                                           <  1:169195/133‑1 (MQ=255)
                                 ggaaggaaAGGTGAGTACGAATCTCATTATTATGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTg                                                                                           <  2:174014/133‑1 (MQ=255)
                                 ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTg                                                                                           <  1:286165/133‑1 (MQ=255)
                                  gaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTgt                                                                                          >  2:472814/1‑133 (MQ=255)
                                      gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtagta                                                                                      <  1:317995/133‑1 (MQ=255)
                                      gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtagta                                                                                      <  2:281754/133‑1 (MQ=255)
                                               gCACGAATTTGATTTTTAAGGACTGTCGCCGGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTGCAATtctct                                                                             <  1:251409/133‑1 (MQ=255)
                                               gCACGAATCTGATTTCTACGGACTGTCGCCTGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtctct                                                                             <  1:135192/133‑1 (MQ=255)
                                                cACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                                                                                <  2:387115/129‑1 (MQ=255)
                                                cACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                                                                                >  1:387115/1‑129 (MQ=255)
                                                cACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTACAAACTGTAGTACAATtctctc                                                                            >  1:204937/1‑133 (MQ=255)
                                                         gATTTTTAAGCAGTGCCGCCAGAGTGCCTTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGaa                                                                   <  2:138035/133‑1 (MQ=255)
                                                          aTTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAAc                                                                  <  2:315410/133‑1 (MQ=255)
                                                             tttAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAgg                                                               <  1:72593/133‑1 (MQ=255)
                                                             tttAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAgg                                                               <  2:29398/133‑1 (MQ=255)
                                                              ttAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGa                                                              <  1:315446/133‑1 (MQ=255)
                                                                    aGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAat                                                        >  1:449733/1‑133 (MQ=255)
                                                                    aGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAat                                                        >  1:505759/1‑133 (MQ=255)
                                                                      tgtCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAatat                                                      >  2:405634/1‑133 (MQ=255)
                                                                            cAGAGTGCCGTGATGCAACGGGCATTGAGGCTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                <  1:132127/133‑1 (MQ=255)
                                                                            cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                <  1:87884/133‑1 (MQ=255)
                                                                            cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                <  1:472814/133‑1 (MQ=255)
                                                                            cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                <  2:275363/133‑1 (MQ=255)
                                                                            cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                <  2:482347/133‑1 (MQ=255)
                                                                            cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                <  2:98320/133‑1 (MQ=255)
                                                                            cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTACAGTTTGAACAGGAAAGAATATGCTATg                                                <  2:445163/133‑1 (MQ=255)
                                                                              gagTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGGAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGat                                              <  2:53417/133‑2 (MQ=255)
                                                                                        aTGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATAtt                                    <  1:349290/133‑1 (MQ=255)
                                                                                                  tATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtggtg                          >  2:316925/1‑133 (MQ=255)
                                                                                                      gAGGGTATATAAAAGAACCTCATTGGGAGCCCAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGc                      <  1:175992/133‑1 (MQ=255)
                                                                                                                aaaaGAACATTAATGCGAGCACAATGATAAAAAAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCa            <  2:155355/133‑1 (MQ=255)
                                                                                                                   aGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCagag         <  2:505759/133‑1 (MQ=255)
                                                                                                                   aGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCGTGGAGGTGCAATACTTGCagag         >  1:206380/1‑133 (MQ=255)
                                                                                                                       cATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACTCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCa     >  2:386658/1‑133 (MQ=255)
                                                                                                                          cAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTg  >  2:7211/1‑133 (MQ=255)
                                                                                                                                  |                                                                                                                            
GGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTG  >  NC_000913/568109‑568363

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