Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,076,774 T→C 19.2% W51R (TGG→CGG)  purL → phosphoribosylformylglycinamidine synthase subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,076,7740TC19.2% 59.5 / 8.5 32W51R (TGG→CGG) purLphosphoribosylformylglycinamidine synthase subunit II
Reads supporting (aligned to +/- strand):  ref base T (12/13);  new base C (3/3);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATCGAACCAAGCGTTGAAGAAATTAAACTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCGTAGTATTTAAA  >  CP000730/1076636‑1076911
                                                                                                                                          |                                                                                                                                         
aTCGAACCAAACGATGAAGAAACTAAACTTGAAAAAGTAGATCCAGATATGGGATTAAGTGAACAAGAATATGAAAAAGATTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTgg                                                                                                                                         <  2:285593/141‑1 (MQ=255)
               gaagaaATTAAACTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTATTCTGTTATGTGGAGTGAACATTGctct                                                                                                                          >  1:249074/1‑141 (MQ=255)
                         aaCTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTACTCTTATAAACAtt                                                                                                                >  1:33111/1‑141 (MQ=255)
                          aCTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCGCTTATAAACATTc                                                                                                               >  1:284293/1‑141 (MQ=255)
                                aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAAcc                                                                                                         <  2:326103/141‑1 (MQ=255)
                                aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAAcc                                                                                                         >  2:79487/1‑141 (MQ=255)
                                aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAAcc                                                                                                         >  2:320932/1‑141 (MQ=255)
                                aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAAcc                                                                                                         >  2:28420/1‑141 (MQ=255)
                                     tataTCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGtttt                                                                                                    >  1:652/1‑141 (MQ=255)
                                          cAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAg                                                                                               <  2:3971/141‑1 (MQ=255)
                                          aaaGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAg                                                                                               <  2:3975/140‑1 (MQ=255)
                                               tatGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGTCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAAtt                                                                                          <  2:354207/141‑1 (MQ=255)
                                                atGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAAttt                                                                                         <  2:130164/141‑1 (MQ=255)
                                                atGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAAttt                                                                                         <  1:215333/141‑1 (MQ=255)
                                                 tGGGACTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACACAATCAGGTATCTTTTCTGATATCGGCAGTGAACAAAGCTCTTATAAACATTCTAAACCGCTTTTCACGCAATTTc                                                                                        >  2:247063/1‑141 (MQ=255)
                                                         aGTGAACAAGAAGATGAAAAAGTTTGCGATAGTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCa                                                                                <  1:146348/141‑1 (MQ=255)
                                                               cAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGAc                                                                          >  1:213828/1‑141 (MQ=255)
                                                                      atGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGc                                                                   >  2:159550/1‑141 (MQ=255)
                                                                      atGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGc                                                                   <  1:120089/141‑1 (MQ=255)
                                                                             aGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATggg                                                            >  2:90999/1‑141 (MQ=255)
                                                                                   cGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTTAATTTATTTCCTACGTCAGGTGGCCATGTGCTTATGGGGCCTgg                                                      >  2:89716/1‑141 (MQ=255)
                                                                                      tatTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGa                                                   <  2:139700/141‑1 (MQ=255)
                                                                                                aGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTCTGGGGCCTGGTGAAGGTGCAggg                                         <  2:256221/141‑1 (MQ=255)
                                                                                                    aaCCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAg                                     <  1:28420/141‑1 (MQ=255)
                                                                                                    aaCCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAg                                     <  1:79487/141‑1 (MQ=255)
                                                                                                    aaCCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTCCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAg                                     <  1:320932/141‑1 (MQ=255)
                                                                                                       cTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCg                                  >  2:115992/1‑141 (MQ=255)
                                                                                                                      cAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAAGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTATTCGATATAGGTGatcatc                   >  2:152505/1‑141 (MQ=255)
                                                                                                                      aaGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATc                   <  1:308289/140‑1 (MQ=255)
                                                                                                                             cTTTTCTGTTATGTGGAGTGAACGTTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCgg            >  2:256472/1‑140 (MQ=255)
                                                                                                                              ttttCTTTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCgta           <  2:325688/141‑1 (MQ=255)
                                                                                                                                       aTGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCGTAGTATTTaaa  >  1:337946/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
ATCGAACCAAGCGTTGAAGAAATTAAACTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCGTAGTATTTAAA  >  CP000730/1076636‑1076911

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: