Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913958,2130CT73.7% 55.4 / 24.2 39S187S (AGC→AGT) serC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base T (14/14);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATC  >  NC_000913/958080‑958342
                                                                                                                                     |                                                                                                                                 
ctgctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTa                                                                                                                              >  2:454399/1‑139 (MQ=255)
ctgctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTa                                                                                                                              <  1:286263/139‑1 (MQ=255)
ctgctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTa                                                                                                                              >  1:83615/1‑139 (MQ=255)
ctgctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTa                                                                                                                              >  1:400824/1‑139 (MQ=255)
ctcctTATATGCATTATTGCCCGAATGAAACCATCGATGGTAACGCCACCGACGCAACGCCAGACTTCGGCGCGGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTTACGTGAGGCGtca                                                                                                                              <  2:303862/136‑3 (MQ=255)
    ttATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTGTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGgt                                                                                                                          <  1:637970/139‑1 (MQ=255)
    ttATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGgt                                                                                                                          >  1:565161/1‑139 (MQ=255)
        atGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTaa                                                                                                                      >  2:641093/1‑139 (MQ=255)
               attGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGc                                                                                                               >  2:38101/1‑139 (MQ=255)
                                  cGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATc                                                                                            >  2:88863/1‑139 (MQ=255)
                                      ggTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGcc                                                                                        >  2:152497/1‑139 (MQ=255)
                                               aTCGACGAAACGCCAGACTTGGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTggac                                                                               >  2:420660/1‑137 (MQ=255)
                                                       aaCGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGgt                                                                                                                          >  1:274973/1‑88 (MQ=255)
                                                       aaCGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGgt                                                                                                                          <  2:274973/88‑1 (MQ=255)
                                                                   cGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCAtcgttcgt                                                           >  1:212031/1‑139 (MQ=255)
                                                                   cGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCAtcgttcgt                                                           >  1:687238/1‑139 (MQ=255)
                                                                       gcAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAg                                                       <  1:707119/139‑1 (MQ=255)
                                                                            tgtgGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTg                                                  <  1:454399/139‑1 (MQ=255)
                                                                                  ggtCGCCGCTGACTTCTCTTCTACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGc                                            <  1:420660/139‑1 (MQ=255)
                                                                                   gtCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCa                                           >  2:685554/1‑139 (MQ=255)
                                                                                     cgccgcTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCaaa                                         >  1:30095/1‑139 (MQ=255)
                                                                                        cgcTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTa                                                                                                                              >  1:179274/1‑51 (MQ=255)
                                                                                        cgcTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTa                                                                                                                              <  2:179274/51‑1 (MQ=255)
                                                                                              cTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGc                                            >  1:331226/1‑127 (MQ=255)
                                                                                              cTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGc                                            <  2:331226/127‑1 (MQ=255)
                                                                                               ttCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATcg                               <  1:38101/139‑1 (MQ=255)
                                                                                                  tctTCAACCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCgt                            >  2:365986/1‑139 (MQ=255)
                                                                                                    ttCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCgtgt                          >  1:570603/1‑139 (MQ=255)
                                                                                                       aaCCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCg                       >  1:644566/1‑139 (MQ=255)
                                                                                                        aCCATTCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGt                      <  2:212031/139‑1 (MQ=255)
                                                                                                            ttCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGAt                  >  1:485532/1‑139 (MQ=255)
                                                                                                             tCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGAtt                 <  1:641093/139‑1 (MQ=255)
                                                                                                             tCTTTCCCGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGAtt                 <  2:485532/139‑1 (MQ=255)
                                                                                                             tCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGAtt                 <  1:88863/139‑1 (MQ=255)
                                                                                                                  cccGTCCGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGa            <  2:565161/139‑1 (MQ=255)
                                                                                                                  cccGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGa            <  2:83615/139‑1 (MQ=255)
                                                                                                                       ccGATTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGattatt       >  2:444097/1‑139 (MQ=255)
                                                                                                                          aTTGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAAACGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCa    <  1:444097/139‑1 (MQ=255)
                                                                                                                            tGACGTCAGTCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATc  <  2:432445/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                 
CTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATC  >  NC_000913/958080‑958342

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