Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 892,869 G→A *301* (TAG→TAA rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913892,8690GA95.8% 148.5 / ‑0.7 48*301* (TAG→TAArimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/24);  minor base T (1/1);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

GGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCT  >  NC_000913/892737‑893003
                                                                                                                                    |                                                                                                                                      
ggTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                                  <  1:85026/139‑1 (MQ=255)
     tGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATcaca                                                                                                                             <  1:627685/139‑1 (MQ=255)
          gTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATGCGCTGGATCGACCGCCACGCTACGACAGAATATGGCCTGAAAACGGGGGGTTGATCGCAATCACATTAcc                                                                                                                        <  1:1011908/139‑2 (MQ=255)
           tGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                                          <  2:826117/136‑1 (MQ=255)
           tGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                                          >  1:826117/1‑136 (MQ=255)
                           gCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTg                                                                                                       <  1:772265/139‑1 (MQ=255)
                           gCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTg                                                                                                       <  1:294074/139‑1 (MQ=255)
                                aaggaaTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                                  <  1:1023784/139‑1 (MQ=255)
                                     aTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTATTCGCAATCACATTa                                                                                                                          <  1:166058/110‑1 (MQ=255)
                                     aTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTATTCGCAATCACATTa                                                                                                                          >  2:166058/1‑110 (MQ=255)
                                              accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTg                                                                                                       >  2:447708/1‑120 (MQ=255)
                                              accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTg                                                                                                       <  1:447708/120‑1 (MQ=255)
                                              accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCtgt                                                                                    >  2:1317272/1‑139 (MQ=255)
                                                   cGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                    <  1:261631/86‑1 (MQ=255)
                                                   cGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                    >  2:261631/1‑86 (MQ=255)
                                                               cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                    <  2:1010013/74‑1 (MQ=255)
                                                               cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                    >  1:1010013/1‑74 (MQ=255)
                                                               cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                                  <  2:907006/76‑1 (MQ=255)
                                                               cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                                  >  1:907006/1‑76 (MQ=255)
                                                                gcgGGTAAACTGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                                  <  1:1105953/139‑1 (MQ=255)
                                                                gcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                                  >  1:556837/1‑139 (MQ=255)
                                                                           gATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTAc                                                       <  2:556837/139‑1 (MQ=255)
                                                                            aTCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                    <  1:717655/61‑1 (MQ=255)
                                                                            aTCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                    >  2:717655/1‑61 (MQ=255)
                                                                              ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGcc                                                                                 >  2:641037/1‑110 (MQ=255)
                                                                              ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGcc                                                                                 <  1:641037/110‑1 (MQ=255)
                                                                              ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCTGGTCTTTTTATCGaa                                                                  <  2:675470/125‑1 (MQ=255)
                                                                              ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                                  >  1:675470/1‑125 (MQ=255)
                                                                                 cTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACaa                                                     <  1:1267024/135‑1 (MQ=255)
                                                                                 cTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACaa                                                     >  2:1267024/1‑135 (MQ=255)
                                                                                   ggATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAg                                                                                        <  1:135610/98‑1 (MQ=255)
                                                                                   ggATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAg                                                                                        >  2:135610/1‑98 (MQ=255)
                                                                                    gATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtttt                                                                                                <  1:574156/89‑1 (MQ=255)
                                                                                    gATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtttt                                                                                                >  2:574156/1‑89 (MQ=255)
                                                                                                                tGCCTGAAAACGGGTGGTTCATCACAATCACATGACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                                        <  1:797639/101‑1 (MQ=255)
                                                                                                                tGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATa                  >  1:1280481/1‑139 (MQ=255)
                                                                                                                tGCCTGAAAACGGGTGGTTAATCACAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                                        >  2:797639/1‑101 (MQ=255)
                                                                                                                   cTGAAAACGCGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGc               >  1:1259661/1‑139 (MQ=255)
                                                                                                                     gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCg                                                                                            >  2:1113845/1‑60 (MQ=255)
                                                                                                                     gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCg                                                                                            <  1:1113845/60‑1 (MQ=255)
                                                                                                                     gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                                        <  1:57145/96‑1 (MQ=255)
                                                                                                                     gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                                        >  2:57145/1‑96 (MQ=255)
                                                                                                                     gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTg             >  2:568507/1‑139 (MQ=255)
                                                                                                                     gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTg             >  1:718523/1‑139 (MQ=255)
                                                                                                                            ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTgc                                                                                                      <  1:395920/43‑1 (MQ=255)
                                                                                                                            ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTgc                                                                                                      >  2:395920/1‑43 (MQ=255)
                                                                                                                             gtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATg     <  1:1317272/139‑1 (MQ=255)
                                                                                                                                gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCt  >  1:465275/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCT  >  NC_000913/892737‑893003

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