Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 217,036 T→A 15.4% intergenic (‑203/+23) argD ← / ← brnQ ornithine‑‑oxo‑acid transaminase/possible LIVCS family branched chain amino acid:cation symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730217,0360TA15.4% 94.6 / 10.6 38intergenic (‑203/+23)argD/brnQornithine‑‑oxo‑acid transaminase/possible LIVCS family branched chain amino acid:cation symporter
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base A (3/3);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAGATAAGGTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAATGTCTGATAATGGTAAATATT  >  CP000730/216899‑217171
                                                                                                                                         |                                                                                                                                       
ataAAGATAAGGTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAAc                                                                                                                                      >  2:193169/1‑141 (MQ=255)
     gATAAGGTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATg                                                                                                                                 <  1:135323/141‑1 (MQ=255)
     gATAAGGTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATg                                                                                                                                 <  1:135331/141‑1 (MQ=255)
       tagggTTGTGTGTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTAATATGTATGTATACGATTGCAGTTTAATGTAAAAAAAATAAAAAATGCCACAAGCATATGTTTAACATATGAc                                                                                                                               <  2:314857/138‑1 (MQ=255)
         aGGTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgcttg                                                                                                                             <  2:140475/141‑5 (MQ=255)
           gTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAGAATGCCATAAGCATATGTTTAACATATGACtgtg                                                                                                                           >  2:254948/1‑141 (MQ=255)
            ttGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTgg                                                                                                                          <  1:145601/141‑1 (MQ=255)
              gtgtTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCa                                                                                                                        >  2:120656/1‑141 (MQ=255)
                    gTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtat                                                                                                                  <  1:89524/141‑1 (MQ=255)
                    gTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtat                                                                                                                  >  1:20140/1‑141 (MQ=255)
                          aTGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCt                                                                                                            >  2:149259/1‑141 (MQ=255)
                                  tGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAATATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAt                                                                                                    <  2:327899/141‑1 (MQ=255)
                                          gAAAGTACTAATATTTTTGGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTATTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTGGTGTAATTTTACGAt                                                                                            >  2:78457/1‑141 (MQ=255)
                                            aaGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTTTGAATAAAAAATGCCACAAGCCTATGTTTAACATATGACTGTGGAATAGTATTATGCTTGTGGAATTTTACGATGc                                                                                          >  1:248874/1‑141 (MQ=255)
                                            aaGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGc                                                                                          <  2:120413/141‑1 (MQ=255)
                                             aGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                                                                         <  2:226037/141‑1 (MQ=255)
                                              gTACTAATATTTTTTTTTTTTTTTTATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTTTTTTTTTTTAATTTTACGATGCtt                                                                                        >  2:106740/1‑141 (MQ=255)
                                                 ttAATATTTTTTGTGGTTAATCATTGATAGGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAt                                                                                     <  2:294890/140‑1 (MQ=255)
                                                            tgtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataaataa                                                                          <  1:120656/141‑1 (MQ=255)
                                                                  tAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGc                                                                    <  2:95356/141‑1 (MQ=255)
                                                                     tCATTGATATGTATGGAGACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGa                                                                 >  2:150426/1‑141 (MQ=255)
                                                                                         gATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTg                                             <  1:283763/141‑1 (MQ=255)
                                                                                            tGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTGTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTgaga                                          >  1:190352/1‑141 (MQ=255)
                                                                                             gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaa                                         >  1:246909/1‑141 (MQ=255)
                                                                                                  ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAAcc                                    >  1:184456/1‑141 (MQ=255)
                                                                                                   tAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAAccc                                   >  1:92068/1‑141 (MQ=255)
                                                                                                       gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGcc                               <  1:149259/141‑1 (MQ=255)
                                                                                                          gtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCAt                            >  1:92361/1‑141 (MQ=255)
                                                                                                                aTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAAt                      <  2:288254/141‑1 (MQ=255)
                                                                                                                aTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAAt                      <  1:226835/141‑1 (MQ=255)
                                                                                                                    aaaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCAt                            <  2:7461/131‑1 (MQ=255)
                                                                                                                    aaaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCAt                            >  1:7461/1‑131 (MQ=255)
                                                                                                                       aTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAATGTCTGAt               >  2:91345/1‑141 (MQ=255)
                                                                                                                               gtcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAATGTCTGATAATGGTaa       >  2:136297/3‑141 (MQ=255)
                                                                                                                                gCATCTGTTTACCATATGAATGTGGCATAGTATTACGCTTGAGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGCAACAAATCATACTGATAAAAACCCTGCCCATACCAATGTCTGAAAATGGTaaa      <  2:80756/141‑1 (MQ=255)
                                                                                                                                 cATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCTATGTCCGATAATGGCaaag     <  2:243137/141‑2 (MQ=255)
                                                                                                                                 cATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAATGTCTGATAATGGTAAat     >  2:30441/1‑141 (MQ=255)
                                                                                                                                    atGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAATGTCTGATAATGGTAAATAtt  >  2:237686/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                       
ATAAAGATAAGGTTGTGTTTGTATTAATGTGATTTGTAAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAACCCAGCCCATACCAATGTCTGATAATGGTAAATATT  >  CP000730/216899‑217171

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