Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,330,494 T→C 9.2% intergenic (+19/‑288) PP_5SE → / → PP_1157 5S ribosomal RNA/acetolactate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,330,4940TC 9.2% 144.9 / 3.4 52intergenic (+19/‑288)PP_5SE/PP_11575S ribosomal RNA/acetolactate synthase
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base C (3/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAG  >  NC_002947/1330388‑1330642
                                                                                                          |                                                                                                                                                    
gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                           <  1:717784/150‑1 (MQ=31)
gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                           <  1:398768/150‑1 (MQ=31)
  tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTgg                                                                                                         <  1:1247565/150‑1 (MQ=32)
    ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTTTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                       >  2:27468/1‑150 (MQ=21)
          cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTc                                                                                                 >  2:1443786/1‑150 (MQ=38)
              cTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCgaga                                                                                             >  2:728609/1‑150 (MQ=39)
               tacgtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGagag                                                                                            >  2:1426767/4‑150 (MQ=255)
               tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGagag                                                                                            >  1:1625463/1‑150 (MQ=255)
                cAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGTGc                                                                                           <  1:876422/150‑1 (MQ=255)
                     gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                      <  2:817681/150‑1 (MQ=255)
                      tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGAt                                                                                     <  2:610651/150‑1 (MQ=255)
                      tGAAACGATGCATCGCAGATGGTAGTGTGGGGTTTCACCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAACTGATCGAGGGCGGGGAt                                                                                     >  2:148583/1‑150 (MQ=14)
                       gAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGGCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATg                                                                                    <  2:1613820/150‑1 (MQ=255)
                            gATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCAcc                                                                               <  1:1443786/150‑1 (MQ=255)
                                        aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGagag                                                                                            >  2:408853/1‑125 (MQ=255)
                                        aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGagag                                                                                            <  1:408853/125‑1 (MQ=255)
                                           gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGtt                                                                >  2:1030035/1‑150 (MQ=255)
                                               tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                            <  2:1043378/150‑1 (MQ=255)
                                                tgaggatttCCCCATGTGAGCGTCGGTCATCGTTTTCACTCATTTCTCCCTTCCCCTACCACCACACGCAGGTACTGTTTTTGTCTTTGTTTGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                           <  2:1251878/144‑1 (MQ=255)
                                                      tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAc                                                     <  2:1625463/150‑1 (MQ=255)
                                                           ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTc                                                               >  2:1222495/1‑135 (MQ=255)
                                                           ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTc                                                               <  1:1222495/135‑1 (MQ=255)
                                                             aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTcc                                              >  1:809287/1‑150 (MQ=255)
                                                               gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTa                                            <  1:37797/150‑1 (MQ=255)
                                                                      tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGa                                     <  1:1636592/150‑1 (MQ=255)
                                                                        gggcATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAGTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATg                                   <  2:162238/147‑1 (MQ=255)
                                                                          tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                          >  1:1136325/1‑93 (MQ=255)
                                                                          tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                          <  2:1136325/93‑1 (MQ=255)
                                                                          tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGt                                 >  1:754117/1‑150 (MQ=255)
                                                                               gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                            >  1:646138/1‑150 (MQ=255)
                                                                                tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                      >  2:390599/1‑91 (MQ=255)
                                                                                tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                      <  1:390599/91‑1 (MQ=255)
                                                                                  aaGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGcc                                                                      <  2:1240853/105‑1 (MQ=255)
                                                                                  aaGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGcc                                                                      >  1:1240853/1‑105 (MQ=255)
                                                                                  aaGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATg                                   <  1:847672/140‑1 (MQ=255)
                                                                                  aaGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATg                                   >  2:847672/1‑140 (MQ=255)
                                                                                    gATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAg                       >  1:1432004/1‑150 (MQ=255)
                                                                                      ttCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGc                     <  2:809287/150‑1 (MQ=255)
                                                                                          tttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                           <  1:437803/108‑1 (MQ=255)
                                                                                          tttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                           >  2:437803/1‑108 (MQ=255)
                                                                                          tttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGa                 >  2:659388/1‑150 (MQ=255)
                                                                                           ttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGaa                <  1:1030035/150‑1 (MQ=255)
                                                                                            tCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAAt               >  2:1034051/1‑150 (MQ=255)
                                                                                             cGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATg              <  1:156990/150‑1 (MQ=255)
                                                                                               cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc            <  1:27468/150‑1 (MQ=255)
                                                                                                aaaaCCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGcc           <  1:1034051/150‑1 (MQ=255)
                                                                                                    ccccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTGGCGGCAAAACCCTCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCaa       >  2:862799/1‑150 (MQ=255)
                                                                                                     cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                        >  2:1457478/1‑68 (MQ=18)
                                                                                                     cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                        <  1:1457478/68‑1 (MQ=18)
                                                                                                     cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc      >  1:913785/1‑150 (MQ=38)
                                                                                                      ccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACag                                         >  2:1055233/1‑114 (MQ=37)
                                                                                                         aTCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgcag  <  2:754117/150‑1 (MQ=255)
                                                                                                          |                                                                                                                                                    
GATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAG  >  NC_002947/1330388‑1330642

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