Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA3,817,9490AC20.6% 57.1 / 12.2 34I64L (ATC→CTC) dinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base A (14/13);  new base C (5/2);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTT  >  NC_000913_3_pae_tpiA/3817813‑3818086
                                                                                                                                        |                                                                                                                                         
aGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATgg                                                                                                                               <  1:246993/149‑1 (MQ=255)
 gCATAAAAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGATAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGt                                                                                                                              <  2:106773/149‑1 (MQ=255)
      aCAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATGCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCaaaa                                                                                                                         <  1:137866/149‑1 (MQ=255)
         aTGGTTAGCAATACGACTTGGGAAAGTCCTGGGTTATTCAGGGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGGCTGTGAGAACAGTGGTCACAGAATTGATTGCCATTTCGAGGAGATGCTCGATATGGTTAAAATTg                                                                                                                      <  2:485760/149‑1 (MQ=255)
                  aaGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACTCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGCTCCTCGATATGGTCAAAATTGGCTCAAATg                                                                                                             >  1:97583/1‑149 (MQ=255)
                           tGGGAAACTACTGGGCTATTCAGAGTAACGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAagag                                                                                                    <  2:132874/149‑1 (MQ=255)
                            gggAAAGCACTGGGGTATACAGAGTATCGTCACTGTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGc                                                                                                   <  2:271454/149‑1 (MQ=255)
                                      tGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCAAAAAGAGCATTAAAAGAc                                                                                         >  1:460039/1‑149 (MQ=255)
                                        ggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACAt                                                                                       >  2:312133/1‑149 (MQ=255)
                                         cttATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATc                                                                                      <  2:216011/148‑1 (MQ=255)
                                               gagagTATAGTCACTTTATACCTGTATTAACGCGGGCAAGAGAGGCGGGTGAGAAGAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTActc                                                                                <  2:353603/148‑1 (MQ=255)
                                                                   ccTGTATTAACGCGCGCCAACGAAGCCTGTGAGAACAGTGGTCACACAAGTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCACATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTCTGCCTGTTAcca                                                            >  2:120041/1‑148 (MQ=255)
                                                                   ccTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCt                                                            <  1:312133/149‑1 (MQ=255)
                                                                       tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgta                                                        <  2:170926/149‑1 (MQ=255)
                                                                                    cAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCg                                           >  1:142931/1‑149 (MQ=255)
                                                                                            ccTGTGAGAACAGTGGTCACACAATCGATGCAAGGATCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGca                                   >  1:246266/1‑148 (MQ=255)
                                                                                              tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaa                                 >  1:175752/1‑149 (MQ=255)
                                                                                              tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaa                                 >  2:163663/1‑149 (MQ=255)
                                                                                               gtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                >  2:429300/1‑149 (MQ=255)
                                                                                                  agaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                             >  1:336009/1‑149 (MQ=255)
                                                                                                  agaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                             >  1:336001/1‑149 (MQ=255)
                                                                                                  agaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                             <  1:206325/149‑1 (MQ=255)
                                                                                                   gaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCgg                            <  1:228763/149‑1 (MQ=255)
                                                                                                   gaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATACTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCgg                            >  1:349772/1‑149 (MQ=255)
                                                                                                            tCACACAATTGAGGACACTTCCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAAGGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCgg                   <  1:400380/149‑1 (MQ=255)
                                                                                                            tCACACAATTGACGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCgg                   <  1:400373/149‑1 (MQ=255)
                                                                                                                    attttgaCCATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggcag           >  2:172603/5‑149 (MQ=39)
                                                                                                                     ttttgaCCATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTACTTGCGGCAGGGCAGa          >  2:314790/4‑149 (MQ=38)
                                                                                                                         gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCCGCCTTa      >  2:42974/1‑149 (MQ=255)
                                                                                                                            cATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAttt   >  2:337652/1‑149 (MQ=255)
                                                                                                                            cATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAttt   >  2:402132/1‑149 (MQ=255)
                                                                                                                             aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtttt  >  1:304096/1‑149 (MQ=255)
                                                                                                                                tCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                          <  2:102578/90‑1 (MQ=39)
                                                                                                                                tCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                          >  1:102578/1‑90 (MQ=39)
                                                                                                                                        |                                                                                                                                         
AGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTT  >  NC_000913_3_pae_tpiA/3817813‑3818086

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: ‑