Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,7610GA21.8% 111.1 / 23.2 55A154V (GCG→GTG) allSDNA‑binding transcriptional activator AllS
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (5/7);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTGGGCGA  >  NC_000913/531624‑531899
                                                                                                                                         |                                                                                                                                          
cAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGcc                                                                                                                                          >  1:690967/1‑140 (MQ=255)
 aaTATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCt                                                                                                                                         <  1:549333/140‑1 (MQ=255)
      tGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGtt                                                                                                                                    >  1:737779/1‑140 (MQ=255)
      tGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGtt                                                                                                                                    <  1:518966/140‑1 (MQ=255)
      tGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGtt                                                                                                                                    >  1:395952/1‑140 (MQ=255)
           gccgGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGt                                                                                                                               <  2:261308/140‑1 (MQ=255)
             cgGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGa                                                                                                                             >  1:625433/1‑140 (MQ=255)
                       cgcAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGa                                                                                                                             <  2:480927/130‑1 (MQ=255)
                       cgcAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGa                                                                                                                             >  1:480927/1‑130 (MQ=255)
                              gcgcTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGaa                                                                                                            >  2:379747/1‑140 (MQ=255)
                                 cTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGaaaaa                                                                                                         >  1:473826/1‑140 (MQ=255)
                                 cTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGaaaaa                                                                                                         >  2:660585/1‑140 (MQ=255)
                                  ttCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAAc                                                                                                        <  2:474552/140‑1 (MQ=255)
                                  ttCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAAc                                                                                                        <  2:473826/140‑1 (MQ=255)
                                    cTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCt                                                                                                      <  2:395952/140‑1 (MQ=255)
                                    cTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCt                                                                                                      <  2:737779/140‑1 (MQ=255)
                                    cTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCt                                                                                                      <  2:437441/140‑1 (MQ=255)
                                         aGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAAc                                                                                                        >  1:432954/1‑133 (MQ=255)
                                         aGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAAc                                                                                                        <  2:432954/133‑1 (MQ=255)
                                            ggCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACaa                                                                                              >  2:110836/1‑140 (MQ=255)
                                              ctctTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATa                                                                                            >  1:134119/1‑140 (MQ=255)
                                              ctctTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATa                                                                                            >  2:560021/1‑140 (MQ=255)
                                                          cGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGa                                                                                <  1:642819/140‑1 (MQ=255)
                                                          cGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGa                                                                                <  2:195926/140‑1 (MQ=255)
                                                           gCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGAc                                                                               >  2:433533/1‑140 (MQ=255)
                                                           gCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGAc                                                                               >  2:649574/1‑140 (MQ=255)
                                                           gCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGAc                                                                               <  1:560021/140‑1 (MQ=255)
                                                           gCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGAc                                                                               <  2:876181/140‑1 (MQ=255)
                                                            ccAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACg                                                                              <  1:379747/140‑1 (MQ=255)
                                                            ccAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACg                                                                              <  1:392805/140‑1 (MQ=255)
                                                             cAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTAGAATAGCGAGTCCCAGACGc                                                                             <  1:902506/140‑1 (MQ=255)
                                                             cAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGc                                                                             <  2:134119/140‑1 (MQ=255)
                                                                   atgatCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCatat                                                                       >  1:700061/1‑140 (MQ=255)
                                                                            cgccATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTc                                                              >  2:392230/1‑140 (MQ=255)
                                                                             gccATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACg                                                                              >  1:816104/1‑123 (MQ=255)
                                                                             gccATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACg                                                                              <  2:816104/123‑1 (MQ=255)
                                                                                 tGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGgag                                                         <  1:210155/140‑1 (MQ=255)
                                                                                        gcgcCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGaa                                                  >  2:526370/1‑140 (MQ=255)
                                                                                         cgcCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGaaa                                                 <  1:392230/140‑1 (MQ=255)
                                                                                          gcCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAAc                                                <  1:660585/140‑1 (MQ=255)
                                                                                              ttGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTggg                                            <  1:611916/140‑1 (MQ=255)
                                                                                                      aTCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAgggg                                    >  1:426237/1‑140 (MQ=255)
                                                                                                           aaGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACg                               <  1:526370/140‑1 (MQ=255)
                                                                                                                 tCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTaa                                         <  1:539706/124‑1 (MQ=255)
                                                                                                                 tCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTaa                                         >  2:539706/1‑124 (MQ=255)
                                                                                                                 tCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCAtt                         >  2:342127/1‑140 (MQ=255)
                                                                                                                 tCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCAtt                         >  1:757905/1‑140 (MQ=255)
                                                                                                                     gACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGc                     <  2:700061/140‑1 (MQ=255)
                                                                                                                         aaaGGTATTTGCCAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTAGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACg                 >  2:392514/1‑140 (MQ=255)
                                                                                                                          aaGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGc                >  1:216611/1‑140 (MQ=255)
                                                                                                                                tttGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGc                <  2:391985/134‑1 (MQ=255)
                                                                                                                                tttGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGc                >  1:391985/1‑134 (MQ=255)
                                                                                                                                 ttGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAAc         <  2:757905/140‑1 (MQ=255)
                                                                                                                                     cAGCACCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTAGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTggg     <  1:392514/140‑1 (MQ=255)
                                                                                                                                      aGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTGGGc    >  1:995800/1‑140 (MQ=255)
                                                                                                                                       gcgcCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTGGGCg   <  1:649574/140‑1 (MQ=255)
                                                                                                                                        cgcCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTGGGCGa  >  1:290217/1‑140 (MQ=255)
                                                                                                                                        cgcCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTGCGCGa  >  1:596359/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTATTTGCCAGCGCCTCAGTTCCCGTGACGCCGATAGCCAGCGAAAAACCTTCGTACAATAGCGAGTCCCAGACGCCCATATAGATTTGTCGGGAGATGTGAAACTGGGTAAAGGGGTAACGCTCATTCAGCCACGCCAGCAACTGGGCGA  >  NC_000913/531624‑531899

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