Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,692,9510TC41.9% 54.0 / 24.5 43I417V (ATC→GTC) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (24/1);  new base C (0/18);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 3.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCC  >  NC_007779/2692852‑2693051
                                                                                                   |                                                                                                    
cGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATc                                                                                                     <  3:424216/101‑1 (MQ=255)
cGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATc                                                                                                     >  1:1157210/1‑101 (MQ=255)
 gCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATcc                                                                                                    >  1:967669/1‑101 (MQ=255)
  ccAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATccc                                                                                                   >  1:430059/1‑101 (MQ=255)
    aaGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATcccgc                                                                                                 >  3:379958/1‑101 (MQ=255)
                  ttGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATc                                                                                   >  1:1205055/1‑101 (MQ=255)
                   tGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCg                                                                                  >  3:298606/1‑101 (MQ=255)
                   tGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCg                                                                                  >  1:604933/1‑101 (MQ=255)
                   cgCCGGCCCGCCGAGACCGCCCAGCTCCGCCCCGACGTTGTTCTCGCTTTTTTGTGCGTGATCGGCGCAAATGTTGCCGACCCCGCCCGCCAGCATGATCg                                                                                  <  1:942654/100‑1 (MQ=255)
                     ccGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGt                                                                                >  1:607614/1‑101 (MQ=255)
                     ccGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGt                                                                                >  3:792396/1‑101 (MQ=255)
                         gcgccggagaCCGCCCAGCTTCGCCCCGATGTTGTTCTCGCTTTTTTGTAGTGGTTGGGCGCGTGTTTTGCCTACCCCGCCCGCCAGCATGATCGGTTtgt                                                                            <  1:1481276/95‑1 (MQ=255)
                         gccgccGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTtgt                                                                            >  3:876496/1‑101 (MQ=255)
                         cccgcggAGACCGTCCAGCTTCGCTCCGGCGTTTTTCTCTCTTTTTTGTACGTGATGGGCGCAAATTTTGCCGACCCCGCCCGCCAGCATGATCGGTTtgt                                                                            <  1:726913/100‑1 (MQ=255)
                               gagaACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAAc                                                                      >  3:510788/1‑101 (MQ=255)
                                aaatctaCCAGCTCCACACCGACGTTTATCTCATTTTTTTATACGTAATGGCCACTAATGTTCCCGACCCCGCCCGCCAGCATGATCGGTTTGTGATAAcc                                                                     <  1:241118/95‑1 (MQ=255)
                                 gaACGTCCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACcg                                                                    >  1:366748/1‑101 (MQ=255)
                                  aaCGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACcgc                                                                   >  3:1499747/1‑101 (MQ=255)
                                    cGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCgcgc                                                                 >  1:363772/1‑101 (MQ=255)
                                     gACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCa                                                                >  3:1086490/1‑101 (MQ=255)
                                        cAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGct                                                             >  1:889845/1‑101 (MQ=255)
                                        cAGCTCCGCACCGCCGTTGTTCTCGCTTTTTTGTGGGGGATGGGCGCGAATGTTGCCGCCCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGct                                                             <  3:173818/101‑1 (MQ=255)
                                            ttCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTc                                                         >  3:1522254/1‑101 (MQ=255)
                                               gCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGcc                                                      >  3:680852/1‑101 (MQ=255)
                                               cccccgttgttgttCTCTCTTTTTTGTAGGGGAGCGGCGCGATTGTTCCCGCCCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGcc                                                      <  3:1192825/93‑1 (MQ=255)
                                               cccccGACGTTGATCTCCCTTTTTTGTACGTGATGGGCGCAAATGTTGCCGCCCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGcc                                                      <  3:1282272/98‑1 (MQ=255)
                                                      cGTTGTTCTCGCTTTTTTGTGCGGGAGCGGCGCGAATGTTGCCGACCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTgg                                               <  3:1483091/101‑1 (MQ=255)
                                                      cGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTgg                                               >  1:734837/1‑101 (MQ=255)
                                                        ttGATCTCGCTTTTTTTTGCGTGGTGGGGGGGATTTTTCCCGACCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCt                                             <  3:191197/101‑1 (MQ=255)
                                                         ttttctcCCTTTTTTGTGGGGGAGCGGGGCGAATGTCCCCGCCCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTg                                            <  1:377076/98‑1 (MQ=255)
                                                           aTCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGtt                                          >  1:380252/1‑101 (MQ=255)
                                                              tcgtttttttGTCGGTGGGCGGCGCGATTGTTGCCGCCCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCAc                                       <  3:1384416/96‑1 (MQ=255)
                                                                     ttgggaCGTGGCCGGCGCGAGTGTGGCCGACCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTttct                                <  1:1008818/96‑1 (MQ=255)
                                                                       tgtgCGTGAGCGGCGCGATTGTGCCCGCCCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTttcttc                              <  3:1045681/101‑1 (MQ=255)
                                                                       tggaCGTGATCGGCGCGATTGTTGCGGCCCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTttcttc                              <  3:444902/98‑1 (MQ=255)
                                                                           ggTGAGCGGCGGGAGTGTCCCCGCCCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATaa                          <  3:159999/100‑1 (MQ=255)
                                                                           cGTGATCGGCGCGAATGTTGCCGATCCCGCCCCCCAGCATGATCGGTTTGTGATAACCGGGCGCTTCTTCCCCGTTGTGGCTGTTACCTTTTTTTTCATaa                          >  1:458747/1‑101 (MQ=255)
                                                                                 cGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACgg                    >  1:1547736/1‑101 (MQ=255)
                                                                                  gccgcgAATTTTGCCCCCCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGa                   <  1:79756/99‑1 (MQ=255)
                                                                                      ggagtgtgGCCGCCCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTa               <  3:184072/97‑1 (MQ=255)
                                                                                        aaTGTTCCCGACCCCCCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGc             <  3:977194/101‑1 (MQ=255)
                                                                                                  aTCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGc   >  3:1050384/1‑101 (MQ=255)
                                                                                                  aTCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGc   >  3:796069/1‑101 (MQ=255)
                                                                                                   tCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGcc  >  1:191015/1‑101 (MQ=255)
                                                                                                   tCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGcc  >  1:1121657/1‑101 (MQ=255)
                                                                                                   |                                                                                                    
CGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCC  >  NC_007779/2692852‑2693051

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑