Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,185 G→C 5.4% A45P (GCG→CCG) ‡ ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1850GC 5.4% 2120.5 / 80.9 589A45P (GCG→CCG) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (294/263);  new base C (22/10);  total (316/273)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 600 total aligned reads displayed.

TACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGT  >  CP009273/4264119‑4264252
                                                                  |                                                                   
tACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGc                                                                    <  1:18698016/68‑1 (MQ=255)
 aCATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCg                                                                   >  1:11991341/1‑68 (MQ=255)
      ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                               >  1:11631891/1‑67 (MQ=255)
      ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                               >  1:39509157/1‑67 (MQ=255)
      ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                               >  1:16502393/1‑67 (MQ=255)
       cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                              <  1:22897044/67‑1 (MQ=255)
        aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                            >  1:13013365/1‑68 (MQ=255)
          aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                          >  1:27596945/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                         >  1:32195816/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                         >  1:32754415/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                         <  1:3336862/68‑1 (MQ=255)
            ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                        >  1:12166247/1‑68 (MQ=255)
             gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                       >  1:39314784/1‑68 (MQ=255)
             gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                       >  1:17504521/1‑68 (MQ=255)
              tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                       <  1:20779772/67‑1 (MQ=255)
              tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                       <  1:25661915/67‑1 (MQ=255)
              tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                       <  1:32405502/67‑1 (MQ=255)
               ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                     <  1:19519549/68‑1 (MQ=255)
               ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                     >  1:27219856/1‑68 (MQ=255)
                gcgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                     <  1:20330645/67‑1 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                   <  1:15833883/67‑1 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                  >  1:38620603/1‑68 (MQ=255)
                   cAGTTGCCAACATTACGCTGGCTACTTACGAATCCTGGCGTGATATAGCGACCGGCGAGATGAAAGaa                                                 <  1:19925193/68‑1 (MQ=255)
                    aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                >  1:37429647/1‑68 (MQ=255)
                    aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                >  1:20697663/1‑68 (MQ=255)
                    aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                >  1:13675886/1‑68 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                              <  1:16916219/67‑1 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                              <  1:13321311/67‑1 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                             <  1:14783848/68‑1 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                             <  1:2358398/68‑1 (MQ=255)
                         ccAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                              >  1:6177100/1‑65 (MQ=255)
                          cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                           <  1:15691488/67‑1 (MQ=255)
                          cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                           <  1:8107808/67‑1 (MQ=255)
                          cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTg                                          >  1:37222847/1‑68 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                         <  1:35558852/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                         <  1:10247858/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                         <  1:31732302/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                         <  1:6433443/68‑1 (MQ=255)
                            aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                        >  1:29236712/1‑68 (MQ=255)
                            aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                        >  1:2832752/1‑68 (MQ=255)
                            aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                        >  1:7428330/1‑68 (MQ=255)
                              aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                      >  1:26344063/1‑68 (MQ=255)
                                tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGCGATGAAAGAACAGACTGAATGGc                                    <  1:5022056/68‑1 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:17325328/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:21030707/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:4405805/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:22397958/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:34960628/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:35269315/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:8740906/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:8980834/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:22088606/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:25403906/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:7121466/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:9444169/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:37843650/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:5151179/1‑68 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:29385484/1‑67 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:7750946/1‑67 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                 >  1:10545012/1‑67 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                               <  1:11001690/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                               <  1:34409498/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                               <  1:30625252/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                               <  1:30131060/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                               <  1:28705492/68‑1 (MQ=255)
                                          aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                          <  1:13870915/68‑1 (MQ=255)
                                          aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                          <  1:1433833/68‑1 (MQ=255)
                                             tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                          >  1:24104619/1‑65 (MQ=255)
                                             tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                          >  1:16109189/1‑65 (MQ=255)
                                             tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                          >  1:31318927/1‑65 (MQ=255)
                                             tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                       <  1:19373263/68‑1 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                     <  1:21360490/68‑1 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                     <  1:25003106/68‑1 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                     <  1:21721735/68‑1 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                     >  1:30290558/1‑68 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                     <  1:17849076/68‑1 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                  >  1:4801264/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                  >  1:4199158/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                  >  1:33673883/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                  >  1:4041357/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                  >  1:30912126/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGGTGTGCTGTTc                   <  1:15475318/67‑1 (MQ=255)
                                                   agcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCgg                 <  1:9607880/66‑1 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                >  1:38864747/2‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                >  1:10078067/2‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                >  1:29848451/2‑68 (MQ=255)
                                                    ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                >  1:36830940/1‑68 (MQ=255)
                                                    ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                >  1:26699917/1‑68 (MQ=255)
                                                         cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa              >  1:22263136/1‑65 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc            >  1:34774578/1‑65 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc            <  1:12701015/67‑3 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc            <  1:27994499/67‑3 (MQ=255)
                                                          gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg          <  1:35796866/68‑1 (MQ=255)
                                                          gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg          >  1:23062370/1‑68 (MQ=255)
                                                          gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt           <  1:26002165/67‑1 (MQ=255)
                                                             aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa       >  1:18972905/1‑68 (MQ=255)
                                                              tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg      >  1:18574536/1‑68 (MQ=255)
                                                                 aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg   >  1:2452551/1‑68 (MQ=255)
                                                                  gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGGAAACTGGCAGAAGt  >  1:36425123/1‑68 (MQ=255)
                                                                  gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt  >  1:5514144/1‑68 (MQ=255)
                                                                  |                                                                   
TACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGT  >  CP009273/4264119‑4264252

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