Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,774,505 C→G 5.7% intergenic (+128/‑70) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,774,5050CG 5.7% 2289.5 / 18.6 658intergenic (+128/‑70)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (529/79);  new base G (2/47);  total (532/126)
Fisher's exact test for biased strand distribution p-value = 4.81e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 682 total aligned reads displayed.

ATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGC  >  CP009273/3774438‑3774571
                                                                   |                                                                  
atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                    >  1:33965641/1‑68 (MQ=255)
atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                    >  1:25537179/1‑68 (MQ=255)
atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                    >  1:8109165/1‑68 (MQ=255)
  gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                  <  1:13446101/68‑1 (MQ=255)
   taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                 >  1:23099923/1‑68 (MQ=255)
   taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                 >  1:36699000/1‑68 (MQ=255)
   taCTCCCTCACTCTGAAACGACACCGCACTCTCTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                 >  1:1002291/1‑68 (MQ=255)
    aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                <  1:22113809/68‑1 (MQ=255)
    aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                 <  1:35332082/67‑1 (MQ=255)
    aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                 <  1:2421921/67‑1 (MQ=255)
     cTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggga                                                               >  1:4550965/1‑67 (MQ=255)
     cTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:7420120/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggggg                                                               >  1:5287226/1‑67 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggggg                                                               >  1:33357185/1‑67 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:3479233/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:1945653/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:21462745/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:24411492/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:3277909/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:23571276/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:35896833/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:25098523/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:30078659/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:26429274/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:26894446/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:27553709/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:28355218/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:38245723/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:7702275/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:6383349/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:1293258/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:13963928/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:14314680/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:38497806/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:19176530/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:1547186/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:16482356/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:9334720/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:10588769/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:1708677/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:37003589/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:36792806/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:1856780/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                               >  1:9805448/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                              >  1:16159078/1‑68 (MQ=255)
     cTCCCTCACTCTGAAACGACACCGCACTC‑‑TTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:11580383/1‑68 (MQ=37)
      tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                              >  1:29876258/1‑68 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:32544907/1‑68 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:29029506/1‑68 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             <  1:33661303/68‑1 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:3808077/1‑68 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:2804489/1‑68 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:20908519/1‑68 (MQ=255)
       cccTCACTCTGAAACGACACCGCACTCTTTTTTACTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                             >  1:36291082/1‑68 (MQ=255)
          tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                          >  1:27210293/1‑68 (MQ=255)
           cACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                         <  1:23778531/68‑1 (MQ=255)
            aCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                                        >  1:3767779/1‑68 (MQ=255)
              tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                      >  1:3135200/1‑68 (MQ=255)
              tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                      >  1:32190308/1‑68 (MQ=255)
              tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                      >  1:11982404/1‑68 (MQ=255)
                   aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                 >  1:1241657/1‑68 (MQ=255)
                    aCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                 >  1:31275187/1‑67 (MQ=255)
                    aCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                 >  1:11030789/1‑67 (MQ=255)
                    aCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                 >  1:903915/1‑67 (MQ=255)
                      gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                               >  1:20201938/1‑67 (MQ=255)
                        cacCCCACTGTTTTTTTCCCCCCCCCCCCGCCGGGTAGTGGGCCGGGGTGTGGGGAAAAGGCCGCACt                                            <  1:25326850/68‑1 (MQ=255)
                                  tttttttCCCCCCCCCCCCGCGGGGGCGGGGGCGGGGGGGCGGGCACAAGGCCGCACTGCTCACCAAt                                  <  1:34414613/68‑1 (MQ=255)
                                      tttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                >  1:10249443/1‑66 (MQ=255)
                                      tttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                >  1:4910923/1‑66 (MQ=255)
                                       tgctcCGTCCCCCCTCGGGGTAGTGGGCGGGGGTTGGGGTAAAAGGCCGCACTGCTCACCAATTTTCt                             <  1:16706760/66‑1 (MQ=255)
                                        tctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                             >  1:5743922/1‑67 (MQ=255)
                                        tctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                             >  1:6711668/1‑67 (MQ=255)
                                        tctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                             >  1:17402026/1‑67 (MQ=255)
                                            ccTCGCCCCTCCGGGGAGAGGGGGGGGGTGCGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTAt                        <  1:6926284/68‑1 (MQ=37)
                                                 ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                    >  1:1759994/1‑67 (MQ=255)
                                                  cccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                   >  1:2868504/1‑67 (MQ=255)
                                                  cccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                   >  1:10808967/1‑67 (MQ=255)
                                                    cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                 >  1:1522081/1‑67 (MQ=255)
                                                     tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc               >  1:22363919/1‑68 (MQ=255)
                                                     cggggggcgggggcgGGGGGGTGGGCACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc               <  1:23320753/65‑1 (MQ=255)
                                                       gggggTGTGGGCGGGGGTGCGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg              <  1:1130063/66‑1 (MQ=255)
                                                       cGGGGTGTGGGCGGGGGTGTGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc             <  1:38918087/68‑1 (MQ=255)
                                                       cGGGGTGCGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc             <  1:35508992/68‑1 (MQ=255)
                                                         gggAGTGGGCGGGGGTGTGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa            <  1:30792360/67‑1 (MQ=255)
                                                           gagaGGGCCGGGGTGAGGGGAAAAGGCCGCGCTGCTCACCAATTTTCTGCAATCCTGCCCCCGCACTa         >  1:8639708/1‑68 (MQ=39)
                                                            tgcgggCCGGGGTGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa         <  1:18253408/64‑1 (MQ=255)
                                                            cggggggCGGGGTGTGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa         <  1:14868002/64‑1 (MQ=255)
                                                            agtgggCGGGGGTTCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa         <  1:35002125/64‑1 (MQ=255)
                                                            agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        >  1:13731191/1‑68 (MQ=255)
                                                            agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        >  1:6039407/1‑68 (MQ=255)
                                                            agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        >  1:20592059/1‑68 (MQ=255)
                                                            agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        >  1:4031045/1‑68 (MQ=255)
                                                            agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        >  1:30392726/1‑68 (MQ=255)
                                                            agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        >  1:15704538/1‑68 (MQ=255)
                                                             gtgggCCGGGGGGCGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg       <  1:8832079/66‑1 (MQ=255)
                                                             gcgggcgggggggcggggACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa        <  1:22850121/65‑1 (MQ=255)
                                                             gcgggCCGGGGTGAGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg       <  1:31804898/66‑1 (MQ=255)
                                                             gagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg       >  1:25868289/1‑68 (MQ=255)
                                                             gagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg       >  1:18930212/1‑68 (MQ=255)
                                                                 gCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                               >  1:29450715/1‑40 (MQ=255)
                                                                   |                                                                  
ATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGC  >  CP009273/3774438‑3774571

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: