Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,672,170 T→C 12.1% T82A (ACG→GCG)  U712_13760 ← Sporulation initiation phosphotransferase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,672,1700TC12.1% 144.5 / 13.0 58T82A (ACG→GCG) U712_13760Sporulation initiation phosphotransferase B
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base C (3/4);  total (29/29)
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.43e-01

GACTGCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCTTTAATCAGCTGCAG  >  CP006881/2672039‑2672308
                                                                                                                                   |                                                                                                                                          
gACTGCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAAATTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGaaa                                                                                                                                      >  6:248870/1‑138 (MQ=255)
    gCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGCCAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAAACa                                                                                                                                  >  5:80560/1‑138 (MQ=255)
    gCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCGTGAATTTCTCCGAGAACTTCATATTCTAGCGTCATATTATGGGTTTTCCAATTAAACGTAAGAAAATCa                                                                                                                                  >  5:16730/1‑138 (MQ=255)
    gCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCAAAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCa                                                                                                                                  >  2:48799/1‑138 (MQ=255)
    gCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTCATCATACACCGTCACATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGTTTTTGCAATTTAACGCAAGAAAATCa                                                                                                                                  >  8:340759/1‑138 (MQ=255)
    gCTTGATCAAACAGATGAAACAGCTTTCTAATCACTTTCGCCAGCATTTGATCATAAGACGACAAATCCTTAATTTCTCAGAGAACTTCATATTCAAGAGTAATAAAATGGGTTTTCCAATTAAACGTAAGAAAAtaa                                                                                                                                  >  6:164237/1‑136 (MQ=255)
     cTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCa                                                                                                                                  <  6:348322/137‑1 (MQ=255)
     cTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCa                                                                                                                                  <  2:114146/137‑1 (MQ=255)
     cTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCa                                                                                                                                  <  4:256345/137‑1 (MQ=255)
        gATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAATCGCAAGAAAATCAAACGc                                                                                                                             >  8:121531/1‑139 (MQ=255)
          tCAAACAGATGAAACAGCTTTCTCATCAGTTTCTCCATCTTTTGATCATAAGCCGACAAATCGTTAATTTCTCGGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTATGAAAATCAATCGcc                                                                                                                            <  5:207971/138‑1 (MQ=255)
                  aTGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGgt                                                                                                                   >  8:89197/1‑139 (MQ=255)
                      aaCAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGttt                                                                                                               <  6:91358/139‑1 (MQ=255)
                        ctgCTTTCTCCTCTGTTTCGCCAGCGTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCTTCATATAATGGTTTTTCCAATTAAACGTAAGAAAATCAATCGCCAAATGCGGTGTTTTc                                                                                                             <  2:302288/137‑1 (MQ=255)
                        cAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTc                                                                                                             <  2:140537/139‑1 (MQ=255)
                               ctcATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGttt                                                                                                       <  7:124194/138‑1 (MQ=255)
                                  atcaGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTga                                                                                                     <  8:206251/137‑1 (MQ=255)
                                     agattCGCCAGCTTTTGATCATAAGCGGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACTTAAGAAAATCAAACGCCAACTGCGGTGGTTTCAGGTTTGAGAGCt                                                                                                <  5:129176/136‑1 (MQ=255)
                                     aGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCt                                                                                                <  5:248870/139‑1 (MQ=255)
                                       tttCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTACGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCttt                                                                                              >  2:183169/1‑139 (MQ=255)
                                        ttCGCCAGCTTTTTATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCttt                                                                                              <  3:330621/138‑1 (MQ=255)
                                            ccAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTAAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTc                                                                                         >  6:234755/1‑139 (MQ=255)
                                             cAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCg                                                                                        <  8:76148/139‑1 (MQ=255)
                                                    tGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCttt                                                                                  <  7:62411/138‑1 (MQ=255)
                                                    tGATCATAAGCCCACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCttt                                                                                  <  1:311742/138‑1 (MQ=255)
                                                     gATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTg                                                                                 <  2:336904/138‑1 (MQ=255)
                                                     gATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGCAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTg                                                                                 <  5:146311/138‑1 (MQ=255)
                                                     gATCATAAGCCGAAAAATCCTTAATTTCTCCGAGAACTTCATATTCAATCGTCATAAATTGGGTTTACCAATTAAACGCAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTg                                                                                 <  7:320422/138‑1 (MQ=255)
                                                          tAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCta                                                                           >  2:252907/1‑139 (MQ=255)
                                                          aaaGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCATTTAAACGTAAGAAAATCAAACCCCAAATGCGGTCTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCta                                                                           >  5:64621/2‑139 (MQ=255)
                                                                     aTCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACAATttct                                                                >  4:330659/1‑139 (MQ=255)
                                                                           aaTTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGGTTTCCAATTAAACGTAAGAAAAGCACACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGGCTATAACCATTTCTTCAAt                                                           <  4:302231/138‑1 (MQ=255)
                                                                                     agaACTTCATATTCAAGCGTCATATAATGTGGTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTGTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGa                                                <  2:318954/139‑1 (MQ=255)
                                                                                        aCTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGCAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGaca                                              <  7:121531/138‑1 (MQ=255)
                                                                                            cATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTGTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACggg                                          >  1:286391/1‑137 (MQ=255)
                                                                                             atatTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCa                                        <  5:164237/139‑1 (MQ=255)
                                                                                             atatTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCa                                        <  1:252907/139‑1 (MQ=255)
                                                                                                ttCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATa                                      <  3:40471/138‑1 (MQ=255)
                                                                                                   aaGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTc                                  >  4:184378/1‑139 (MQ=255)
                                                                                                    aGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGc                                                                                >  2:296846/1‑92 (MQ=255)
                                                                                                    aGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGc                                                                                <  1:296846/92‑1 (MQ=255)
                                                                                                    aGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCt                                 >  3:96087/1‑139 (MQ=255)
                                                                                                        tCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGt                               >  2:356335/1‑137 (MQ=255)
                                                                                                         cATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTga                                                                                                     >  4:388345/1‑66 (MQ=255)
                                                                                                         cATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTga                                                                                                     <  3:388345/66‑1 (MQ=255)
                                                                                                          atatAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGc                           >  6:144234/1‑139 (MQ=255)
                                                                                                                tGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGt                     >  1:210177/1‑139 (MQ=255)
                                                                                                                tGGGTTTTCCAATTAAACGCAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGt                     <  1:342200/139‑1 (MQ=255)
                                                                                                                  ggTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACCCGGTCATACTTCTGTAAGCTTAAGttt                   >  5:234667/1‑139 (MQ=255)
                                                                                                                  ggTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGttt                   >  1:353545/1‑139 (MQ=255)
                                                                                                                     tttCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCt                >  7:47235/1‑139 (MQ=255)
                                                                                                                      ttCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGaca                                              <  6:309972/108‑1 (MQ=255)
                                                                                                                      ttCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGaca                                              >  5:309972/1‑108 (MQ=255)
                                                                                                                      ttCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCtt               >  1:187738/1‑139 (MQ=255)
                                                                                                                      ttCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCtt               >  7:52904/1‑139 (MQ=255)
                                                                                                                      ttCCAATTAAACGCAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTACGCTTAAGTTTCCtt               >  3:289227/1‑139 (MQ=255)
                                                                                                                       tCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACCCGGTCATACTTCTGTAAGCTTAAGTTTCCtt               >  8:271525/1‑138 (MQ=255)
                                                                                                                               aaCGTAAGAAAATCACACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCTTTAATCAGc       <  4:32947/138‑1 (MQ=255)
                                                                                                                                   tAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCTTTAATCAGCTGCAg  >  2:336014/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                          
GACTGCTTGATCAAACAGATGAAACAGCTTTCTCATCAGTTTCGCCAGCTTTTGATCATAAGCCGACAAATCCTTAATTTCTCCGAGAACTTCATATTCAAGCGTCATATAATGGGTTTTCCAATTAAACGTAAGAAAATCAAACGCCAAATGCGGTGTTTTCAGGTTTGAGAGCTTTGATTCGTGCTTTGCGTCTATAACCATTTCTTCAATCATTTCAAAGACACGGTCATACTTCTGTAAGCTTAAGTTTCCTTTAATCAGCTGCAG  >  CP006881/2672039‑2672308

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