Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,731,726 T→C 8.5% N361D (AAT→GAT)  sucB ← 2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,731,7260TC 8.5% 131.3 / 2.8 47N361D (AAT→GAT) sucB2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (2/2);  total (24/23)
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.19e-01

AGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCTTCGATGGACAGTTTGC  >  NC_002947/4731577‑4731875
                                                                                                                                                     |                                                                                                                                                     
agagGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATAACGACCTGGCCAtt                                                                                                                                                       >  1:189823/1‑150 (MQ=255)
   gggTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGAt                                                                                                                                                    <  2:454425/150‑1 (MQ=255)
   gggTCAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGAt                                                                                                                                                    <  1:1073842/150‑1 (MQ=255)
    ggTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATg                                                                                                                                                   >  2:806508/1‑150 (MQ=255)
    ggTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTCACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATg                                                                                                                                                   <  1:806508/150‑1 (MQ=255)
            gTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATCAATGGCCaggtc                                                                                                                                           <  1:1055154/150‑5 (MQ=255)
                        gCGGGACGGCTCTCCCAGCAGGCTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCAGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCGGgctgag                                                                                                                               >  2:198941/1‑149 (MQ=255)
                                 aTCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCAt                                                                                                                      <  2:71741/150‑1 (MQ=255)
                                  ttttCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCAATCGGGCGCTGGATGCagttgtctatt                                                                                                                     >  2:1086038/3‑145 (MQ=255)
                                  tCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATg                                                                                                                     >  1:964428/1‑150 (MQ=255)
                                    ttCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGCAAGGCTGCCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGcc                                                                                                                   >  1:345921/1‑150 (MQ=255)
                                          cagGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAAt                                                                                                             <  1:383812/150‑1 (MQ=255)
                                              ttCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATCAATGGCCaggtc                                                                                                                                           <  1:1126526/116‑5 (MQ=255)
                                              ttCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATCAATGGCCaggtc                                                                                                                                           >  2:1126526/1‑112 (MQ=255)
                                                      aGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTgcgg                                                                                                 <  2:824353/150‑1 (MQ=255)
                                                         gACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGgatga                                                                                                                                <  1:191282/116‑1 (MQ=255)
                                                         gACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGgatga                                                                                                                                >  2:191282/1‑116 (MQ=255)
                                                         gACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGcggcgg                                                                                              <  1:729074/150‑1 (MQ=255)
                                                          aCCAGGAAGGTTACCGCTTCCTTGTCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCggg                                                                                             <  2:145676/150‑1 (MQ=255)
                                                          aCCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCggg                                                                                             >  1:493207/1‑150 (MQ=255)
                                                           ccAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGt                                                                                            <  2:658243/150‑1 (MQ=255)
                                                              ggAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGa                                                                                         <  2:117916/150‑1 (MQ=255)
                                                                      aCCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCggg                                                                                 <  2:964079/150‑1 (MQ=255)
                                                                      aCCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCCATGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCggg                                                                                 <  1:198941/150‑1 (MQ=255)
                                                                        cGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGt                                                                               <  2:92465/150‑1 (MQ=255)
                                                                             ccTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGAca                                                                          >  1:693566/1‑150 (MQ=255)
                                                                              cTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGAcat                                                                         >  2:1094435/1‑150 (MQ=255)
                                                                                        aTCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATg                                                                                                            <  1:266369/105‑1 (MQ=255)
                                                                                        aTCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATg                                                                                                            >  2:266369/1‑105 (MQ=255)
                                                                                        aTCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCg                                                               >  1:530216/1‑150 (MQ=255)
                                                                                         tCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGa                                                              >  1:1054982/1‑150 (MQ=255)
                                                                                          cAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGaa                                                             <  2:693566/150‑1 (MQ=255)
                                                                                                  gATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCAc                                                     >  2:713954/1‑150 (MQ=255)
                                                                                                    tCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCg                                                   >  1:1107501/1‑150 (MQ=255)
                                                                                                        aCGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTtgg                                               >  1:432746/1‑150 (MQ=255)
                                                                                                          gACNGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTtggtg                                             >  2:701814/1‑150 (MQ=255)
                                                                                                          gACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTtggtg                                             >  1:221813/1‑150 (MQ=255)
                                                                                                                      tACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTg                                 >  2:726674/1‑150 (MQ=255)
                                                                                                                          tcatcGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCAc                             >  1:590382/1‑150 (MQ=255)
                                                                                                                            atcGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCg                           >  2:114390/1‑150 (MQ=255)
                                                                                                                                 gCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCAt                      <  2:380099/150‑1 (MQ=255)
                                                                                                                                    aaTCACGACCGGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGAGCGGGGTCGACATCATCGCACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCAtctc                   <  2:345921/150‑1 (MQ=255)
                                                                                                                                      tCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCtt                 <  2:189823/150‑1 (MQ=255)
                                                                                                                                         cGACCTGGCCATTGATGGTCATCGGGCGCTGTATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCTTCGa              <  2:141013/150‑1 (MQ=255)
                                                                                                                                           cgatGGCCATCAATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCTTCGATg            >  1:119745/4‑150 (MQ=255)
                                                                                                                                             cTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGAcatcat                                                                      >  2:793023/1‑90 (MQ=255)
                                                                                                                                             cTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGAcatcat                                                                      <  1:793023/90‑1 (MQ=255)
                                                                                                                                               ggCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCAACATCATCGAACCGAAGGGGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCTTCGATGGACa        >  1:860262/1‑150 (MQ=255)
                                                                                                                                                     tGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCTTCGATGGACAGTTTGc  <  1:215047/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAGGTACATCATCGGGCGAATCACGACCTGGCCATTGATGGCCATCGGGCGCTGGATGATGTTGTGCATGCCGAGAATGGCGGCCTGCGGCGGGTTGACGATCGGGGTCGACATCATCGAACCGAAGGTGCCACCGTTGGTGATGGTGAAGGTGCCACCGGTCATCTCTTCGATGGACAGTTTGC  >  NC_002947/4731577‑4731875

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: ‑