Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,419,720 T→C 37.5% D183G (GAT→GGT)  hisJ ← histidine/lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,419,7200TC37.5% 35.1 / 45.7 53D183G (GAT→GGT) hisJhistidine/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (13/20);  new base C (14/6);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 4.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCG  >  CP009273/2419577‑2419865
                                                                                                                                               |                                                                                                                                                 
cTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTgg                                                                                                                                              >  1:422135/1‑149 (MQ=255)
  tCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGaa                                                                                                                                            <  1:75063/149‑1 (MQ=255)
    gcgcAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAAcg                                                                                                                                          <  2:360782/149‑1 (MQ=255)
     cgcAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAAcgc                                                                                                                                         >  1:409400/1‑149 (MQ=255)
      gcAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAAcgcg                                                                                                                                        <  1:450832/149‑1 (MQ=255)
           tCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAATCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATc                                                                                                                                   <  2:377886/149‑1 (MQ=255)
              ttATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAAt                                                                                                                                >  1:393933/1‑149 (MQ=255)
                      tttACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCCTCACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCgg                                                                                                                        <  2:290676/149‑1 (MQ=255)
                      tttACGCAGGCCCACGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCgg                                                                                                                        <  2:393647/149‑1 (MQ=255)
                      ctgccGCAGCCGCATGGCGGTCTCGATGCCCACCAGTTTTCACTGTTTAACATAGGGGCGACCGGAGTTGTCGTTTTCACCGACTGGTTGTTGGATGAAAGCTTCGCAGGTAGAGACCTCATCCTGGAACGCGGCATCAATACGGCCgg                                                                                                                        <  2:290809/145‑1 (MQ=255)
                          cGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAgt                                                                                                                    <  2:190787/149‑1 (MQ=255)
                                 ccATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAgtcaggtca                                                                                                             >  2:397221/1‑149 (MQ=255)
                                        ggaccccACGCCAAACCGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATaaa                                                                                                      <  1:392364/146‑1 (MQ=255)
                                         gTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTGGGTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAAt                                                                                                     >  2:444072/1‑149 (MQ=255)
                                              tACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCCGGTCAGAATAAAtgttgt                                                                                                >  1:205563/1‑149 (MQ=255)
                                              tACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAAtgttgt                                                                                                <  2:422135/149‑1 (MQ=255)
                                               aCGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTc                                                                                               <  1:230743/149‑1 (MQ=255)
                                                  ccAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGCCAGAGTAAATGTTGTCCTg                                                                                            >  1:410475/1‑149 (MQ=255)
                                                        aGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTg                                                                                      <  2:225256/149‑1 (MQ=255)
                                                              tcctCTTTAACAGCCGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACga                                                                                <  1:160437/146‑1 (MQ=255)
                                                                aTCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACgaga                                                                              >  1:104176/1‑149 (MQ=255)
                                                                 tCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGACACCTTCGCTGGCAGCGACCTCACCCTGGAACGCCGCATCCATCCGTCCGGCAGTCAGGTCCGAATAAAGGTGGGCCTGCCCCTGATACGAGAc                                                                             <  2:274554/149‑1 (MQ=255)
                                                                    ttAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGAt                                                                          >  2:196558/1‑149 (MQ=255)
                                                                    ttAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGAt                                                                          >  1:84236/1‑149 (MQ=255)
                                                                     tAACAGACGGGCCACCGACTTTGGAATCTTTGCCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGAtt                                                                         <  2:223142/149‑1 (MQ=255)
                                                                      aaCAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGAttt                                                                        <  2:393933/149‑1 (MQ=255)
                                                                            cGGGCCACCGAATTTGTAATCTTTCCCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGc                                                                  <  1:61758/149‑1 (MQ=255)
                                                                                ccACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCttt                                                              <  2:389083/149‑1 (MQ=255)
                                                                                  aCCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTg                                                            <  1:167805/149‑1 (MQ=255)
                                                                                           gTAATCTTTACCGACAGGTTGTTTGAGGAAACCTGAGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTCGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAAt                                                   <  1:261467/149‑1 (MQ=255)
                                                                                            tAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAAGGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATg                                                  >  2:259414/1‑149 (MQ=255)
                                                                                            tAATCTTTACCGACAGGTTCTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCGGGATCGAGGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATg                                                  >  1:340783/1‑149 (MQ=255)
                                                                                               tCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTc                                               >  2:294812/1‑149 (MQ=255)
                                                                                               tCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTc                                               >  2:387383/1‑149 (MQ=255)
                                                                                                 tttACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCAt                                             <  1:360780/149‑1 (MQ=255)
                                                                                                    aCCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTAc                                          <  1:355466/149‑1 (MQ=255)
                                                                                                      cGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGCCGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCg                                        >  2:196631/1‑149 (MQ=255)
                                                                                                               gttTGAGGAACCCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTcc                               <  2:55908/149‑1 (MQ=255)
                                                                                                                tttGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCt                              >  1:120240/1‑149 (MQ=255)
                                                                                                                tttGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCt                              >  1:77980/1‑149 (MQ=255)
                                                                                                                tttGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCt                              >  1:167375/1‑149 (MQ=255)
                                                                                                                   gAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTggg                           >  1:328723/1‑149 (MQ=255)
                                                                                                                     ggaagctctCGCTGGCAGCAGCCTCCCCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGCTACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGTACGTCTCCTGggtg                         <  2:345621/141‑1 (MQ=255)
                                                                                                                      gAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGggtgg                        >  1:274656/1‑149 (MQ=255)
                                                                                                                       aaaCCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGggtggt                       >  1:433219/1‑149 (MQ=255)
                                                                                                                           cTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGccc                   >  2:286200/1‑149 (MQ=255)
                                                                                                                                   gCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAAtgc           >  1:70926/1‑147 (MQ=255)
                                                                                                                                     aGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCTTACTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGc         >  1:340735/1‑149 (MQ=255)
                                                                                                                                     aGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGc         <  1:140526/149‑1 (MQ=255)
                                                                                                                                        gACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTTCAATACGCCTa      >  2:348248/1‑149 (MQ=255)
                                                                                                                                         aCCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGTTGCCCAATGTTCATTACCGAACGTCTCCTGGGGGGTGCCCTGCAATGCGccctc     >  2:421474/1‑146 (MQ=255)
                                                                                                                                          ccTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTAcc    <  1:1047/149‑1 (MQ=255)
                                                                                                                                            tCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCg  >  1:269615/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCG  >  CP009273/2419577‑2419865

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