Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,871 T→G 24.0% intergenic (+33/‑126) hemB → / → proP 5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,313,8710TG24.0% 67.8 / 26.7 46intergenic (+33/‑126)hemB/proP5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease
Reads supporting (aligned to +/- strand):  ref base T (16/19);  new base G (5/6);  total (21/25)
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.11e-01

GTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAGC  >  NC_002947/3313726‑3314017
                                                                                                                                                  |                                                                                                                                                  
ttGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGTCCTTGGTCATGATAGCCGCTGCCATTGCAg                                                                                                                                                <  2:126569/149‑1 (MQ=255)
 tGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGc                                                                                                                                               <  2:561712/150‑1 (MQ=255)
    cATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGc                                                                                                                                            <  1:399904/150‑1 (MQ=255)
    cATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGc                                                                                                                                            <  1:290721/150‑1 (MQ=255)
          cgaAGGCCGCGTAATGGCCCAAAGCGCTGTCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTG‑CCTGCTTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCt                                                                                                                                      <  1:265813/150‑1 (MQ=255)
           gaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttc                                                                                                                                        <  2:1132333/147‑2 (MQ=255)
           gaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttc                                                                                                                                        >  1:1132333/1‑146 (MQ=255)
            aaGGCCGGGTGGTGCGCGAATGCATTGGAGCAAGCAAGCGGGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCttt                                                                                                                                    <  2:615100/150‑1 (MQ=255)
             aGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTtg                                                                                                                                   <  1:1015875/150‑1 (MQ=255)
                         gcgcGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGCGCGCCGTTggg                                                                                                                       >  1:1022181/1‑150 (MQ=255)
                           gcgAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGggtg                                                                                                                     >  2:385603/1‑150 (MQ=255)
                                    ttGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGTACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCt                                                                                                                                      >  2:550546/1‑124 (MQ=255)
                                    ttGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGTACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCt                                                                                                                                      <  1:550546/124‑1 (MQ=255)
                                     tGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCaaa                                                                                                           <  1:1080307/150‑1 (MQ=255)
                                                          tGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAgcg                                                                                      <  2:682135/150‑1 (MQ=255)
                                                              gACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGa                                                                                  >  1:135300/1‑150 (MQ=255)
                                                                        tGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtt                                                                                                                                          >  2:909244/1‑84 (MQ=255)
                                                                        tGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtt                                                                                                                                          <  1:909244/84‑1 (MQ=255)
                                                                          aCCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAgtgt                                                                      >  1:571963/1‑150 (MQ=255)
                                                                                      cATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATc                                                           >  2:241130/1‑150 (MQ=255)
                                                                                       aTGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCt                                                          <  1:857961/150‑1 (MQ=255)
                                                                                       aTGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat                                                                                                                                      <  2:406358/74‑4 (MQ=37)
                                                                                       aTGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat                                                                                                                                      >  1:406358/1‑71 (MQ=37)
                                                                                          attccggCGCGCGAAGGTATATGAAGCCAGGCATTGGTCATGAAAGCCGCAGCCATTGCATCGGCTTTTCTTAGTAATCCGTTGGGTGATGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGAt                                                       <  2:236389/145‑1 (MQ=255)
                                                                                                   cgcgAAGGTATCTGAATCCAGGCCTTGGTCATGAAAGCCGCTTCCATTGCAGCGGCTATTCTTTGTGCGTCGTTCGCGGGTGCCACAAACGCACAAAAAAACAGTAAGCGCTGGATTATTTAGTGTCCTGTTGTATCTGATCCCTTagag                                              >  2:644274/1‑150 (MQ=255)
                                                                                                     cgAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGATGTTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCa                                            <  1:385603/150‑1 (MQ=255)
                                                                                                      gAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGAGAAGACAAAGCCGAAAAAAATAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTAACTGATCCCTTAGAGCAt                                           >  1:142128/1‑150 (MQ=255)
                                                                                                            aTCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCg                                     >  2:1155692/1‑150 (MQ=255)
                                                                                                             tCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTagag                                              <  2:607683/140‑1 (MQ=255)
                                                                                                             tCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTagag                                              >  1:607683/1‑140 (MQ=255)
                                                                                                                gAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAg                                 >  1:1179609/1‑150 (MQ=255)
                                                                                                                      aGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTAATGCCGCAAACGACGCAATAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATAc                           >  1:749150/1‑150 (MQ=255)
                                                                                                                        gCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACgg                         >  2:679160/1‑150 (MQ=255)
                                                                                                                              caaagaaaaGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGa                   >  2:387499/7‑150 (MQ=255)
                                                                                                                                 aTGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAAc                >  1:624068/1‑150 (MQ=255)
                                                                                                                                  tGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgttgt                                                              <  2:770045/103‑1 (MQ=255)
                                                                                                                                  tGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgttgt                                                              >  1:770045/1‑103 (MQ=255)
                                                                                                                                    aaaGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTc             <  2:1179609/150‑1 (MQ=255)
                                                                                                                                    aaaGCCGCTGCCATTGCAGCGGCTTTTCTTTGCGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTc             <  2:1022181/150‑1 (MQ=255)
                                                                                                                                       gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGTTGCCGCAAACGCACAAAAAAACa                                                                                             <  1:337333/67‑1 (MQ=37)
                                                                                                                                       gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACa                                                                                             >  2:337333/1‑67 (MQ=38)
                                                                                                                                       gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCa          >  2:1089206/1‑150 (MQ=255)
                                                                                                                                        ccGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCaa         <  2:571963/150‑1 (MQ=255)
                                                                                                                                            tGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCAGGGACAGGGATACGGGTATGAAACCTCGCAAGaaa     <  1:400621/150‑1 (MQ=255)
                                                                                                                                             gCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGaaaa    <  1:241130/150‑1 (MQ=255)
                                                                                                                                               cATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAGc  <  1:1155692/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAGC  >  NC_002947/3313726‑3314017

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