Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,868 C→A 24.0% intergenic (+30/‑129) 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,8680CA24.0% 86.5 / 23.0 46intergenic (+30/‑129)hemB/proP5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease
Reads supporting (aligned to +/- strand):  ref base C (17/18);  new base A (5/6);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

TGGCCGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAGC  >  NC_002947/3313720‑3314017
                                                                                                                                                     |                                                                                                                                                     
tggccGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCa                                                                                                                                                      >  2:857961/1‑150 (MQ=255)
      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                                              >  1:607683/1‑140 (MQ=255)
                                                                                                                   tCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTagag                                              <  2:607683/140‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
TGGCCGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTAC‑TTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAGC  >  NC_002947/3313720‑3314017

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