Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4860CA37.8% 80.6 / 21.9 45intergenic (+291/‑165)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (24/4);  new base A (0/17);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 5.43e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29336‑29635
                                                                                                                                                         |                                                                                                                                                     
tttttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCcac                                                                                                                                                       >  2:124369/1‑151 (MQ=255)
         ttttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAaataat                                                                                                                                               >  2:661046/1‑150 (MQ=255)
          tttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAataata                                                                                                                                              <  1:1093/150‑1 (MQ=255)
              aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATaaa                                                                                                                                          >  1:577989/1‑150 (MQ=255)
                gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATaaaaa                                                                                                                                        >  2:18715/1‑150 (MQ=255)
                    aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATccc                                                                                                                                   >  2:433988/1‑151 (MQ=255)
                    aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATccc                                                                                                                                   >  1:613654/1‑151 (MQ=255)
                      tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                  <  1:124332/150‑1 (MQ=255)
                                 ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                       >  2:793037/1‑150 (MQ=255)
                                    cctcTTTTTTCTCCTTTTTTTTATTTCATCCAATCTTCTTGCTGCGCAAGCTTTTCCAAAAACGG‑TTT‑AAATTATTTTTTTGTCCCTAAATCCCTTAAAAACATGCA‑TGACCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc                                                                                                                     <  2:139244/146‑1 (MQ=255)
                                    ccATTTTTTTCTTTTTTTTTTTATTTCATGCAATCTTCTTGCTGCGAAAGC‑TTTTCCAAAAACG‑GTTAAAATAATTTTTTTGTCTTTTAATCCCTTTAAAACATGCA‑TGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACt                                                                                                                    <  1:282949/150‑1 (MQ=255)
                                        ttttttttGCTCTTTTTTTTTTAATGCATTTTTTTTGCTCCAAATGCTTTTTCAAAAAAAG‑GTA‑AGATAATTTTTTTGTCGTTTATTCCCGAAAAAAGTTGCT‑TGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                 <  1:169831/149‑1 (MQ=255)
                                          ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                             >  2:52126/1‑151 (MQ=255)
                                                    gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   >  1:279304/1‑151 (MQ=255)
                                                     tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   >  1:754187/1‑150 (MQ=255)
                                                      ttttATTTCATGCATTTTTTTTGCTGCGAAAGCTTTTCCAAAAACAG‑GTT‑AGATATTTTTTTTGCCGTTTATTCCCTAAAAAACATGAA‑TGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                  <  1:784151/150‑1 (MQ=255)
                                                             tCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                          >  2:780930/1‑151 (MQ=255)
                                                               aTGCATTCTTCTTGCGGCG‑AAGGCTTTTACAAAAAAGTGTT‑AAAAAATTTTTTTGACGTTTAGGCCCTAAAAAACGTGTT‑GGAGCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                          <  1:660987/149‑1 (MQ=255)
                                                                tGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                       >  2:231453/1‑151 (MQ=255)
                                                                       ttcttGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  2:132716/1‑151 (MQ=255)
                                                                            gctgcGCAAGCGTTTTCCAGTACAG‑GTT‑AGATGATCTTTTTGTGGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                             >  1:96022/1‑149 (MQ=255)
                                                                              tgcGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                         >  1:359277/1‑151 (MQ=255)
                                                                                   aagcttttcaaaaaaCAG‑TTT‑AAATATTTTTTTTGTCGTTAATTCCCTGAAAAACATTCA‑TGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                     <  1:158349/137‑1 (MQ=255)
                                                                                   aaGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                    >  1:750671/1‑151 (MQ=255)
                                                                                   aaGCGTTTCCAAAAACGG‑GTT‑AAATGATCTTTTTGTCGTTTAATGCCTGAAAAACATGCA‑TGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                      <  2:565673/149‑1 (MQ=255)
                                                                                         tttCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                              >  1:94899/1‑151 (MQ=255)
                                                                                          ttCCAGAACAG‑GTT‑AGATGATCTTTTTGGCGCTTAAAGCCTGTAACACCTGCA‑TAGGCCACCAAATACTCTAAAAAATCCCGCCATTAAGTTAACTGTTAGCGCCCATATCTCCCGAATGCCGCCGTTTGCAACAAATTCGTCGGTAAGc                                                              >  1:529128/1‑150 (MQ=255)
                                                                                            ccAAAACAG‑GTA‑AGTTGATTTTTTGTTCTTTTATGGCCTTAAAAAATTGCA‑GGAACAAAAAAATAAAAAAAAACACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                             <  1:354375/149‑1 (MQ=255)
                                                                                                          aaattatttttttGTCGCTTAATGCTTAAAAAACAGGAA‑TGAGCCAAAAAATAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                                 <  1:106500/141‑1 (MQ=255)
                                                                                                            atgtttttttttGCGCTAAATGCCTTAAAAACAGGTA‑GGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                                 <  2:485951/141‑1 (MQ=255)
                                                                                                             taatttttttGGCGAATATTGCCTTAAAAACAAGCA‑GGACCAAAAAAAAAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                             <  1:708301/145‑1 (MQ=255)
                                                                                                                tttttttGTCATTTATTCCTTGAAAAACATGAA‑TGAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                          <  2:790300/148‑1 (MQ=255)
                                                                                                                tCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGTAGATTTGCATTGATTTACGtc                                         >  2:423739/1‑151 (MQ=255)
                                                                                                                 ttttttGTCGTTTAATCCTTGAAAAACTTGCA‑TGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                         <  2:614755/149‑1 (MQ=255)
                                                                                                                  tttttGTCGTTTAATGCCTGAAAAACATGCA‑TGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                        <  2:776181/150‑1 (MQ=255)
                                                                                                                    tttGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                      >  2:455373/1‑150 (MQ=255)
                                                                                                                      tGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                   >  1:106089/1‑151 (MQ=255)
                                                                                                                                     ttAAAACTTTCATTGACCCAAAAAATAATATAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                      <  1:45468/149‑1 (MQ=255)
                                                                                                                                       aaacaatgct‑tgAGCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                  <  1:11599/146‑1 (MQ=255)
                                                                                                                                          aCAGGCA‑GGACCAACAAAATAATAAAAAAAAACCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                <  1:540004/150‑1 (MQ=255)
                                                                                                                                           catgca‑tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa               >  2:435888/1‑150 (MQ=255)
                                                                                                                                           catgca‑tgAACCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa               <  2:105881/150‑1 (MQ=255)
                                                                                                                                               ca‑tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg           >  1:680790/1‑150 (MQ=255)
                                                                                                                                                a‑tgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga          <  2:297622/150‑1 (MQ=255)
                                                                                                                                                      ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGTTTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga      >  1:511518/1‑149 (MQ=255)
                                                                                                                                                      ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat    >  2:301258/1‑151 (MQ=255)
                                                                                                                                                        acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat  >  2:681467/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAG‑GTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCA‑TGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29336‑29635

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