Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3950CA44.4% 33.4 / 25.9 36S131R (AGC→AGA) ydcRputative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (19/1);  new base A (0/16);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 2.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTT  >  NC_000913/1510249‑1510532
                                                                                                                                                     |                                                                                                                                         
cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCaaaa                                                                                                                                       <  1:340999/151‑1 (MQ=255)
    aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGc                                                                                                                                   >  1:163368/1‑151 (MQ=255)
               gATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGAt                                                                                                                          >  2:624537/1‑149 (MQ=255)
                            atatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGATCGTTCCGCTTGTCTCGGCCTTT‑CCCGACCAGCGACTTT‑TTGCCTTCTAACTACTAAA‑CCGCTCGCGGGCGCAGGACCGCAAAACCGCCACAGCGATGATCGTGAt                                                                                                                 >  2:408216/1‑145 (MQ=255)
                                    gATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGc                                                                                                   >  1:587481/1‑151 (MQ=255)
                                     atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc                                                                                                  >  1:404812/1‑151 (MQ=255)
                                     atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc                                                                                                  >  2:372871/1‑151 (MQ=255)
                                        tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                               >  1:490923/1‑151 (MQ=255)
                                             cAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                          >  1:16514/1‑151 (MQ=255)
                                               gcccAGCCGCGTCCCGTGGGTTGTCCCGTTTCCCCGGGCTTTC‑CCCGCCCCGCGTTTTT‑CCCCCCCCAAAAAACAAAC‑CCGCCCGGGGGCGAAGTTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                         <  1:91398/148‑1 (MQ=255)
                                               gcccAGCCGCGCCCCGTCGGTC‑TTCCTTTTCCCCCCGCCTTTCCCCCCCCCGCGTTTTT‑CCCCCCCCAAAAAACAAAC‑CCCCCCGCGGGCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                         <  1:662267/148‑1 (MQ=255)
                                                 ccggccgcgccccGTGGGTCGTCCCTTTTCCCT‑GGCCTTTCCCCGCCCGGCGCTTTT‑CCCCCCCCCAAAAACAAAC‑CCCCCCGCTGGCGCAGGTAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAg                                                                                       <  1:31341/139‑1 (MQ=255)
                                                 ccagccccaccccgccgGTCGTCCCTTTGCCCCCGCCCTTC‑CCCCCCCCGCGTTTTT‑CCCCCCCCAAACAACAAAC‑CCGCTCGGGGCCGAAGATAAGAAAATCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                          <  2:566616/136‑1 (MQ=255)
                                                              gtcggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt                                                                         >  2:415918/1‑151 (MQ=255)
                                                               tcggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGa                                                                                                              >  1:607887/1‑113 (MQ=255)
                                                                cggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                       >  2:47464/1‑151 (MQ=255)
                                                                          ttttGCCCCGGCCTTTCCCCGCCCCGC‑ACTTTCCCCCCCCCCAAAAAAAAAC‑CCGCTCGGGGGCGAGGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                       <  2:628437/141‑1 (MQ=255)
                                                                           ttggCCGCGCCCTTC‑CCCCCCCCGC‑ATTTT‑TCCCCCCCCAAAAAAAAAACCCTCCCCTTGGCGGAGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtcgt                                                                <  1:517623/144‑1 (MQ=255)
                                                                            ttccccccgCCTTTCCCCGCCCCGCTTTTTT‑CCCCCCCCCAAAAACAAAC‑CCGCTCGCGGGCGAGGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGc                                                                          <  1:706297/129‑1 (MQ=255)
                                                                                   ggCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAC‑CCGCTCGCTGGCGCAGGTAAGCAACACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt                                                     >  2:285066/1‑150 (MQ=255)
                                                                                    gCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc                                                   >  2:233736/1‑151 (MQ=255)
                                                                                      cTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc                                                   >  1:298582/1‑149 (MQ=255)
                                                                                              ccccccgctttttccccccccccAACACCAAAC‑CCGCCCGCGGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCAt                                            <  1:469412/135‑1 (MQ=255)
                                                                                                   cgcactttcccccccTCCAACACAAACC‑CCGCGCGCGGCGGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc                                        <  2:423392/139‑1 (MQ=255)
                                                                                                       cTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc                                 >  1:244062/1‑151 (MQ=255)
                                                                                                         ttccccccccccAACAAAAAAC‑CCCCCCGCGGCCGCAGAAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtct                                  <  2:352438/145‑1 (MQ=255)
                                                                                                         ttccccccccCAAACAACAAAC‑CCGCCCGCGGGCGCAGGAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtc                                   <  1:488018/144‑1 (MQ=255)
                                                                                                             ccccccCAAACAACAAAC‑CCGCCCGTGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                            <  1:520320/151‑1 (MQ=255)
                                                                                                             ccccccCAAACAACAAAC‑CCGCCCGGGGGCGCAGAAAACAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                            <  2:201185/151‑1 (MQ=255)
                                                                                                             ccccccCAAACAAAAAAC‑CCCCCCGCGGGCGAAGAAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAttt                                    <  2:536805/143‑1 (MQ=255)
                                                                                                                ccTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga                          >  2:344070/1‑150 (MQ=255)
                                                                                                                    caacaaCTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                     >  2:379102/1‑151 (MQ=255)
                                                                                                                        aacTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAt                 >  2:12061/1‑151 (MQ=255)
                                                                                                                                        ggtcaaagtaaaGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTGCCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg         <  1:549383/134‑1 (MQ=255)
                                                                                                                                         ggCGCAGGAAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCg    <  1:75446/148‑1 (MQ=255)
                                                                                                                                           cgcAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt  >  1:188617/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTT  >  NC_000913/1510249‑1510532

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