Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,650,833 A→C 11.6% S5A (TCA→GCA)  gcvPB ← glycine dehydrogenase (decarboxylating) subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,650,8330AC11.6% 84.0 / 5.2 42S5A (TCA→GCA) gcvPBglycine dehydrogenase (decarboxylating) subunit 2
Reads supporting (aligned to +/- strand):  ref base A (14/23);  new base C (2/3);  total (16/26)
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.73e-01

AACTTCAGGAAACTCCGCTTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACCTAAATCAAAACCACCGATAATAT  >  CP000730/1650693‑1650971
                                                                                                                                            |                                                                                                                                          
aaCTTCAGGAAACTCCGCTTTGTTTTTACGAATAAATTTATCATCTAAACATGACTTATCAGAATTTGTGAGAATCTCACTTTTTGGTAATGAGTATTCATATCTGCCTTCACTAGATCTTTCTAAAATTAATGGACTTGa                                                                                                                                            >  1:80621/1‑141 (MQ=255)
  ctgcAGGAAACTCCGATTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTTATGGACTTGAtt                                                                                                                                          <  1:243386/138‑1 (MQ=255)
  cTTCAGGAAACTCCGCTTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGAtt                                                                                                                                          <  1:34527/141‑1 (MQ=255)
  cTTCAGGAAACTCCGCTTTGTTTTTACGAATAAATTTATAATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGAtt                                                                                                                                          >  1:66874/1‑141 (MQ=255)
    taaGGAAACTCCGCTTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTa                                                                                                                                        <  2:183912/139‑1 (MQ=255)
      aGGAAACTCCGCTTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACt                                                                                                                                      >  1:101458/1‑141 (MQ=255)
                gCTTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTAATTTCCTAGTCATTTaa                                                                                                                            >  2:97081/1‑141 (MQ=255)
                             gAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCttt                                                                                                               <  2:24179/141‑1 (MQ=255)
                                      tATCCTCTAACTATGACTCAACAGAATTTGTTTTAATATCACTTTTTAGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATGTAACTCACCGGCCTTTTCTACAAAt                                                                                                      <  1:172050/141‑1 (MQ=255)
                                       atcatcTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTaataatt                                                                                                     >  2:267950/1‑135 (MQ=255)
                                       atcatcTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATg                                                                                                     <  1:272566/141‑1 (MQ=255)
                                            cTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTGTCTACTAATGTATCg                                                                                                >  1:109626/1‑141 (MQ=255)
                                                 aaTTACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTc                                                                                           <  2:135724/141‑1 (MQ=255)
                                                    taCTCAATAGAATTTGTTTTAATATCACTTTTTTGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATc                                                                                        <  2:169102/140‑1 (MQ=255)
                                                     aCTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTCCCAATGTCTCTATTTCATCt                                                                                       >  1:17811/1‑141 (MQ=255)
                                                       tCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGCTTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCttt                                                                                     >  1:201806/1‑141 (MQ=255)
                                                       tCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGCTTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCttt                                                                                     >  1:201673/1‑141 (MQ=255)
                                                              aaTTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCtt                                                                              <  2:80621/141‑1 (MQ=255)
                                                                      ttAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGt                                                                      <  1:35232/141‑1 (MQ=255)
                                                                            tCACTTTTTTGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGc                                                                <  2:101458/141‑1 (MQ=255)
                                                                              gcTTGTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTa                                                              <  1:160378/140‑1 (MQ=255)
                                                                              aCTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTa                                                              >  2:17746/1‑141 (MQ=255)
                                                                              aCTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTa                                                              >  2:153311/1‑141 (MQ=255)
                                                                                  tttttttttttttatGCATATCTACCTTCACAAGATCTTTCAAAAATTAATGGAGTTGATTTACTAGTCATTTAACTCACCAACCTTTTCTTTTTTTTTTTTGTTTTCTTCTTTTGTTCTTATTTCATTTACAGCTATTaa                                                          >  2:229033/13‑141 (MQ=255)
                                                                                     ggTAATGAATATGCATATCTTCCTTCACGAAATCTTTCAAAAATTAATTTACTATATTTACTAATCATTTAACTCACCAGCCTTTTCTACAAATTTATCTATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAt                                                       >  2:86497/1‑141 (MQ=255)
                                                                                     ggTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAt                                                       <  2:59118/141‑1 (MQ=255)
                                                                                     ggTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAt                                                       <  1:327492/141‑1 (MQ=255)
                                                                                       tAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAtgt                                                     <  2:170519/141‑1 (MQ=255)
                                                                                        aaTGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAtgtg                                                    <  1:276749/141‑1 (MQ=255)
                                                                                        aaTGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAtgtg                                                    <  1:142198/141‑1 (MQ=255)
                                                                                        aaTGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACAtgtg                                                    <  2:163252/141‑1 (MQ=255)
                                                                                         aTGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGa                                                   >  1:333585/1‑141 (MQ=255)
                                                                                                tGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTaa                                            <  1:196367/141‑1 (MQ=255)
                                                                                                     atCTTCCTGCACGAGAGCTCTAAAAAACTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAgtcg                                       <  1:57192/141‑1 (MQ=255)
                                                                                                      tCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGCTTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAgtcgt                                      <  2:201806/141‑1 (MQ=255)
                                                                                                      tCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGCTTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATCCAGTTACAGCTATTAACATGTGATTTTTAAAgtcgt                                      <  2:201673/141‑1 (MQ=255)
                                                                                                                   gaTCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATTGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAATTCGTCTGAAACTAGAcc                         >  2:278326/1‑141 (MQ=255)
                                                                                                                           aaaaaTTAATGGACTTGCTTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTACTTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACCTAAATCaa                 <  2:86975/141‑1 (MQ=255)
                                                                                                                                   aTGGACGTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACCTAAATCAAAACCACCg         <  2:74921/141‑1 (MQ=255)
                                                                                                                                    tgcaCTTGATTTACTAGTCATTTAACCCACCAGCCTTCTCTACAAATGTATCGATCTCATCCATTGTTATCAATTCAGTTACACCTATCAACATGCGCTTTTACAAGTCGCCTCACACACCACCTAACTCAAAacaccaga        >  2:34303/4‑135 (MQ=255)
                                                                                                                                     ggACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACATAAATCAAAACCACCGaa       >  2:150180/1‑140 (MQ=255)
                                                                                                                                          tGATTTCCTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACCTAAATCAAAACCACCGATAatat  <  1:278698/141‑1 (MQ=255)
                                                                                                                                          tGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACCTAAATCAAAACCACCGATAatat  <  1:97081/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AACTTCAGGAAACTCCGCTTTGTTTTTACGAATAAATTTATCATCTAACAATGACTCAACAGAATTTGTTTTAATATCACTTTTTGGTAATGAATATGCATATCTGCCTTCACGAGATCTTTCAAAAATTAATGGACTTGATTTACTAGTCATTTAACTCACCAGCCTTTTCTACAAATGTATCGATTTCATCTTTTGTTCTTAATTCAGTTACAGCTATTAACATGTGATTTTTAAAGTCGTCTGAAACAACACCTAAATCAAAACCACCGATAATAT  >  CP000730/1650693‑1650971

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