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→AGAydcRputative 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

CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTT  >  NC_000913/1510249‑1510532
                                                                                                                                                     |                                                                                                                                         
cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa                                                                                                                                       <  1:340999/151‑1 (MQ=255)
    aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGc                                                                                                                                   >  1:163368/1‑151 (MQ=255)
               gATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGAt                                                                                                                          >  2:624537/1‑149 (MQ=255)
                            atatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGATCGTTCCGCTTGTCTCGGCCTTTCCCGACCAGCGACTTTTTGCCTTCTAACTACTAAACCGCTCGCGGGCGCAGGACCGCAAAACCGCCACAGCGATGATCGTGAt                                                                                                                 >  2:408216/1‑145 (MQ=255)
                                    gATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGc                                                                                                   >  1:587481/1‑151 (MQ=255)
                                     atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc                                                                                                  >  1:404812/1‑151 (MQ=255)
                                     atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc                                                                                                  >  2:372871/1‑151 (MQ=255)
                                        tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                               >  1:490923/1‑151 (MQ=255)
                                             cAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                          >  1:16514/1‑151 (MQ=255)
                                               gcccAGCCGCGTCCCGTGGGTTGTCCCGTTTCCCCGGGCTTTCCCCGCCCCGCGTTTTTCCCCCCCCAAAAAACAAACCCGCCCGGGGGCGAAGTTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                         <  1:91398/148‑1 (MQ=255)
                                               gcccAGCCGCGCCCCGTCGGTCTTCCTTTTCCCCCCGCCTTTCCCCCCCCCGCGTTTTTCCCCCCCCAAAAAACAAACCCCCCCGCGGGCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                         <  1:662267/148‑1 (MQ=255)
                                                 ccggccgcgccccGTGGGTCGTCCCTTTTCCCTGGCCTTTCCCCGCCCGGCGCTTTTCCCCCCCCCAAAAACAAACCCCCCCGCTGGCGCAGGTAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAg                                                                                       <  1:31341/139‑1 (MQ=255)
                                                 ccagccccaccccgccgGTCGTCCCTTTGCCCCCGCCCTTCCCCCCCCCGCGTTTTTCCCCCCCCAAACAACAAACCCGCTCGGGGCCGAAGATAAGAAAATCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                          <  2:566616/136‑1 (MQ=255)
                                                              gtcggtcgTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt                                                                         >  2:415918/1‑151 (MQ=255)
                                                               tcggtcgTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGa                                                                                                              >  1:607887/1‑113 (MQ=255)
                                                                cggtcgTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                       >  2:47464/1‑151 (MQ=255)
                                                                          ttttGCCCCGGCCTTTCCCCGCCCCGCACTTTCCCCCCCCCCAAAAAAAAACCCGCTCGGGGGCGAGGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                       <  2:628437/141‑1 (MQ=255)
                                                                           ttggCCGCGCCCTTCCCCCCCCCGCATTTTTCCCCCCCCAAAAAAAAAACCCTCCCCTTGGCGGAGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtcgt                                                                <  1:517623/144‑1 (MQ=255)
                                                                            ttccccccgCCTTTCCCCGCCCCGCTTTTTTCCCCCCCCCAAAAACAAACCCGCTCGCGGGCGAGGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGc                                                                          <  1:706297/129‑1 (MQ=255)
                                                                                   ggCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAACCCGCTCGCTGGCGCAGGTAAGCAACACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt                                                     >  2:285066/1‑150 (MQ=255)
                                                                                    gCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc                                                   >  2:233736/1‑151 (MQ=255)
                                                                                      cTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc                                                   >  1:298582/1‑149 (MQ=255)
                                                                                              ccccccgctttttccccccccccAACACCAAACCCGCCCGCGGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCAt                                            <  1:469412/135‑1 (MQ=255)
                                                                                                   cgcactttcccccccTCCAACACAAACCCCGCGCGCGGCGGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc                                        <  2:423392/139‑1 (MQ=255)
                                                                                                       cTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc                                 >  1:244062/1‑151 (MQ=255)
                                                                                                         ttccccccccccAACAAAAAACCCCCCCGCGGCCGCAGAAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtct                                  <  2:352438/145‑1 (MQ=255)
                                                                                                         ttccccccccCAAACAACAAACCCGCCCGCGGGCGCAGGAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtc                                   <  1:488018/144‑1 (MQ=255)
                                                                                                             ccccccCAAACAACAAACCCGCCCGTGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                            <  1:520320/151‑1 (MQ=255)
                                                                                                             ccccccCAAACAACAAACCCGCCCGGGGGCGCAGAAAACAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                            <  2:201185/151‑1 (MQ=255)
                                                                                                             ccccccCAAACAAAAAACCCCCCCGCGGGCGAAGAAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAttt                                    <  2:536805/143‑1 (MQ=255)
                                                                                                                ccTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga                          >  2:344070/1‑150 (MQ=255)
                                                                                                                    caacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                     >  2:379102/1‑151 (MQ=255)
                                                                                                                        aacTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAt                 >  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)
                                                                                                                                                     |                                                                                                                                         
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTT  >  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: