Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,236,491 A→C F214C (TTT→TGT)  folE ← GTP cyclohydrolase I FolE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,236,4910AC98.0% 171.7 / ‑5.2 49F214C (TTT→TGT) folEGTP cyclohydrolase I FolE
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (23/25);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GCGACGCCGCGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAC  >  NZ_CP009273/2236406‑2236575
                                                                                     |                                                                                    
gCGACGCCGCGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTc                                                                                  >  2:1808209/1‑90 (MQ=255)
   aCGCCGCGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTg                                                                               <  2:613227/90‑1 (MQ=255)
   aCGCCGCGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTg                                                                               <  2:296111/90‑1 (MQ=255)
     gccgcGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGgc                                                                             >  1:1227411/1‑90 (MQ=255)
     gccgcGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGgc                                                                             >  1:1720758/1‑90 (MQ=255)
       cgcgAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGgcgc                                                                           <  1:1527981/90‑1 (MQ=255)
           aaCAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTAt                                                                       <  2:1379747/90‑1 (MQ=255)
                       cGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGAtt                                                           >  2:1477364/1‑90 (MQ=255)
                         tGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGAGCGTATTCTGACTGGATTTg                                                         >  2:1901337/1‑90 (MQ=255)
                            cGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAAc                                                      <  1:850462/90‑1 (MQ=255)
                             gTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACa                                                     <  2:1887817/90‑1 (MQ=255)
                               tGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAAt                                                   <  2:1227411/90‑1 (MQ=255)
                               tGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAAt                                                   <  1:1673797/90‑1 (MQ=255)
                               tGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAAt                                                   <  2:813601/90‑1 (MQ=255)
                               tGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAAt                                                   <  1:328838/90‑1 (MQ=255)
                                  gcTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATcca                                                <  2:1195649/90‑1 (MQ=255)
                                  gcTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATcca                                                <  2:1279510/90‑1 (MQ=255)
                                     ccATGGTTCCTGCCTTTTAATCAGTTGTGATTACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATccacca                                             >  1:1841481/1‑90 (MQ=255)
                                          gTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACt                                                                >  2:1650737/1‑66 (MQ=255)
                                          gTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACt                                                                <  1:1650737/66‑1 (MQ=255)
                                          gTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACa                                                     <  1:1986094/77‑1 (MQ=255)
                                          gTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACa                                                     >  2:1986094/1‑77 (MQ=255)
                                             ccTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGt                                     >  2:747257/1‑90 (MQ=255)
                                             ccTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGt                                     >  2:2471347/1‑90 (MQ=255)
                                               tGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCg                                   >  2:249175/1‑90 (MQ=255)
                                               cagcTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCg                                   >  2:1171980/4‑90 (MQ=255)
                                               cagcTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCg                                   >  2:1306783/4‑90 (MQ=255)
                                                gCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGt                                  <  1:533798/90‑1 (MQ=255)
                                                gCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGt                                  <  1:443637/90‑1 (MQ=255)
                                                 ccTTTTAATCAGTTGTGATTACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGtt                                 <  2:1841481/90‑1 (MQ=255)
                                                   ttttAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGt                               >  1:2465741/1‑90 (MQ=255)
                                                    tttAATCAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTc                              >  1:495016/1‑90 (MQ=255)
                                                          cAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCa                        <  1:2251581/90‑1 (MQ=255)
                                                          cAGTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCa                        <  1:747257/90‑1 (MQ=255)
                                                            gTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATcca                                                >  1:1717123/1‑64 (MQ=255)
                                                            gTTGTGATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATcca                                                <  2:1717123/64‑1 (MQ=255)
                                                              tgtgATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAagag                                         <  2:1657543/69‑1 (MQ=255)
                                                              tgtgATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAagag                                         <  2:2190000/69‑1 (MQ=255)
                                                              tgtgATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAagag                                         >  1:1657543/1‑69 (MQ=255)
                                                              tgtgATGACGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAagag                                         >  1:2190000/1‑69 (MQ=255)
                                                                 gatgaCGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACt                      >  2:1649168/1‑85 (MQ=255)
                                                                 gatgaCGCACAGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACt                      <  1:1649168/85‑1 (MQ=255)
                                                                          caGCGCGCAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGa        >  1:752224/1‑90 (MQ=255)
                                                                            gcgcgcAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATg      >  2:23448/1‑90 (MQ=255)
                                                                              gcgcAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAAt                                                   <  2:2184204/43‑1 (MQ=255)
                                                                              gcgcAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAAt                                                   >  1:2184204/1‑43 (MQ=255)
                                                                              gcgcAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGcc    >  1:2024123/1‑90 (MQ=255)
                                                                              gcgcAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGcc    >  2:829624/1‑90 (MQ=255)
                                                                                gcAGACACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAc  <  1:2128856/90‑1 (MQ=255)
                                                                                     |                                                                                    
GCGACGCCGCGAACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAC  >  NZ_CP009273/2236406‑2236575

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