Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,200,881 G→A A51V (GCC→GTC)  xisE ← e14 prophage, putative excisionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,8810GA97.8% 135.6 / ‑6.8 46A51V (GCC→GTC) xisEe14 prophage, putative excisionase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/22);  minor base C (1/0);  total (24/22)
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.85e-01

TACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCA  >  NC_000913/1200746‑1201005
                                                                                                                                       |                                                                                                                            
tactaTAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGt                                                                                                                           >  2:218976/1‑139 (MQ=255)
        aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTACCTCACCTTCCTCACga                                                                                                                   >  1:334827/1‑138 (MQ=255)
        aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                   >  2:464192/1‑139 (MQ=255)
        aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                   >  2:389101/1‑139 (MQ=255)
        aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGc                                                                                                                   >  1:608754/1‑139 (MQ=255)
        aGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCACCAACTAACTCGACGTCCTCACGc                                                                                                                   >  2:617664/1‑139 (MQ=255)
         ggggATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCCAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                  <  2:550314/137‑1 (MQ=255)
         gCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCa                                                                                                                  <  1:486931/139‑1 (MQ=255)
            gATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTc                                                                                                               >  1:296660/1‑139 (MQ=255)
            gATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTc                                                                                                               >  1:46052/1‑139 (MQ=255)
                            tGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTg                                                                                               >  2:678149/1‑139 (MQ=255)
                             gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGt                                                                                              >  1:76221/1‑139 (MQ=255)
                                   ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                    <  2:725229/111‑1 (MQ=255)
                                   ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                    >  1:725229/1‑111 (MQ=255)
                                   ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATg                                                                                        <  1:871644/139‑1 (MQ=255)
                                              cgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACa                                                                             >  2:141440/1‑139 (MQ=255)
                                                gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                           >  1:177298/1‑139 (MQ=255)
                                                   tGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaata                                                                        >  2:281618/1‑139 (MQ=255)
                                                       tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  2:296660/139‑1 (MQ=255)
                                                       tGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  1:617664/139‑1 (MQ=255)
                                                          cATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTg                                                                 <  2:909770/139‑1 (MQ=255)
                                                            ttACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                               <  1:141440/139‑1 (MQ=255)
                                                                        cGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                   <  2:177298/139‑1 (MQ=255)
                                                                        cGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                   <  1:281618/139‑1 (MQ=255)
                                                                                  gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGa                                         <  2:440079/139‑1 (MQ=255)
                                                                                   tAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGac                                        <  2:46052/139‑1 (MQ=255)
                                                                                            aTTTTTCTTGATTACAGGTGTTACGAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgga                               >  2:631518/1‑139 (MQ=255)
                                                                                             tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                    <  1:852923/53‑1 (MQ=255)
                                                                                             tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACg                                                                                                                    >  2:852923/1‑53 (MQ=255)
                                                                                                ttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAt                           <  1:218976/139‑1 (MQ=255)
                                                                                                          cAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATa                                                        <  2:445168/100‑1 (MQ=255)
                                                                                                          cAGGTGTTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATa                                                        >  1:445168/1‑100 (MQ=255)
                                                                                                               gtTACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttttt            >  2:674042/1‑139 (MQ=255)
                                                                                                                ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCga                                                                                                    <  1:80975/50‑2 (MQ=255)
                                                                                                                ttACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCga                                                                                                    >  2:80975/1‑49 (MQ=255)
                                                                                                                 tACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                               <  2:517940/86‑1 (MQ=255)
                                                                                                                 tACCAACTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                               >  1:517940/1‑86 (MQ=255)
                                                                                                                       cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc    >  1:305312/1‑139 (MQ=255)
                                                                                                                       cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc    >  1:446469/1‑139 (MQ=255)
                                                                                                                       cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc    <  2:334827/139‑1 (MQ=255)
                                                                                                                       cTCCCCAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc    <  1:759961/139‑1 (MQ=255)
                                                                                                                         ccccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa  <  2:363383/139‑1 (MQ=255)
                                                                                                                          cccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    >  2:782154/1‑72 (MQ=255)
                                                                                                                          cccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                    <  1:782154/72‑1 (MQ=255)
                                                                                                                          cccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACtt                      <  2:642916/118‑1 (MQ=255)
                                                                                                                          cccAACTAACTCGACGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACtt                      >  1:642916/1‑118 (MQ=255)
                                                                                                                                       |                                                                                                                            
TACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCA  >  NC_000913/1200746‑1201005

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