Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,734,098 A→G 11.4% intergenic (‑68/+89) proY ← / ← hutH proline (histidine) APC transporter/histidine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,734,0980AG11.4% 85.3 / 3.4 35intergenic (‑68/+89)proY/hutHproline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (16/15);  new base G (2/2);  total (18/17)
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.97e-01

CATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAG  >  NC_002947/5733974‑5734246
                                                                                                                            |                                                                                                                                                    
cATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCGTGCCCGcgaagg                                                                                                                                          >  2:279063/1‑132 (MQ=255)
cATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCGTGCCCGcgaagg                                                                                                                                          <  1:279063/137‑6 (MQ=255)
   gAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtttt                                                                                                                          <  1:168536/150‑1 (MQ=255)
     aaCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGc                                                                                                                        >  1:294068/1‑150 (MQ=255)
               gCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCCCCTACTTCt                                                                                                              >  1:257921/1‑150 (MQ=255)
                 gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                            <  2:393117/150‑1 (MQ=255)
                 gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                            <  1:555379/150‑1 (MQ=255)
                   gCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTgg                                                                                                          >  2:131799/1‑150 (MQ=255)
                                       aCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGTGGGCGCGCCCGCGAAg                                                                                      >  2:404914/1‑150 (MQ=255)
                                         cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAggg                                                                                    >  1:98114/1‑150 (MQ=255)
                                               cTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAgg                                                                                     >  2:370156/1‑143 (MQ=255)
                                               cTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAgg                                                                                     <  1:370156/143‑1 (MQ=255)
                                                            aCGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAggg                                                                 >  2:217555/1‑150 (MQ=255)
                                                            aCGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTAGCGGGCGTGCCCGc                                                                                                                                               <  2:480892/72‑1 (MQ=38)
                                                            aCGTCCTTGAATTGTTGTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTAGCGGGCGTGCCCGc                                                                                                                                               >  1:480892/1‑72 (MQ=38)
                                                              gTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGc                                                               >  2:148196/1‑150 (MQ=255)
                                                                ccTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGt                                                             <  2:257921/150‑1 (MQ=255)
                                                                       ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGc                                                      >  1:118829/1‑150 (MQ=255)
                                                                             ttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGTGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                <  1:404914/150‑1 (MQ=255)
                                                                                tgGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTa                                             <  1:217555/150‑1 (MQ=255)
                                                                                     gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                                                  <  2:209193/60‑1 (MQ=255)
                                                                                     gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                                                  >  1:209193/1‑60 (MQ=255)
                                                                                     gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAccc                                        >  1:425223/1‑150 (MQ=255)
                                                                                     gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAccc                                        <  2:294068/150‑1 (MQ=255)
                                                                                            tCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGgca                                 >  1:427781/1‑150 (MQ=255)
                                                                                              aGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGgcagc                               >  2:262877/1‑150 (MQ=255)
                                                                                                        tgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCa                     <  2:81569/150‑1 (MQ=255)
                                                                                                          cgcAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGa                   <  1:148196/150‑1 (MQ=255)
                                                                                                             aGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTCCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGCAGGGGCGCCCAGCGGCCCCCGGGGCGTTAAAGGCTTGGCAGTACACCCCCCCGCAGCAAGCCGGTGCGACtc                >  2:604729/1‑149 (MQ=255)
                                                                                                              gCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTAc               <  2:425223/150‑1 (MQ=255)
                                                                                                                 ccTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                            <  2:25381/54‑1 (MQ=255)
                                                                                                                 ccTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                            >  1:25381/1‑54 (MQ=255)
                                                                                                                  cTTCGCGGGCACGCCCGCTCCCACAGGGGTTGGGGTTTTGCACCTACTTCTGTGGGAGCGGGCTCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCtt           >  2:470668/1‑150 (MQ=255)
                                                                                                                     cgcgGGCACGCCCGCTCCCACGGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGc        <  1:604729/150‑1 (MQ=255)
                                                                                                                           cACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAg  <  2:427781/150‑1 (MQ=255)
                                                                                                                            |                                                                                                                                                    
CATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAG  >  NC_002947/5733974‑5734246

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