Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,199,764 (T)5→6 coding (1043/1128 nt) intE ← e14 prophage, putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,199,7591.T97.4% 115.0 / NA 38L350I (CTA→ATA) intEe14 prophage, putative integrase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (18/19);  minor base C (0/1);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTCATCGTCAA  >  NC_000913/1199622‑1199891
                                                                                                                                          |                                                                                                                                    
tttAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAagtagt                                                                                                                                     <  1:118771/139‑1 (MQ=255)
       tttATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttgt                                                                                                                              <  1:825781/139‑4 (MQ=255)
                  cGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAAccc                                                                                                                    >  2:1406819/1‑139 (MQ=255)
                    gTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCt                                                                                                                  <  2:1345168/139‑1 (MQ=255)
                    gTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCt                                                                                                                  <  1:1406819/139‑1 (MQ=255)
                               aaTGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttg                                                                                                                               <  1:672191/114‑3 (MQ=255)
                               aaTGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttg                                                                                                                               >  2:672191/1‑112 (MQ=255)
                                aTGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACa                                                                                                      <  1:768969/139‑1 (MQ=255)
                                                cTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttgtg                                                                                                                             >  1:1504344/1‑95 (MQ=255)
                                                cTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttgtg                                                                                                                             <  2:1504344/99‑5 (MQ=255)
                                                       gTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGtttt                                                                                                                                  <  2:484285/87‑1 (MQ=255)
                                                       gTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGtttt                                                                                                                                  >  1:484285/1‑87 (MQ=255)
                                                          cacaACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCa                                                                            >  2:1056615/1‑139 (MQ=255)
                                                          cacaACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCa                                                                            >  1:1120830/1‑139 (MQ=255)
                                                          cacaACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttgt                                                                                                                              >  1:987257/1‑85 (MQ=255)
                                                          cacaACGCTACTTTGCTCCATCCTTAACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGCTTTTtgt                                                                                                                               <  2:987257/88‑1 (MQ=255)
                                                                       tGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAg                                                               >  2:290626/1‑139 (MQ=255)
                                                                           ccATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGa                                                           <  2:1120830/139‑1 (MQ=255)
                                                                              tCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAAccc                                                                                                                    <  1:913182/79‑1 (MQ=255)
                                                                              tCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAAccc                                                                                                                    >  2:913182/1‑79 (MQ=255)
                                                                                         gATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGt                                             <  1:290626/139‑1 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTaa                                                                                                                       >  1:709584/1‑64 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTaa                                                                                                                       <  2:709584/64‑1 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGtt                                                                                                >  1:243937/1‑87 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGtt                                                                                                <  2:243937/87‑1 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTc                                            >  1:1334980/1‑139 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttgt                                                                                                                              <  2:1277802/56‑4 (MQ=255)
                                                                                          atcatcATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGttttttgt                                                                                                                              >  1:1277802/1‑53 (MQ=255)
                                                                                              tcatGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAAtt                                        <  2:1334980/139‑1 (MQ=255)
                                                                                               catGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTc                                       >  2:1041553/1‑139 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCt                                                                                                                         <  2:77760/51‑1 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCt                                                                                                                         >  1:77760/1‑51 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTta                                                                                                          >  2:735989/1‑66 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTta                                                                                                          <  1:735989/66‑1 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCg                                                                  <  1:325773/106‑1 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCg                                                                  >  2:325773/1‑106 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTg                                                            >  2:1225278/1‑112 (MQ=255)
                                                                                                     aaCGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTg                                                            <  1:1225278/112‑1 (MQ=255)
                                                                                                       cGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTaaa                               >  1:1196337/1‑139 (MQ=255)
                                                                                                         aTCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCt                                                                                                                  <  2:1483602/54‑1 (MQ=255)
                                                                                                         aTCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCt                                                                                                                  >  1:1483602/1‑54 (MQ=255)
                                                                                                         aTCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATa                                                                                          >  1:662329/1‑78 (MQ=255)
                                                                                                         aTCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATa                                                                                          <  2:662329/78‑1 (MQ=255)
                                                                                                                 gttgttgtGTTCTATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATTCAGACGTTCAGATAAAGACCTTTGCTCATGGAATGGCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCtt                     <  1:979163/139‑1 (MQ=255)
                                                                                                                            ttATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTc          >  1:156143/1‑139 (MQ=255)
                                                                                                                             tATGACCAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTCa         <  2:1196337/139‑1 (MQ=255)
                                                                                                                                   cAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTCATCGTCa   <  1:1056615/139‑1 (MQ=255)
                                                                                                                                   cAAGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTCATCGTCa   <  1:1041553/139‑1 (MQ=255)
                                                                                                                                    aaGTAGTTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTCATCGTCaa  >  1:416416/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                    
TTTAGTTTTTATCATGTGCGGTGAGAAATTCAATGTGGCGTTGAGATGCTTAAAGGTTCACAACGCTACTTTGCTCCATCCTTTACCTCGATCATCATGATAACGATCGGTTTGTTGTTGTGTTTTATGACCAAGTAGTTTTTGTGTGTCTAACCCCTGTTCTTTATACAGACGTTCAGATAAAGACCTTTGCTCATGGAATGTCGCAGGTGAACCCTCTCCCCAGTCAATTCTTGCTAAATCTCTCGCTTTACTAAAATTCATCGTCAA  >  NC_000913/1199622‑1199891

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: