Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,381,445 T→G 20.8% intergenic (+31/+56) oprF → / ← cobA porin F/uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,381,4450TG20.8% 50.8 / 12.1 33intergenic (+31/+56)oprF/cobAporin F/uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (12/14);  new base G (3/4);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

CGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCA  >  NC_002947/2381301‑2381585
                                                                                                                                                |                                                                                                                                            
cGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGcc                                                                                                                                         <  1:360999/150‑1 (MQ=255)
   gtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCggg                                                                                                                                      <  1:706465/150‑1 (MQ=255)
    ttGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGt                                                                                                                                     <  1:251680/150‑1 (MQ=255)
          cGGTTGGCTCCGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGTTTTTca                                                                                                                               <  1:650451/150‑2 (MQ=255)
              tGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGACGAAGTGCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGc                                                                                                                           <  1:1037008/150‑1 (MQ=255)
                  tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAACCCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                       <  2:741934/150‑1 (MQ=255)
                            cccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGCACCACCCCGCCTCTGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg                                                                                                             <  1:834144/150‑1 (MQ=255)
                               gTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGAAACCATGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTACATTAGAAACCCGGCCTCGGCCGGGTTTTCCTTTGCCTGGGTTTTAACCGGcc                                                                                                          >  1:822315/1‑150 (MQ=255)
                                         tGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCg                                                                                                >  1:459479/1‑150 (MQ=255)
                                                  cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCgg                                                                                       >  2:939754/1‑150 (MQ=255)
                                                                      tCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATcc                                                                   >  1:576929/1‑150 (MQ=255)
                                                                             tCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                     >  1:518567/1‑109 (MQ=255)
                                                                             tCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                     <  2:518567/109‑1 (MQ=255)
                                                                                     aaGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgct                                                    >  2:843681/1‑150 (MQ=255)
                                                                                         tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGGCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTTCGATTGCGATCCGCCCCATAGCCTGCTGCTc                                                >  1:330383/1‑150 (MQ=255)
                                                                                         tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCAAGAAAAACCCGGCCGAGGCCGGGtt                                                                                                                                    >  1:678136/1‑66 (MQ=25)
                                                                                         tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCAAGAAAAACCCGGCCGAGGCCGGGtt                                                                                                                                    <  2:678136/66‑1 (MQ=25)
                                                                                                     ccAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCct                                                       <  1:764533/131‑1 (MQ=255)
                                                                                                     ccAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCct                                                       >  2:764533/1‑131 (MQ=255)
                                                                                                                      tttGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATc                             >  1:655463/1‑140 (MQ=255)
                                                                                                                      tttGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATc                             <  2:655463/140‑1 (MQ=255)
                                                                                                                      tttGAAGCTTCATGAAAAACCCGGCCTCGACCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAAGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAAGATTg                   >  1:1042107/1‑150 (MQ=255)
                                                                                                                       ttGAAGCTTAATGGAACACGTGCCCGCGGCCGGGTTTTCCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGc                  <  2:822315/150‑1 (MQ=255)
                                                                                                                        tGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGcc                 <  2:848763/150‑1 (MQ=255)
                                                                                                                               tCATGAAAACCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAAAGGGAACTGGAGCCCGCTCCTGCGGTCGGCTCGTGATTGCGATCATCCCCATAGTCTTCTGCTCCCCAACTACCTCCCCAATCACCAGGATTGCCGGGCTAc          <  2:1042107/150‑1 (MQ=255)
                                                                                                                                    aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCt            <  1:1107943/143‑1 (MQ=255)
                                                                                                                                    aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCt            >  2:1107943/1‑143 (MQ=255)
                                                                                                                                    aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa     <  2:459479/150‑1 (MQ=255)
                                                                                                                                    aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa     <  2:576929/150‑1 (MQ=255)
                                                                                                                                       aaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTa           <  2:1124784/141‑1 (MQ=255)
                                                                                                                                       aaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTa           >  1:1124784/1‑141 (MQ=255)
                                                                                                                                       aaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCCGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaacca  >  2:1159465/1‑150 (MQ=255)
                                                                                                                                       aaCCCGGCCGAGACCGGGTTTTTCTTTGCCTGGGTTATAACCGGCCCTGAAGGCCGCTCCTGCGGTCATCCTGCACTTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCCCCTCCCCTATCCCCAGGATTGCCGGGCTACCcaaacca  >  2:181064/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                            
CGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCA  >  NC_002947/2381301‑2381585

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