Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,868 C→A 18.2% 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,8680CA18.2% 102.9 / 13.8 44intergenic (+30/‑129)hemB/proP5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease
Reads supporting (aligned to +/- strand):  ref base C (14/22);  new base A (4/4);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

TGGCCGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAG  >  NC_002947/3313720‑3314016
                                                                                                                                                    |                                                                                                                                                    
tggccGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCa                                                                                                                                                     >  1:568653/1‑150 (MQ=255)
       tGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGc                                                                                                                                              <  1:550874/150‑1 (MQ=255)
           aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCt                                                                                                                                          >  2:122007/1‑150 (MQ=255)
             cgacgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCttt                                                                                                                                        <  1:252193/150‑1 (MQ=255)
              gacgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCtttt                                                                                                                                       <  1:1260800/150‑1 (MQ=255)
                cgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat                                                                                                                                     <  1:46844/150‑4 (MQ=255)
                cgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat                                                                                                                                     <  2:87310/150‑4 (MQ=255)
                    ggCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTtgt                                                                                                                                 <  2:878666/150‑1 (MQ=255)
                       cGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTgcg                                                                                                                              >  2:198886/1‑150 (MQ=255)
                           tggtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTACGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGTGCCGt                                                                                                                          <  2:23165/150‑1 (MQ=255)
                                    aaaGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGgtggtg                                                                                                                 >  2:680207/1‑150 (MQ=255)
                                     aaGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGc                                                                                                                <  2:540123/150‑1 (MQ=255)
                                         aTTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCa                                                                                                            >  1:491542/1‑150 (MQ=255)
                                            ggCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAAc                                                                                                         <  1:2444/150‑1 (MQ=255)
                                             gcgcAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACg                                                                                                        <  2:239798/150‑1 (MQ=255)
                                              cgcAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGc                                                                                                       >  1:1251225/1‑150 (MQ=255)
                                                       gCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTCCGCCGTTGGGTGGTGCCGCAAACGCACaaaaaaa                                                                                              >  1:514068/1‑150 (MQ=255)
                                                         gTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACa                                                                                            <  2:514068/150‑1 (MQ=255)
                                                          tGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAg                                                                                           <  2:568653/150‑1 (MQ=255)
                                                           gCTGGTGCAGACCTGATCCTGACCTACTTCGCCATAGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGt                                                                                          <  1:680207/150‑1 (MQ=255)
                                                            cTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTa                                                                                         >  2:533387/1‑150 (MQ=255)
                                                                                aCCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAgtgt                                                                     >  2:991775/1‑150 (MQ=255)
                                                                                      ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                               <  1:484103/150‑1 (MQ=255)
                                                                                      ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                               <  2:491542/150‑1 (MQ=255)
                                                                                      ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTAATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                               <  2:1115661/150‑1 (MQ=255)
                                                                                             ttgACGTGGCGCGCGAAGGTATCTGCAGCCAGGTGTTGGTCACGAAAGCCGCTGCCATTGCAGCGGCGGTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTg                                                        <  1:536778/148‑1 (MQ=255)
                                                                                                aCCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATc                                                     >  2:508718/1‑150 (MQ=255)
                                                                                                         gcgAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGc                                            >  1:952921/1‑150 (MQ=255)
                                                                                                            aaGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCAtt                                         >  2:50705/1‑150 (MQ=255)
                                                                                                             aGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTc                                        <  1:448079/150‑1 (MQ=255)
                                                                                                               gTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCAt                                      <  1:396565/150‑1 (MQ=255)
                                                                                                                tATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATc                                     <  1:50705/150‑1 (MQ=255)
                                                                                                                    tGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACa                                 <  1:991775/150‑1 (MQ=255)
                                                                                                                          cAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATa                           >  2:927481/1‑150 (MQ=255)
                                                                                                                           aGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATAc                          <  2:952921/150‑1 (MQ=255)
                                                                                                                               cTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACggat                      <  1:533387/150‑3 (MQ=255)
                                                                                                                                   gTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGa                  >  2:1196035/1‑150 (MQ=255)
                                                                                                                                    aaagaaaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa                                                >  1:934079/6‑119 (MQ=39)
                                                                                                                                    aaagaaaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa                                                <  2:934079/114‑1 (MQ=39)
                                                                                                                                            gCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAAc                                                                                             >  2:519559/1‑66 (MQ=255)
                                                                                                                                            gCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAAc                                                                                             <  1:519559/66‑1 (MQ=255)
                                                                                                                                            gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCa         >  1:1088519/1‑150 (MQ=255)
                                                                                                                                             ccGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAggg                              <  2:11010/128‑1 (MQ=255)
                                                                                                                                             ccGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAggg                              >  1:11010/1‑128 (MQ=255)
                                                                                                                                                   ccATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAg  <  1:508718/150‑1 (MQ=255)
                                                                                                                                                   ccATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAg  >  1:26079/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCCGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAG  >  NC_002947/3313720‑3314016

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