Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,781,031 (C)5→4 intergenic (+14/‑184) 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,0270C.97.5% 164.1 / ‑1.8 40intergenic (+10/‑188)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (16/23);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGC  >  NC_000913/3780890‑3781144
                                                                                                                                         |                                                                                                                     
tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                     <  1:166468/140‑1 (MQ=255)
  aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccct                                                                                                                   <  1:1096358/140‑2 (MQ=255)
  aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccct                                                                                                                   <  1:470942/140‑2 (MQ=255)
  aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccct                                                                                                                   <  1:397302/140‑2 (MQ=255)
       ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCAccc                                                                                                             >  1:358324/1‑140 (MQ=255)
       ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCAccc                                                                                                             <  1:234750/140‑1 (MQ=255)
        gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCt                                                                                                            <  2:282031/140‑1 (MQ=255)
          gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTaa                                                                                                          <  2:99659/140‑1 (MQ=255)
                gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCctct                                                                                                    <  2:39519/140‑1 (MQ=255)
                gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCctct                                                                                                    <  2:358324/140‑1 (MQ=255)
                  tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTcc                                                                                                  >  1:124523/1‑140 (MQ=255)
                     cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTc                                                                                               <  2:996888/140‑1 (MQ=255)
                        gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                            >  1:526851/1‑140 (MQ=255)
                              aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                      >  2:190121/1‑140 (MQ=255)
                                      aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                              <  2:304263/140‑1 (MQ=255)
                                        cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCg                                                                            <  2:209548/140‑1 (MQ=255)
                                          aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                          <  1:38241/140‑1 (MQ=255)
                                           aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                         <  1:355122/140‑1 (MQ=255)
                                           aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                         <  1:231843/140‑1 (MQ=255)
                                           aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                         <  1:194295/140‑1 (MQ=255)
                                             ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                       <  2:618354/140‑1 (MQ=255)
                                                cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                    >  2:211870/1‑140 (MQ=255)
                                                    aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                <  2:1012484/140‑1 (MQ=255)
                                                          gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                      <  1:1037510/128‑1 (MQ=255)
                                                          gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                      >  2:1037510/1‑128 (MQ=255)
                                                          gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                          <  1:1061991/140‑1 (MQ=255)
                                                                    gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaa                                                 >  2:507000/1‑140 (MQ=255)
                                                                      gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                              >  1:614256/1‑140 (MQ=255)
                                                                        tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGac                                            >  1:922001/1‑140 (MQ=255)
                                                                        tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGac                                            >  2:406114/1‑140 (MQ=255)
                                                                         gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                           >  1:31052/1‑140 (MQ=255)
                                                                              aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                      >  2:843699/1‑140 (MQ=255)
                                                                                     gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                               >  2:801993/1‑140 (MQ=255)
                                                                                           ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                         >  1:187731/1‑140 (MQ=255)
                                                                                                tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                    <  2:922001/140‑1 (MQ=255)
                                                                                                tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                    >  2:652632/1‑140 (MQ=255)
                                                                                                tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                    <  2:702190/140‑1 (MQ=255)
                                                                                                  gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                  >  2:292208/1‑140 (MQ=255)
                                                                                                  gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                  <  1:843699/140‑1 (MQ=255)
                                                                                                   cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                 <  1:190121/140‑1 (MQ=255)
                                                                                                   cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                 <  1:652632/140‑1 (MQ=255)
                                                                                                           gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg         >  2:706975/1‑140 (MQ=255)
                                                                                                            cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga        <  1:292208/140‑1 (MQ=255)
                                                                                                                  ggAATGCGGCATAGTCGTTTGCGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc  <  2:526851/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                     
TGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGC  >  NC_000913/3780890‑3781144

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