Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,875,792 A→G 11.1% L45P (CTG→CCG)  PP_2528 ← O‑acetylhomoserine (thiol)‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,875,7920AG11.1% 89.0 / 3.7 36L45P (CTG→CCG) PP_2528O‑acetylhomoserine (thiol)‑lyase
Reads supporting (aligned to +/- strand):  ref base A (14/18);  new base G (2/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TCTGGATGGCGTAGGTGATAGCCGCCATGCCCGAGGCTACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCATGGATGGCAAGCGTTTCCAGCTTCATGCAATCGTTCC  >  NC_002947/2875652‑2875936
                                                                                                                                            |                                                                                                                                                
tCTGGATGGCGTAGGTGATAGCCGCCATGCCCGAGGCTACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACa                                                                                                                                         >  1:556413/1‑150 (MQ=255)
                             cccGAGGCTACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTcgaag                                                                                                            >  2:239605/1‑150 (MQ=255)
                                     tACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAg                                                                                                                                                 <  1:121703/105‑1 (MQ=255)
                                     tACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAg                                                                                                                                                 >  2:121703/1‑105 (MQ=255)
                                     tACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAagga                                                                                                    >  1:639398/1‑150 (MQ=255)
                                              cgcgcCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTg                                                                                           <  2:214914/150‑1 (MQ=255)
                                                cgcCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGt                                                                                         <  1:179891/150‑1 (MQ=255)
                                                 gcCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAggt                                                                                                 <  2:54156/141‑1 (MQ=255)
                                                 gcCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAggt                                                                                                 >  1:54156/1‑141 (MQ=255)
                                                 gcCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTgg                                                                                          <  1:688403/148‑1 (MQ=255)
                                                 gcCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTgg                                                                                          >  2:688403/1‑148 (MQ=255)
                                                    aaCCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTggg                                                                                                                       <  2:248655/116‑1 (MQ=255)
                                                    aaCCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTggg                                                                                                                       >  1:248655/1‑116 (MQ=255)
                                                    aaCCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGAt                                                                                     <  1:239605/150‑1 (MQ=255)
                                                     aCCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATc                                                                                    <  2:556413/150‑1 (MQ=255)
                                                             ttCCAGGGCGGCCATGCGCTGCTCGAGCCCGTCGTTGGTGGGGTTAATGACGCGCGAATAGATGTAAGCAGCCAAATTAAGGTTACACAGGTGGGTGCTGTGGTGGGTGTCGTGGAAGGTGAAGGAGGAGGAAGGGTAGCTCGGTACGGc                                                                            >  2:130681/1‑150 (MQ=255)
                                                             ttCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCCGGGCTCTGTCGTCGCTCTCCTAGTTCTAAGGAGGTGGTCTGGTAGATCGGTACGGc                                                                            <  1:130681/150‑1 (MQ=255)
                                                                       gcCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCACCCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCACCCGCCTtg                                                                  >  1:320474/1‑150 (MQ=255)
                                                                                   tCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTggtc                                                            <  2:145090/144‑1 (MQ=255)
                                                                                   tCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTggtc                                                            >  1:145090/1‑144 (MQ=255)
                                                                                     gAGCACGGCGGTGGTGGAGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGTTGGATGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGtcgggtcggg                                                    <  2:292602/150‑1 (MQ=255)
                                                                                              gtTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTggtcaggg                                                        <  2:189245/137‑2 (MQ=255)
                                                                                              gtTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTggtcaggg                                                        >  1:189245/1‑136 (MQ=255)
                                                                                                   tggGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCat                                      <  2:639398/150‑1 (MQ=255)
                                                                                                    ggggTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatg                                     <  1:744610/150‑1 (MQ=255)
                                                                                                      ggTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatgga                                   <  1:16780/150‑1 (MQ=255)
                                                                                                        ttCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatggatg                                 >  2:408292/1‑150 (MQ=255)
                                                                                                             gatgCGCGAATAGATGTTGCCGGCAACCTTCGGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatggat                                  <  1:504773/144‑1 (MQ=255)
                                                                                                             gatgCGCGAATAGATGTTGCCGGCAACCTTCGGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatggat                                  <  1:504690/144‑1 (MQ=255)
                                                                                                             gatgCGCGAATAGATGTTGCCGGCAACCTTCGGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatggat                                  >  2:504690/1‑144 (MQ=255)
                                                                                                             gatgCGCGAATAGATGTTGCCGGCAACCTTCGGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCatggat                                  >  2:504773/1‑144 (MQ=255)
                                                                                                              atgCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGccc                                          <  2:249850/135‑1 (MQ=255)
                                                                                                              atgCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGccc                                          >  1:249850/1‑135 (MQ=255)
                                                                                                                       tAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCATGGATGGCAAGCGTTTCCAGc                  <  1:408292/150‑1 (MQ=255)
                                                                                                                          aTGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGGAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACGGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCATGGCTGGCAAGCGTTTCCAGCTTc               >  2:165674/1‑150 (MQ=255)
                                                                                                                                       ccTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCATGGATGGCAAGCGTTTCCAGCTTCATGCAATCGTTcc  <  2:738016/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
TCTGGATGGCGTAGGTGATAGCCGCCATGCCCGAGGCTACGGCCAGCGCGCCAACCCCGCCTTCCAGGGCGGCCATGCGCTGCTCGAGCACGTCGTTGGTGGGGTTCATGATGCGCGAATAGATGTTGCCGGCAACCTTCAGGTCAAACAGGTCGGCGCCGTGCTGGGTGTCGTCGAAGGCGAAGGAGGTGGTCTGGTAGATCGGTACGGCCACCGCCTTGGTGGTCGGGTCGGGGCTGAAGCCCGCATGGATGGCAAGCGTTTCCAGCTTCATGCAATCGTTCC  >  NC_002947/2875652‑2875936

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: