Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,073,800 C→A 10.5% S136Y (TCT→TAT)  purE → 5‑(carboxyamino)imidazole ribonucleotide mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,073,8000CA10.5% 169.2 / 4.5 53S136Y (TCT→TAT) purE5‑(carboxyamino)imidazole ribonucleotide mutase
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base A (2/4);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAAGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAGATGATGGC  >  CP000730/1073664‑1073936
                                                                                                                                        |                                                                                                                                        
ttAAAGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTa                                                                                                                                      >  1:685830/1‑141 (MQ=255)
 tAAAGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAg                                                                                                                                     >  1:117337/1‑141 (MQ=255)
  aaaGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGt                                                                                                                                    <  1:513009/141‑1 (MQ=255)
  aaaGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGt                                                                                                                                    >  2:134658/1‑141 (MQ=255)
  aaaGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTACTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGt                                                                                                                                    >  1:651399/1‑141 (MQ=255)
          taGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCTTGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGtataatgaa                                                                                                                            >  1:695952/1‑133 (MQ=255)
            aaTTCTTTAAAATCAATTGTACAAATGCCAGGAGGTATACCGGTTGCAACGACTACAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATATTAAGTATTCAAAATACAAAAATAGATGAAAAACTa                                                                                                                          <  2:272904/140‑1 (MQ=255)
              ttCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGtttttcaactaaa                                                                                                                        <  2:246089/141‑12 (MQ=255)
                   tattatCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGt                                                                                                                   <  2:181405/141‑1 (MQ=255)
                   tattatCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGt                                                                                                                   <  2:180358/141‑1 (MQ=255)
                    attatCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTAATAAAAATGAAAAACTAAATCAGTa                                                                                                                  <  2:610610/141‑1 (MQ=255)
                          aaTTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATc                                                                                                            >  2:132396/1‑141 (MQ=255)
                                tCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTATTTAGtttttcaactaaag                                                                                                                       <  2:507209/124‑13 (MQ=255)
                                tCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTATTTAGtttttcaactaaag                                                                                                                       <  2:507352/124‑13 (MQ=255)
                                tCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTATTTAGtttttcaactaaag                                                                                                                       >  1:507352/1‑112 (MQ=255)
                                tCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTATTTAGtttttcaactaaag                                                                                                                       >  1:507209/1‑112 (MQ=255)
                                 cAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTa                                                                                                     >  1:218729/1‑141 (MQ=255)
                                    aTGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAAtt                                                                                                  >  1:104070/1‑141 (MQ=255)
                                            aggTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGt                                                                                          <  1:666094/141‑1 (MQ=255)
                                                      gTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATg                                                                                <  1:118039/141‑1 (MQ=255)
                                                        tGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCa                                                                              >  1:504751/1‑141 (MQ=255)
                                                                cTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAAATCAGAAAGAGGAGGACATGCAAAATGAAc                                                                      >  1:83423/1‑141 (MQ=255)
                                                                          tGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTTAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAAt                                                            >  1:398406/1‑141 (MQ=255)
                                                                            cagcagGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTa                                                          <  1:484776/141‑1 (MQ=255)
                                                                             agcagGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTaa                                                         >  2:79590/1‑141 (MQ=255)
                                                                                       aaaaaCGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGc                                               >  1:490013/1‑141 (MQ=255)
                                                                                              cAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTAtt                                        <  2:442143/141‑1 (MQ=255)
                                                                                              cAGGTATACGTGCAGCAAGAATGGTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGCGGACATGAAAAATGAACCTCCATAAATTATTGTTTGCTGTGACTAtt                                        <  2:490704/141‑1 (MQ=255)
                                                                                               aGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTg                                       >  1:281333/1‑141 (MQ=255)
                                                                                                  tataCTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACATCAATAAATTAAAGTTTGGTGCGACTATTGGCa                                    >  1:518537/1‑141 (MQ=255)
                                                                                                    taCTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCatt                                  <  1:283140/141‑1 (MQ=255)
                                                                                                       ttGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCattatt                               <  1:360542/141‑1 (MQ=255)
                                                                                                            gcaAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATtggtgg                          >  1:252369/1‑141 (MQ=255)
                                                                                                            gcaAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATtggtgg                          <  1:133906/141‑1 (MQ=255)
                                                                                                              aaGAATGTTAAGTATTCAAAATCCTTCTGTATTTACAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAACGTTTGGTGCGACTATTGGCATTATtggtggtg                        <  2:295693/141‑1 (MQ=255)
                                                                                                               aGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATtggtggtgg                       >  2:355535/1‑141 (MQ=255)
                                                                                                                gAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTggtggtggt                      >  1:416438/1‑141 (MQ=255)
                                                                                                                 aaTGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTc                     <  2:5084/141‑1 (MQ=255)
                                                                                                                  aTGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCa                    >  2:699498/1‑141 (MQ=255)
                                                                                                                    gTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTg                                       <  2:144713/120‑1 (MQ=255)
                                                                                                                    gTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTg                                       >  1:144713/1‑120 (MQ=255)
                                                                                                                      tAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCtt                <  1:288425/141‑1 (MQ=255)
                                                                                                                      tAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCtt                >  2:309565/1‑141 (MQ=255)
                                                                                                                       aaGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTg               >  1:587811/1‑141 (MQ=255)
                                                                                                                       aaGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTg               <  2:134701/141‑1 (MQ=255)
                                                                                                                       aaGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTg               >  1:131395/1‑141 (MQ=255)
                                                                                                                        aGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTgg              <  2:226098/141‑1 (MQ=255)
                                                                                                                        aGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTgg              <  2:56621/141‑1 (MQ=255)
                                                                                                                            ttCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAg          >  2:216052/1‑141 (MQ=255)
                                                                                                                              cAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAgat        >  2:245965/1‑141 (MQ=255)
                                                                                                                                 aaTCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAgatgat     <  1:166015/141‑1 (MQ=255)
                                                                                                                                 aaTCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAgatgat     >  2:325935/1‑141 (MQ=255)
                                                                                                                                    ccTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAGATGATGGc  <  2:187576/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TTAAAGGGTATAGATTCTTTATTATCAATTGTTCAAATGCCAGGAGGTATTCCGGTTGCAACGACTGCAATTGGTGCAGCAGGTGCTAAAAACGCAGGTATACTTGCAGCAAGAATGTTAAGTATTCAAAATCCTTCTTTAGTTGAAAAACTAAATCAGTATGAATCTTCGTTAATTCAAAAAGTGGAGGACATGCAAAATGAACTTCAATAAATTAAAGTTTGGTGCGACTATTGGCATTATTGGTGGTGGTCAGCTTGGAAAGATGATGGC  >  CP000730/1073664‑1073936

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