Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 12.3% intergenic (+34/‑164) lldD → / → trmL L‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0510CG12.3% 93.5 / 2.8 36intergenic (+34/‑164)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (14/17);  new base G (2/3);  total (16/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

TGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCT  >  NC_000913/3780912‑3781184
                                                                                                                                           |                                                                                                                                     
tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                    <  1:246665/143‑1 (MQ=255)
    gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                <  2:70260/143‑1 (MQ=255)
     ggttaTCGGTCAGCGCAATTACGTGAGTTTCGCTGGTGCATGGGCTGGGTAAAGAGTTGCCGGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCTCACCCTAACCCTCTCCCTGATGGagag                                                                                                                               <  2:156482/139‑1 (MQ=255)
      gAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAgg                                                                                                                              >  1:96075/1‑143 (MQ=255)
       aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCAGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAggg                                                                                                                             <  1:259206/143‑1 (MQ=255)
          tCGATCAGCGAAATGACGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                                          <  1:201393/143‑1 (MQ=255)
                      aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                              >  1:283150/1‑143 (MQ=255)
                      aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                              >  1:517550/1‑143 (MQ=255)
                       cgaCGCAAGAGTCGCTGGTGCCGGGGCCTGGTAAAGAGGTGCCTGCGGCGCTGGCGCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGAACCGTTCGGCGCtgt                                                                                                             <  2:436192/141‑1 (MQ=255)
                             aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAc                                                                                                       >  1:297608/1‑143 (MQ=255)
                                            aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGGGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                        <  2:206978/143‑1 (MQ=255)
                                                 cTGGGTAAAGAGCGGCCTGGGGCACAGGCTCCCTGTGCGCCCGCGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                   <  2:247173/143‑1 (MQ=37)
                                                    ggTAAGGAGTTGCCGGCGGCACTGGCCCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCg                                                                                <  2:411557/143‑1 (MQ=255)
                                                    ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACTCTGCAACGCCACCg                                                                                >  1:482819/1‑143 (MQ=255)
                                                         agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGCCCGTCCGGCGCTGTCTGTACTCCCTCACTCTGAAACGACACCGCAcgc                                                                           >  2:384076/1‑141 (MQ=255)
                                                                       cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                             <  1:271748/143‑1 (MQ=255)
                                                                            gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCGcccc                                                        >  1:304594/1‑143 (MQ=255)
                                                                              tCCCACGGCGAAAGGGAATGCGGCATCGTCGTTTGCGCCCCCGCCCCCTAACCCTCGCCCTCAGGGAGAGGGGACCGTCCGGCGCTGTATGTACTCCCTCTCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTc                                                      <  2:120799/143‑1 (MQ=255)
                                                                                ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCTGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                    <  2:100647/143‑1 (MQ=255)
                                                                                 cATGGCGAAAGGGCATGCGGCATAGTCGTCTGCGCCCCCTCGCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCGGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                   <  2:283150/143‑1 (MQ=255)
                                                                                    gtcGAAGGGGGGTGCGGCATAGGCGGTGGTGCCCCCGCAGCCGAGCCCTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTGCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                <  2:384032/141‑1 (MQ=255)
                                                                                     gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggtgg                                               >  1:259495/1‑140 (MQ=255)
                                                                                          aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGACACCACACCGCACTCGTTTTTTCTTCCTCGCCCCTGCGGGGAGAggac                                          >  2:68169/1‑141 (MQ=255)
                                                                                              aaTGCGGCATAGTCGTTTTCGCCCCCTCACCCTCACCCTCTCCCTCAGGGAGGGGGGACCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                      <  2:401408/143‑1 (MQ=255)
                                                                                                  cGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGCTCGGCGCTGTATATACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGCGGCGACAGCCggggcca                                  >  2:90447/1‑140 (MQ=39)
                                                                                                  cGGCATAGTCGTTTGCGCCCCCTATACCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                  <  1:260333/143‑1 (MQ=255)
                                                                                                    gCATAGTCGTTTGCGCCCCCTTACCCTAAGCCTGTTCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                <  1:58548/143‑1 (MQ=255)
                                                                                                    gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                <  1:428180/143‑1 (MQ=255)
                                                                                                     cATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACCCTGAAAAGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgcg                               >  2:503372/1‑141 (MQ=255)
                                                                                                      aTAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                              <  2:32609/143‑1 (MQ=255)
                                                                                                        aGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                            >  1:420394/1‑143 (MQ=255)
                                                                                                         gTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGGGGGGaaa                           >  1:399276/1‑143 (MQ=255)
                                                                                                           cgccggcgcCCCGTCAGCCCATCCATCCCCCTCAGGGCGAGGCGACCGTTCGGTGCTGTACGTACTCCCGAACTCTGAGCGGACACCGGACTCTTTTTTTTTACCTCGCCCGTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                         <  1:172810/138‑1 (MQ=255)
                                                                                                             tttGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                       >  1:181310/1‑143 (MQ=255)
                                                                                                                                  ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt  >  1:41377/1‑143 (MQ=255)
                                                                                                                                  ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGGGGGCCGGGGCGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt  >  1:197692/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                     
TGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCT  >  NC_000913/3780912‑3781184

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