Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 134,393 G→A L12L (TTG→TTA sbnA → possible pyridoxal‑phosphate dependent enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730134,3930GA96.7% 102.7 / NA 30L12L (TTG→TTAsbnApossible pyridoxal‑phosphate dependent enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/19);  minor base C (1/0);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTGTTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCAAGTAC  >  CP000730/134257‑134516
                                                                                                                                        |                                                                                                                           
aaTGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAg                                                                                                                         <  1:307005/141‑1 (MQ=255)
aaTGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAg                                                                                                                         <  1:379630/141‑1 (MQ=255)
  tGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGAt                                                                                                                       <  1:346155/141‑1 (MQ=255)
    cATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTc                                                                                                                     <  1:403465/141‑1 (MQ=255)
      tAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTAATTATTAGATTCTg                                                                                                                   >  1:305926/1‑141 (MQ=255)
                     tAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTa                                                                                                    <  1:126412/141‑1 (MQ=255)
                               aTAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACt                                                                                          <  1:127869/141‑1 (MQ=255)
                                    gagaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatc                                                                                     <  1:228557/141‑1 (MQ=255)
                                    gagaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatc                                                                                     <  1:431808/141‑1 (MQ=255)
                                     agaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatca                                                                                    <  1:410440/141‑1 (MQ=255)
                                     agaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatca                                                                                    <  1:380746/141‑1 (MQ=255)
                                       aaTCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAAc                                                                                  <  1:96343/141‑1 (MQ=255)
                                                  aTTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGaaa                                                                       >  1:350984/1‑141 (MQ=255)
                                                     aGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACAt                                                                    <  1:346424/141‑1 (MQ=255)
                                                        tAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATGGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGaa                                                                 >  1:253380/1‑141 (MQ=255)
                                                                   attatAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGAGTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCaaa                                                      >  1:206765/1‑141 (MQ=255)
                                                                          acttttAAATTTGAGGAGGAAGCGCTGTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTagag                                               <  1:172287/139‑1 (MQ=255)
                                                                               tAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTCTTAGATTCTGTAGGGCAAACACCTATGGTGTTACGTCATCAACTATTTCCGAAACATGAAGGGGTTGCAAAGTTAGAGTatat                                          >  1:385519/1‑141 (MQ=255)
                                                                                aaaTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATg                                         >  1:374317/1‑141 (MQ=255)
                                                                                aaaTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATg                                         <  1:389190/141‑1 (MQ=255)
                                                                                aaaTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCACCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATg                                         <  1:140154/141‑1 (MQ=255)
                                                                                  aTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCGAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGaa                                       >  1:315420/1‑141 (MQ=255)
                                                                                               gcgcTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGgcagc                          >  1:157018/1‑141 (MQ=255)
                                                                                               gcgcTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACAACTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGgcagc                          >  1:439828/1‑141 (MQ=255)
                                                                                                  cTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATg                       <  1:304005/141‑1 (MQ=255)
                                                                                                       gattgaAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGa                  <  1:256246/141‑1 (MQ=255)
                                                                                                         ttgaAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAGCTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGGAAGATc                >  1:84566/1‑141 (MQ=255)
                                                                                                         ttgaAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATc                >  1:171650/1‑141 (MQ=255)
                                                                                                                  gTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCa       <  1:480715/141‑1 (MQ=255)
                                                                                                                       gCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTAAAATACCCATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCAAGTAc  <  1:488726/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                           
AATGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTGTTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCAAGTAC  >  CP000730/134257‑134516

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