Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,481,631 G→A 7.7% A320A (GCG→GCA)  purH → Bifunctional purine biosynthesis protein PurH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,481,6310GA 7.7% 168.7 / 3.1 52A320A (GCG→GCA) purHBifunctional purine biosynthesis protein PurH
Reads supporting (aligned to +/- strand):  ref base G (33/15);  new base A (2/2);  total (35/17)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

AAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACGTGGTGGCTGCCAAA  >  ACB143/5481499‑5481763
                                                                                                                                    |                                                                                                                                    
aaaGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTc                                                                                                                                   >  1:218702/1‑136 (MQ=255)
aaaGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTc                                                                                                                                   >  1:64121/1‑136 (MQ=255)
aaaGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTc                                                                                                                                   >  4:284778/1‑136 (MQ=255)
aaaGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCATTc                                                                                                                                   >  5:192960/1‑136 (MQ=255)
aaaGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCATTc                                                                                                                                   >  5:192958/1‑136 (MQ=255)
  agagCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTcgg                                                                                                                                 >  4:138694/1‑136 (MQ=255)
         cGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTGCGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGcat                                                                                                                           >  4:160737/1‑135 (MQ=255)
         cGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGca                                                                                                                            <  3:126886/134‑1 (MQ=255)
            cAAGCCGCCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACAGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACAGACACCGAATCGGCGTTCGGCGGCaacatc                                                                                                                       >  1:299539/1‑136 (MQ=255)
              aGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGc                                                                                                                     <  3:111801/136‑1 (MQ=255)
              aGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGGACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGc                                                                                                                     <  2:11214/136‑1 (MQ=255)
                 cGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCAtt                                                                                                                  >  7:41315/1‑136 (MQ=255)
                   gCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCa                                                                                                                <  5:293856/136‑1 (MQ=255)
                      tGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAAcc                                                                                                             >  8:71607/1‑136 (MQ=255)
                        cGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGt                                                                                                           <  6:172047/136‑1 (MQ=255)
                                 cAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGAt                                                                                                  >  7:6629/1‑136 (MQ=255)
                                    gCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGc                                                                                               <  8:148273/136‑1 (MQ=255)
                                     cACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCg                                                                                              <  3:211799/136‑1 (MQ=255)
                                     cACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCg                                                                                              <  1:41234/136‑1 (MQ=255)
                                          cAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAAcc                                                                                         <  5:133973/136‑1 (MQ=255)
                                           aaCCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAACCCg                                                                                        >  7:89480/1‑136 (MQ=255)
                                                     gcgTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCAt                                                                              >  3:209081/1‑136 (MQ=255)
                                                     gcgTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCAt                                                                              >  4:10402/1‑136 (MQ=255)
                                                     gcgTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCAt                                                                              >  7:140335/1‑136 (MQ=255)
                                                     gcgTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCAt                                                                              >  4:6211/1‑136 (MQ=255)
                                                     gcgTGGCAGTAGTAACGGAAGACGAAGGCGGGATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCAt                                                                              >  8:222073/1‑136 (MQ=255)
                                                           cAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGa                                                                        >  5:69222/1‑136 (MQ=255)
                                                                   ccGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGc                                                                                                                     <  1:274021/83‑1 (MQ=255)
                                                                   ccGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGc                                                                                                                     >  2:274021/1‑83 (MQ=255)
                                                                      gAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTtcg                                                             >  6:262076/1‑136 (MQ=255)
                                                                         gACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTtcgtc                                                           >  4:291946/1‑135 (MQ=255)
                                                                              aGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGt                                                      >  5:1972/1‑135 (MQ=255)
                                                                                                      ccTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGTCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTAGAAGTGCTCATCGCCCCGAAAATCTccccg                             >  6:210128/1‑134 (MQ=255)
                                                                                                      ccTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAg                             >  8:183023/1‑136 (MQ=255)
                                                                                                            cTACGCGACCGACACCGACTCGGCGTTCGGCGGCATCATCGCAAACAACCGTGAACAGGAAGCCGAAACCGCCAAGGCCATAGACGACCGCCAGTTCCACCAACTGAACATCGCCCCCAAAAACACCCAGgccgcc                       >  6:63293/1‑136 (MQ=255)
                                                                                                            cTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAACCCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGgccgcc                       >  3:53393/1‑136 (MQ=255)
                                                                                                            cTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGgccgcc                       >  4:239625/1‑136 (MQ=255)
                                                                                                            cTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGgccgcc                       >  2:165371/1‑136 (MQ=255)
                                                                                                            cTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGgccgcc                       <  2:269154/136‑1 (MQ=255)
                                                                                                            cTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGgccgcc                       >  3:30724/1‑136 (MQ=255)
                                                                                                            cTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGgccgcc                       >  7:183037/1‑136 (MQ=255)
                                                                                                                    ccgacACCGTATAGGCGGTCGGCGGAAGCTTCGCATTCAACGGTGAACTGGATGGCGAAACAGCCAAGGCCATTGTCGACCGACTGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGACCGCGACgt               <  8:53205/136‑1 (MQ=255)
                                                                                                                      gacACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACgtgg             <  8:140335/136‑1 (MQ=255)
                                                                                                                      gacACCGAATCGGCATTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACgtgg             <  6:192960/136‑1 (MQ=255)
                                                                                                                      gacACCGAATCGGCATTCGGCGGCATCATCGCATTCAACCGTCAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACgtgg             <  6:192958/136‑1 (MQ=255)
                                                                                                                       acacCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAACCCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACgtggt            >  6:21492/1‑136 (MQ=255)
                                                                                                                       acacCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACgtggt            >  6:278486/1‑136 (MQ=255)
                                                                                                                       acacCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGAGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACttggt            >  6:87355/1‑136 (MQ=255)
                                                                                                                         acCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACGtggtgg          <  6:267519/136‑1 (MQ=255)
                                                                                                                          ccGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACGTGGTGGc         >  6:122540/1‑136 (MQ=255)
                                                                                                                                 ggCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACGTGGTGGCTGCCaaa  >  2:20697/1‑136 (MQ=255)
                                                                                                                                 ggCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACGTGGTGGCTGCCaa   <  8:61581/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
AAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCATCGCATTCAACCGTGAACTGGATGGCGAAACCGCCAAGGCCATTGTCGACCGCCAGTTCGTCGAAGTGATCATCGCCCCGAAAATCTCCCAGGCCGCCCGCGACGTGGTGGCTGCCAAA  >  ACB143/5481499‑5481763

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