Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,638 G→A 7.7% D94N (GAC→AAC)  alaC → aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,6380GA 7.7% 149.6 / 2.7 52D94N (GAC→AAC) alaCaminotransferase
Reads supporting (aligned to +/- strand):  ref base G (28/20);  new base A (2/2);  total (30/22)
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.55e-01

ACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGATTCACATCTACGGCGCAGTGATCGC  >  NC_002947/955492‑955783
                                                                                                                                                   |                                                                                                                                                   
aCCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGa                                                                                                                                                    >  6:178067/1‑148 (MQ=255)
  cccGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACAt                                                                                                                                                 >  4:9700/1‑149 (MQ=255)
  cccGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACAt                                                                                                                                                 <  5:178067/149‑1 (MQ=255)
     gATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGAc                                                                                                                                                   >  1:155418/1‑144 (MQ=255)
     gATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGAc                                                                                                                                                   <  2:155418/144‑1 (MQ=255)
         gcgcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGAccc                                                                                                                                           <  6:29326/148‑1 (MQ=255)
           gcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCgg                                                                                                                                         >  3:129693/1‑148 (MQ=255)
                 ccgccgCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCg                                                                                                                                   >  6:57169/1‑148 (MQ=255)
                  cgccgcACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGaa                                                                                                                                 >  4:87628/1‑149 (MQ=255)
                                 agaAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATc                                                                                                                  <  2:6787/149‑1 (MQ=255)
                                  gaAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑g                                                                                                                <  5:57169/149‑1 (MQ=255)
                                     gCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCt                                                                                                            >  7:208949/1‑149 (MQ=255)
                                                tCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGc                                                                                                 <  3:9700/149‑1 (MQ=255)
                                                    ccAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑C‑GCTCGAAAGAGGGCCTgg                                                                                             >  5:31515/1‑149 (MQ=255)
                                                           gAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCt                                                                                      >  7:32813/1‑149 (MQ=255)
                                                           gAAGACACCCACGGCTACTCCACGTCGCACGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCt                                                                                      >  1:157865/1‑149 (MQ=255)
                                                            aaGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCt                                                                                      >  8:222799/1‑148 (MQ=255)
                                                             aGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCACAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCt                                                                                      >  1:73505/1‑148 (MQ=255)
                                                                 acCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCAACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCt                                                                                >  7:253297/1‑149 (MQ=255)
                                                                 acCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCAACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCt                                                                                >  1:45248/1‑149 (MQ=255)
                                                                  cccACGGCTACTGGACGTCGCGCGGCATCCCGCGGCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTg                                                                               <  6:31515/149‑1 (MQ=255)
                                                                  cccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAagag                                                                                                    >  6:153898/1‑128 (MQ=255)
                                                                  cccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAagag                                                                                                    <  5:153898/128‑1 (MQ=255)
                                                                  cccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTg                                                                               <  3:87628/149‑1 (MQ=255)
                                                                   ccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTgg                                                                              <  1:204134/149‑1 (MQ=255)
                                                                   ccACGGCTACTCCACGTCGCGCGGCATCCCGAGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGACAGCGAAGCCATCGTCACCATC‑G‑CCTCGAAAGAGAGCCTGGACCACATCATGCTgg                                                                              >  5:97915/1‑149 (MQ=255)
                                                                       ggCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCAc                                                                          >  2:59645/1‑149 (MQ=255)
                                                                        gCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCAcc                                                                         <  7:222799/149‑1 (MQ=255)
                                                                          tACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑GCGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCa                                                                           <  1:180053/146‑1 (MQ=255)
                                                                          tACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCa                                                                           >  2:180053/1‑145 (MQ=255)
                                                                               cACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTAACAATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTAGAc                                                                   >  7:167354/1‑148 (MQ=255)
                                                                                   gcgcgcgGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGG‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAgg                                                               <  1:59645/148‑1 (MQ=255)
                                                                                                     cTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGcc                                            >  4:112964/1‑149 (MQ=255)
                                                                                                                    tc‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCa                                    <  3:265357/142‑1 (MQ=255)
                                                                                                                    tc‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCa                                    >  4:265357/1‑142 (MQ=255)
                                                                                                                    tc‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTAccc                             >  5:112630/1‑149 (MQ=255)
                                                                                                                       cAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCt                                                  <  3:165720/126‑1 (MQ=255)
                                                                                                                       cAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCt                                                  >  4:165720/1‑126 (MQ=255)
                                                                                                                       cAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGa                           <  8:208949/149‑1 (MQ=255)
                                                                                                                        aaCTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGAt                          >  8:92964/1‑149 (MQ=255)
                                                                                                                         aCTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGAtt                         <  8:32813/149‑1 (MQ=255)
                                                                                                                          cTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTgg                                                                              <  7:117168/95‑1 (MQ=255)
                                                                                                                          cTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTgg                                                                              >  8:117168/1‑95 (MQ=255)
                                                                                                                          cTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGATTc                        >  8:41140/1‑149 (MQ=255)
                                                                                                                               aCAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGAc                                                                   <  3:53865/101‑1 (MQ=255)
                                                                                                                               aCAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGAc                                                                   >  4:53865/1‑101 (MQ=255)
                                                                                                                                        gCTACGAGGTCAACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGATTCACATCTACGGCGCa          <  2:45248/149‑1 (MQ=255)
                                                                                                                                        gCTACGAGGTCAACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGATTCACATCTACGGCGCa          <  8:253297/149‑1 (MQ=255)
                                                                                                                                         cTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTAGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCAACCCGATTCACATCTACGGCGCAg         >  7:60972/1‑149 (MQ=255)
                                                                                                                                            cGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATAc                                                        >  6:22269/1‑99 (MQ=255)
                                                                                                                                            cGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATAc                                                        <  5:22269/99‑1 (MQ=255)
                                                                                                                                                gTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGATTCACATCTACGGCGCAGTGATCGc  >  2:40152/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTC‑CAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATC‑G‑GCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGCGATACGGTGCTGGTGCCCAACCCCAGCTACCCGATTCACATCTACGGCGCAGTGATCGC  >  NC_002947/955492‑955783

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