Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,688,221 C→T 10.3% L99L (CTG→TTG)  PP_3253 → carboxylate‑amine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,688,2210CT10.3% 96.7 / 3.9 39L99L (CTG→TTG) PP_3253carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base T (2/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

AGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAC  >  NC_002947/3688078‑3688361
                                                                                                                                               |                                                                                                                                            
aGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCtgctgc                                                                                                                                        <  2:297385/150‑1 (MQ=255)
  aTCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAg                                                                                                                                      >  2:603720/1‑150 (MQ=255)
    cGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                    <  1:603720/150‑1 (MQ=255)
          ggCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCATCCGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGGGCTTCTGCAGCAGccacgc                                                                                                                             <  1:685766/150‑5 (MQ=255)
                          aCCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGTCTTCTGCAGAAGCCCGCGGCACAAGGGCACtat                                                                                                              <  2:690238/150‑1 (MQ=255)
                          aCCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACtat                                                                                                              <  1:285523/150‑1 (MQ=255)
                           ccAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACtata                                                                                                             >  1:563860/1‑150 (MQ=255)
                                cTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCATGCGCCGCGTGGCTGCTGCAGACGCCCGCGGCACAAGGGCACTATAagcag                                                                                                        >  2:283491/1‑150 (MQ=255)
                                                      tcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAg                                                                                                                      >  1:417957/1‑114 (MQ=255)
                                                      tcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAg                                                                                                                      <  2:417957/114‑1 (MQ=255)
                                                       ccAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCAtgt                                                                                 >  2:297489/1‑150 (MQ=255)
                                                             cAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                           <  2:563860/150‑1 (MQ=255)
                                                              aaCCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCa                                                                                                                  <  1:460072/110‑1 (MQ=255)
                                                              aaCCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCa                                                                                                                  >  2:460072/1‑110 (MQ=255)
                                                              aaCCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAgg                                                                          >  2:341729/1‑150 (MQ=255)
                                                               actgTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGc                                                                         <  2:293870/147‑1 (MQ=255)
                                                                        ggCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGc                                                                <  1:341729/150‑1 (MQ=255)
                                                                         gCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAgg                                                                                                                     >  2:391639/1‑96 (MQ=255)
                                                                         gCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAgg                                                                                                                     <  1:391639/96‑1 (MQ=255)
                                                                                 tGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCTCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGcc                                                       >  1:258454/1‑150 (MQ=255)
                                                                                        ggCCGACAAGGGAATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAtgt                                                <  2:258454/150‑1 (MQ=255)
                                                                                            gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGTTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAagcag                                                                                                        <  2:652439/90‑1 (MQ=255)
                                                                                            gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGTTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAagcag                                                                                                        <  2:256213/90‑1 (MQ=255)
                                                                                            gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGTTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAagcag                                                                                                        >  1:652439/1‑90 (MQ=255)
                                                                                            gAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGTTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAagcag                                                                                                        >  1:256213/1‑90 (MQ=255)
                                                                                                    gATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTg                                        <  2:421104/146‑1 (MQ=255)
                                                                                                    gATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTg                                        >  1:421104/1‑146 (MQ=255)
                                                                                                         ggCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTATCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgc                               >  1:406208/1‑150 (MQ=255)
                                                                                                             aTACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGc                           >  1:319084/1‑150 (MQ=255)
                                                                                                               aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTc                                                                                                  >  2:473706/1‑77 (MQ=255)
                                                                                                               aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTc                                                                                                  <  1:473706/77‑1 (MQ=255)
                                                                                                                 gCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgt                       <  1:23568/150‑1 (MQ=255)
                                                                                                                   tGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                     <  2:319084/150‑1 (MQ=255)
                                                                                                                            ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGc                                                        <  2:450657/106‑1 (MQ=255)
                                                                                                                            ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGc                                                        >  1:450657/1‑106 (MQ=255)
                                                                                                                                cccTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGGAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTg        >  1:299856/1‑150 (MQ=255)
                                                                                                                                  cTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTg                                   >  1:573815/1‑121 (MQ=255)
                                                                                                                                  cTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTg                                   <  2:573815/121‑1 (MQ=255)
                                                                                                                                      gccgcGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATTGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGTAACTGATCAAc  <  1:297489/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
AGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCCCGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAC  >  NC_002947/3688078‑3688361

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