Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,627,469 T→C 13.3% K104R (AAG→AGG)  U712_18705 ← Sporulation initiation phosphotransferase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,627,4690TC13.3% 148.0 / 16.5 61K104R (AAG→AGG) U712_18705Sporulation initiation phosphotransferase F
Reads supporting (aligned to +/- strand):  ref base T (20/32);  new base C (2/6);  total (22/39)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-02

ATAACTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATGACCTTCATCCGTTTTAAGATTTCGATTCCGTC  >  CP006881/3627336‑3627599
                                                                                                                                     |                                                                                                                                  
aTAACTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCa                                                                                                                               <  3:202207/139‑1 (MQ=255)
aTAACTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCa                                                                                                                               <  1:176168/139‑1 (MQ=255)
 tAACTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCa                                                                                                                               >  4:50954/1‑138 (MQ=255)
  aaCTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCaaa                                                                                                                             <  6:248382/139‑1 (MQ=255)
   aCTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCaaa                                                                                                                             <  2:288154/138‑1 (MQ=255)
             cAACTTTTCGACAAATCGGCAACATAATGACATGCTTCATTTGTTTCTTCTTTTTGTCAGATAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCagag                                                                                                                  <  8:206163/139‑1 (MQ=255)
             cAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCaaa                                                                                                                             >  4:350362/1‑128 (MQ=255)
             cAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCaaa                                                                                                                             <  3:350362/128‑1 (MQ=255)
             cAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGGCGATGTCAAACGGCTTGGCAAAGTGCGTCagag                                                                                                                  >  3:39920/1‑139 (MQ=255)
               aCTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGAGTCTCTTATTTAGTCGATGTCAAAAGGCTTGGCAAAGTGCGTCAGAgcg                                                                                                                >  7:200491/1‑139 (MQ=255)
                  tttCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGcc                                                                                                              >  3:327163/1‑138 (MQ=255)
                     cGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCaa                                                                                                           <  3:363275/138‑1 (MQ=255)
                          aaTCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCAAAGTGCGTCAGAGCGCCCAATTCCtt                                                                                                     <  2:310242/139‑1 (MQ=255)
                                cAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCTGGGGCAGTTATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTc                                                                                               <  4:310948/139‑1 (MQ=255)
                                    aTAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTg                                                                                            >  7:253117/1‑138 (MQ=255)
                                     tAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGa                                                                                          <  4:327163/139‑1 (MQ=255)
                                     tAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGa                                                                                          <  4:116502/139‑1 (MQ=255)
                                       aTGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATc                                                                                        <  6:350022/139‑1 (MQ=255)
                                       aTGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATc                                                                                        <  6:86948/139‑1 (MQ=255)
                                           cATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGt                                                                                    <  3:128422/139‑1 (MQ=255)
                                             tGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTAGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATACCTTCGATTCCTGGATCATGTc                                                                                   >  7:195816/1‑138 (MQ=255)
                                              gaatCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTc                                                                                   >  3:238931/4‑137 (MQ=255)
                                              gATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTc                                                                                   >  4:167842/1‑137 (MQ=255)
                                              gATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTAGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGCTTCCTGGATCATGTc                                                                                   >  5:76974/1‑137 (MQ=255)
                                               aTTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGCTGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCGGGTTCATGTagag                                                                                >  5:245364/1‑139 (MQ=255)
                                               aTTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAg                                                                                <  3:79365/139‑1 (MQ=255)
                                               aTTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAg                                                                                <  4:156072/139‑1 (MQ=255)
                                               aTTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGGCGATGTCAAACGGCTTGGCAAAGTGCGTCCGCGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAg                                                                                <  7:295132/139‑1 (MQ=255)
                                                ttCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGc                                                                               <  4:236341/139‑1 (MQ=255)
                                                  cATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCa                                                                                                                     >  4:144434/1‑99 (MQ=255)
                                                  cATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCa                                                                                                                     <  3:144434/99‑1 (MQ=255)
                                                   aTTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGctc                                                                             <  2:206456/138‑1 (MQ=255)
                                                        tttCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGt                                                                        <  1:203983/138‑1 (MQ=255)
                                                         ttcttcTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTc                                                                                                    >  5:101769/1‑109 (MQ=255)
                                                         ttcttcTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTc                                                                                                    >  5:107213/1‑109 (MQ=255)
                                                         ttcttcTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTc                                                                                                    <  6:107213/109‑1 (MQ=255)
                                                         ttcttcTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTc                                                                                                    <  6:101769/109‑1 (MQ=255)
                                                          tcttctTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCa                                                                                                                               <  8:87398/81‑1 (MQ=255)
                                                          tcttctTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCa                                                                                                                               >  7:87398/1‑81 (MQ=255)
                                                          tcttctTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCa                                                                                                                               <  8:253527/81‑1 (MQ=255)
                                                          tcttctTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCCTGGCa                                                                                                                               >  7:253527/1‑81 (MQ=255)
                                                              cttTTTGTCAGTTAGACTTCTGGGGCATATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGc                                                                    <  7:330064/136‑1 (MQ=255)
                                                                 tttGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCa                                                               <  1:257794/138‑1 (MQ=255)
                                                                 tttGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCa                                                               <  8:287894/138‑1 (MQ=255)
                                                                 tttGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCa                                                               <  7:206072/138‑1 (MQ=255)
                                                                  ttGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTACGTCGATGTCAAACGGCTTGGCAAAGTGCCTCAGAGCGCCCAAATCCTTCCCTTCCGGGATCATGTCGGGCTCTGCGTATGCCGTc                                                                >  7:134712/1‑136 (MQ=255)
                                                                       aGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATaa                                                         <  3:167842/138‑1 (MQ=255)
                                                                           aGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCGTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGAc                                                     <  6:123390/138‑1 (MQ=255)
                                                                              cTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCgg                                                 <  3:190754/139‑1 (MQ=255)
                                                                                         gATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGc                                                                               >  2:33788/1‑98 (MQ=255)
                                                                                         gATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGc                                                                               <  1:33788/98‑1 (MQ=255)
                                                                                         gATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGt                                       >  8:275019/1‑138 (MQ=255)
                                                                                           tatTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCaa                                    <  7:270968/139‑1 (MQ=255)
                                                                                           tatTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCaa                                    <  8:253117/139‑1 (MQ=255)
                                                                                             ttttttGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATg                                  <  3:94786/139‑1 (MQ=255)
                                                                                                        cGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCa                                                                                                                     <  3:4882/45‑1 (MQ=255)
                                                                                                        cGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCa                                                                                                                     >  4:4882/1‑45 (MQ=255)
                                                                                                                     gtcgATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATGACCTTCATCCGTTTTAAGATTTc           >  2:137603/1‑138 (MQ=255)
                                                                                                                     gtcgATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATGACCTTCATCCGTTTTAAGATTTc           >  6:357063/1‑138 (MQ=255)
                                                                                                                        gATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATGACCTTCATCCGTTTTAAGATTTCGAtt       <  1:137603/139‑1 (MQ=255)
                                                                                                                              aaccGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCGGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATGACCTTCCTCCGTTTTAAGATTTCGATTCCGTc  <  5:302205/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
ATAACTAGGATACCAACTTTTCGACAAATCGGCAACATAATGACATGATTCATTTGTTTCTTCTTTTTGTCAGTTAGACTTCAGGGGCAGATATTTTTTGACGGCGTCTCTGATTTCGTCGATGTCAAACGGCTTGGCAAAGTGCGTCAGAGCGCCCAATTCCTTCGATTCCTGGATCATGTCGAGCTCTCCGTATGCCGTCATGATAATGACCCGGATGTTTTCGTCAATGACCTTCATCCGTTTTAAGATTTCGATTCCGTC  >  CP006881/3627336‑3627599

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