Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,115,408 T→C 13.3% E444G (GAA→GGA)  priA ← Primosome factor n' (replication factor Y)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,115,4080TC13.3% 69.8 / 4.2 30E444G (GAA→GGA) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base C (2/2);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

GGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCAGACCTGCCTGCACCTTCTGACCTTT  >  CP009273/4115266‑4115556
                                                                                                                                              |                                                                                                                                                    
ggggCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCaa                                                                                                                                                <  1:354748/146‑1 (MQ=255)
               ggggACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGt                                                                                                                                 >  2:319737/1‑149 (MQ=255)
                   aCCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGca                                                                                                                             <  2:295094/149‑1 (MQ=255)
                           cGTGGAACCGCAGGAGGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGGGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTTGGCATTCGGCAATCCAGCCACAGTCGTGGCGCAGCAGTg                                                                                                                     >  1:166990/1‑149 (MQ=255)
                              ggAACCGCAGGAAGGGGACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAg                                                                                                                  <  2:51954/149‑1 (MQ=255)
                                  ccGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGgcgc                                                                                                              >  2:341075/1‑149 (MQ=255)
                                  ccGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGgcgc                                                                                                              <  1:169731/149‑1 (MQ=255)
                                  ccGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATCCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGgcgc                                                                                                              >  1:178408/1‑149 (MQ=255)
                                  ccGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATCCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGgcgc                                                                                                              >  1:241227/1‑149 (MQ=255)
                                   cGCAGGAAGGGCACTGGCGCGGCACCGGACGCTAGCTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCCTTCGGCAATCCAGCCACAGGCGTGGCACCGCAGTGCAGgcgcg                                                                                                             >  2:230053/1‑148 (MQ=255)
                                    gCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCaa                                                                                                            <  2:136970/149‑1 (MQ=255)
                                          aGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAc                                                                                                      >  2:274637/1‑149 (MQ=255)
                                                    cgcgGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCGTGATGCAGCGTGAAGCTGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTaag                                                                                            >  2:270306/1‑149 (MQ=255)
                                                         cACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTaagaaaga                                                                                       <  1:131626/149‑1 (MQ=255)
                                                            cGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGaa                                                                                    <  1:96628/149‑1 (MQ=255)
                                                             ggACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAAt                                                                                   <  1:319737/149‑1 (MQ=255)
                                                                                 gtggCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTCGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAACGCGGTTAAGAAAGAGAAGGACCTAGTTCTCAGCCTGTa                                                               >  1:106019/1‑149 (MQ=255)
                                                                                      ggcgcAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGc                                                          <  1:45931/148‑1 (MQ=255)
                                                                                         gcAGATGGGGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACCGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTgg                                                       <  1:179757/149‑1 (MQ=255)
                                                                                          cAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGgc                                                      <  1:217263/149‑1 (MQ=255)
                                                                                                 tgctgcGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACg                                               >  2:40595/1‑149 (MQ=255)
                                                                                                           gATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAAc                                     <  2:20488/149‑1 (MQ=255)
                                                                                                           gATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAAc                                     >  2:326110/1‑149 (MQ=255)
                                                                                                                    tgtAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGGGcc                            >  2:285656/1‑149 (MQ=255)
                                                                                                                     gtagtaGTGATCGCAACGTGGGCATCCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCa                           <  2:241227/149‑1 (MQ=255)
                                                                                                                     gtagtaGTGATCGCAACGTGGGCATCCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCa                           <  2:178408/149‑1 (MQ=255)
                                                                                                                          agtGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCAGANGAGNGNTTANCNCCGGCGCCAGGcct                      >  2:190678/1‑149 (MQ=255)
                                                                                                                           gtGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCAGAcctg                     >  1:132948/1‑149 (MQ=255)
                                                                                                                              aTCGCAACGTGGGCATTCGGCAAACCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCAGAcctgcct                  >  2:188008/1‑149 (MQ=255)
                                                                                                                                  cAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCAGACCTGCCTGCAc              <  2:154808/149‑1 (MQ=255)
                                                                                                                                              tCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCAGACCTGCCTGCACCTTCTGACCttt  <  2:145888/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGATTAACGCCGGAGCCAGACCTGCCTGCACCTTCTGACCTTT  >  CP009273/4115266‑4115556

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