Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,370 Δ1 bp 8.7% intergenic (‑100/+23) pheM ← / ← rplT phenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,3700A. 8.7% 158.2 / 9.6 46intergenic (‑100/+23)pheM/rplTphenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base . (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGAT  >  NC_000913/1799225‑1799512
                                                                                                                                                 |                                                                                                                                               
ttCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACt                                                                                                                                               <  2:707987/148‑1 (MQ=255)
                                             ttgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                 >  3:288935/1‑148 (MQ=255)
                                                          cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                 <  3:183113/135‑1 (MQ=255)
                                                          cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                 >  4:183113/1‑135 (MQ=255)
                                                           ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCgg                                                                                    >  4:571923/1‑147 (MQ=255)
                                                           ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCgg                                                                                    <  3:571923/147‑1 (MQ=255)
                                                            tttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTg                                                                                  >  1:389577/1‑148 (MQ=255)
                                                                  aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTAt                                                                                                                      >  2:795415/1‑106 (MQ=255)
                                                                  aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTAt                                                                                                                      <  1:795415/106‑1 (MQ=255)
                                                                        cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAgct                                                                                                             >  4:838037/1‑109 (MQ=255)
                                                                        cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAgct                                                                                                             <  3:838037/109‑1 (MQ=255)
                                                                          aaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATg                                                                                                                     <  3:177481/99‑1 (MQ=255)
                                                                          aaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATg                                                                                                                     >  4:177481/1‑99 (MQ=255)
                                                                                      tGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTa                                                                                                                       >  4:453848/1‑85 (MQ=255)
                                                                                      tGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTa                                                                                                                       <  3:453848/85‑1 (MQ=255)
                                                                                         aaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTTAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGat                                                        >  1:1031338/1‑145 (MQ=255)
                                                                                         aaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTTAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGat                                                        <  2:1031338/145‑1 (MQ=255)
                                                                                         aaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATc                                                     >  1:609431/1‑148 (MQ=255)
                                                                                                     gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCtt                                         <  3:267350/148‑1 (MQ=255)
                                                                                                       tAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTgg                                                                                                                           <  4:376793/64‑1 (MQ=255)
                                                                                                       tAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTgg                                                                                                                           >  3:376793/1‑64 (MQ=255)
                                                                                                        aaaaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCgg                                                                                    >  2:1129695/1‑102 (MQ=255)
                                                                                                        aaaaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCgg                                                                                    <  1:1129695/102‑1 (MQ=255)
                                                                                                             tgaaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTACAACTGGCTTATGCCAGAGTTGCTTTCGCtt                                                                                                   >  4:960071/3‑82 (MQ=255)
                                                                                                             tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                   <  3:960071/82‑1 (MQ=255)
                                                                                                             tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                >  4:944548/1‑101 (MQ=255)
                                                                                                             tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                <  3:944548/101‑1 (MQ=255)
                                                                                                             tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTAAGGTCGa                                 <  2:609431/148‑1 (MQ=255)
                                                                                                              cAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCaa                                                                                                                               >  1:491660/1‑53 (MQ=255)
                                                                                                              cAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCaa                                                                                                                               <  2:491660/53‑1 (MQ=255)
                                                                                                                aTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                   >  3:417907/1‑79 (MQ=255)
                                                                                                                aTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                   <  4:417907/79‑1 (MQ=255)
                                                                                                                    tttGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                  <  3:151482/76‑1 (MQ=255)
                                                                                                                    tttGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                  >  4:151482/1‑76 (MQ=255)
                                                                                                                               tGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                     >  1:788935/1‑94 (MQ=255)
                                                                                                                               tGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                     <  2:788935/94‑1 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCa            <  4:288935/148‑1 (MQ=255)
                                                                                                                                    cGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                     <  4:160678/89‑1 (MQ=255)
                                                                                                                                    cGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                     >  3:160678/1‑89 (MQ=255)
                                                                                                                                         agagGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                     <  3:959755/84‑1 (MQ=255)
                                                                                                                                         agagGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                     >  4:959755/1‑84 (MQ=255)
                                                                                                                                          gagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAgg                                              >  2:248946/1‑106 (MQ=255)
                                                                                                                                          gagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAgg                                              <  1:248946/106‑1 (MQ=255)
                                                                                                                                           agGGAGCTAGCCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAgg                                              >  2:268976/1‑105 (MQ=255)
                                                                                                                                           agGGAGCTAGCCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAgg                                              <  1:268976/105‑1 (MQ=255)
                                                                                                                                            gggAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATcgat  >  1:1124040/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGAT  >  NC_000913/1799225‑1799512

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