Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,063,627 G→A A50V (GCT→GTT)  cbpA ← curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,063,6270GA97.0% 104.3 / ‑4.0 33A50V (GCT→GTT) cbpAcurved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/12);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAA  >  NC_000913/1063490‑1063761
                                                                                                                                         |                                                                                                                                      
cGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACg                                                                                                                                      >  2:484612/1‑140 (MQ=255)
      aaTCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTCAACACTTCCCAGGCTTCAACGACctct                                                                                                                                >  1:316122/1‑140 (MQ=255)
           tCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTCAACACTTCCCAGGCTTCAACGACCTCTGTGaa                                                                                                                           <  2:316122/140‑1 (MQ=255)
                           gACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATcc                                                                                                           >  2:189273/1‑140 (MQ=255)
                           gACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATcc                                                                                                           >  1:402585/1‑140 (MQ=255)
                            aCCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCg                                                                                                          >  1:391470/1‑140 (MQ=255)
                                    cATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTg                                                                                                  <  1:214027/140‑1 (MQ=255)
                                      tggtggAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCt                                                                                                <  1:287291/140‑1 (MQ=255)
                                           gAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACAt                                                                                           >  2:62317/1‑140 (MQ=255)
                                                   ggTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATgg                                                                                   >  2:244977/1‑140 (MQ=255)
                                                   ggTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATgg                                                                                   >  2:493010/1‑140 (MQ=255)
                                                      tAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTAt                                                                                >  1:115424/1‑140 (MQ=255)
                                                             cGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCggg                                                                         <  1:244977/140‑1 (MQ=255)
                                                              ggATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGc                                                                        <  1:62317/140‑1 (MQ=255)
                                                                 tCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAg                                                                     <  2:210168/140‑1 (MQ=255)
                                                                            ttGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAgg                                                          >  2:425480/1‑140 (MQ=255)
                                                                                  caTCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCt                                                    <  2:402585/140‑1 (MQ=255)
                                                                                     cTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGa                                                 <  1:293580/140‑1 (MQ=255)
                                                                                     cTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGa                                                 <  1:32993/140‑1 (MQ=255)
                                                                                             aCTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTg                                         >  1:213154/1‑140 (MQ=255)
                                                                                                 aGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAt                                     >  1:275088/1‑140 (MQ=255)
                                                                                                   cgcgACGTTGTTCATCACTTAACACTTCCTAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATcg                                   >  2:532230/1‑140 (MQ=255)
                                                                                                   cgcgACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATcg                                   >  1:413805/1‑140 (MQ=255)
                                                                                                   cgcgACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATcg                                   >  1:73079/1‑140 (MQ=255)
                                                                                                   cgcgACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATcg                                   >  2:37947/1‑140 (MQ=255)
                                                                                                    gcgACGTTGTTCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATcgt                                  <  1:369386/140‑1 (MQ=255)
                                                                                                         gttgttCATCACTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATcgtccgtc                             <  1:493010/140‑1 (MQ=255)
                                                                                                               catcaCTTAACACTTCCCAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTc                       <  1:189273/140‑1 (MQ=255)
                                                                                                                              cccAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCg        >  1:155973/1‑140 (MQ=255)
                                                                                                                              cccAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCg        >  1:71515/1‑140 (MQ=255)
                                                                                                                              cccAGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCg        >  1:497696/1‑140 (MQ=255)
                                                                                                                                 aGGCTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCGtaa     >  1:207750/1‑140 (MQ=255)
                                                                                                                                    cTTCAACGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCGtaataa  <  1:466353/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAA  >  NC_000913/1063490‑1063761

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