Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA47.7% 34.8 / 26.9 44intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (18/5);  new base A (0/21);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 3.40e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTG  >  NC_000913/29347‑29641
                                                                                                                                                       |                                                                                                                                                 
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                   >  2:294444/1‑150 (MQ=255)
             tGACCATTTGGTCCTCTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                      <  2:278138/150‑1 (MQ=255)
               aCCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                   >  1:827036/1‑151 (MQ=255)
                 cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                  >  2:332198/1‑150 (MQ=255)
                 cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGACACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                  >  1:183384/1‑150 (MQ=255)
                   tttGGCCCCTTTTTTTCTCCTCTTTTTTATTTCATGCAATCTTTTGGCTGCGAAAGCTTTTCCAAAAAAAGGTTAAATAATTTTTTTGTCGTTTAATGCCTGAAAAACAGGAATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAg                                                                                                                                <  1:333971/150‑1 (MQ=255)
                          cttttttttCTGCTTTTTTTTTTTTCATGCAATCTTTTGGCGGCGCAAGCTTTTCCAAAAACGGTTTAGATAATTTTTTTGTCGCTTAATGCCTGTAAAACATGAAGGAGCCAAAAAAAAATATAAAAAACCCCGCCATTAAGTTGAc                                                                                                                           <  1:408472/145‑1 (MQ=255)
                            cTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                       >  2:381236/1‑150 (MQ=255)
                             tttttcCTGCTCTTTTTTATTCCATGCATTCTTTTTGCTGCGAAAGCTTTTCCAAAAACGGGTAAGATGATCTTTTTGTCGTTTATTGCCTGAAAACCATGCAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                        <  2:229633/148‑1 (MQ=255)
                              tttttcttttgttttatATTTAATCCAATTTTCTTCCTGCGAAGGCTTTTCCAAAAACGGGTTAGATATTTTTTTTGCCGTTAATTCCTTAAAAATCATGAATGGCCCAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAg                                                                                                                     <  1:814948/138‑1 (MQ=255)
                                  tCTGCTAGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCAAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATACTATAAATAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                 >  1:701030/1‑150 (MQ=255)
                                       tctttttttttTCATGCAATTTTCTTGCTGCGAAGGCTTTTCCAAAAACAGTTTAGATAATTTTTTTGTCGTTAAATCCCTGAAAACCATGCATGACCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                              <  2:793981/145‑1 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                          >  1:96398/1‑151 (MQ=255)
                                                  tCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                                >  1:581915/1‑151 (MQ=255)
                                                       cAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                           >  1:167892/1‑151 (MQ=255)
                                                                     cgcAAGCTTTTTCAAAAACAGGTTAGATGATTTTTTTGTCGTTTAATGCTTGAAAAACATGAATGACCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                   <  2:770994/145‑1 (MQ=255)
                                                                       cAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                           >  2:766049/1‑151 (MQ=255)
                                                                        aggctttttCAAGAACAGTTTAGATGATTTTTTTGCCGTTTAATCCTTGAAAAACAGGAATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                <  2:478489/140‑1 (MQ=255)
                                                                              ttcccAAAAAAGGTTAAATGATTTTTTTGTCGTTTAAGCCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                     <  2:405906/147‑1 (MQ=255)
                                                                                       cAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCCCCATTAAGTTGACTTTTAGCGCCCATATCTACAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                             >  1:520838/1‑149 (MQ=255)
                                                                                           ttAGATTATTTTTTTTTCGTTTACTCCTGAAAAAACATGAATGGCCCACAAAATAATATAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                        <  1:133192/150‑1 (MQ=255)
                                                                                             aGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                     >  1:259492/1‑151 (MQ=255)
                                                                                                   tttttttGTCGCTAAATGCCTAAAAAACATGCATGACCAAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                     <  2:642015/143‑1 (MQ=255)
                                                                                                   tttttttGTAGCTTAATCCCTAAAAAACATGCATGACCCAAAAAATAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                  <  1:439147/146‑1 (MQ=255)
                                                                                                   tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                               >  2:104991/1‑151 (MQ=255)
                                                                                                    cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                               >  2:140248/1‑150 (MQ=255)
                                                                                                     tttttGTCGTTAAATGCTTGAAAAACATGAATGCGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                <  2:659803/148‑1 (MQ=255)
                                                                                                        ttGCCGTTTATTGCTTGAAAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                           <  2:784729/150‑1 (MQ=255)
                                                                                                            cccTTAATGCCTAAAAAACAGAAATAACACAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                      <  1:612960/149‑1 (MQ=255)
                                                                                                            cGCTTAATCCCTGAAAAACATGCATGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                        <  2:674125/149‑1 (MQ=255)
                                                                                                             cgctaaacccctaaaaaaCATCCATGACACAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                       <  1:754135/136‑1 (MQ=255)
                                                                                                                  attgcttaaaaaaCATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                 <  2:699209/142‑1 (MQ=255)
                                                                                                                          tAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                          >  2:75682/1‑150 (MQ=255)
                                                                                                                           aaaaCATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                           <  1:349965/148‑1 (MQ=255)
                                                                                                                              aCATGCATGACCCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                      <  1:742940/150‑1 (MQ=255)
                                                                                                                                       gtaacacaAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGATAAGCAGATTTGCATTGATTTACGTAATCATTGTGAATTAATATGCAAATaa                     <  1:833945/139‑1 (MQ=255)
                                                                                                                                           cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat          >  2:270505/1‑150 (MQ=255)
                                                                                                                                           cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat          >  2:119745/1‑150 (MQ=255)
                                                                                                                                           aaaaaaaTAATAAAAAAAAGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATGTACGTCATCATTGTGAAGTAATATGCAAATAAAGTGAGTGAata         <  2:623607/148‑1 (MQ=255)
                                                                                                                                              aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc      >  2:451291/1‑151 (MQ=255)
                                                                                                                                               aaaaaaTAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct     <  1:19767/151‑1 (MQ=255)
                                                                                                                                                 ataataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct   <  1:332161/151‑1 (MQ=255)
                                                                                                                                                  taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg  >  2:202398/1‑151 (MQ=255)
                                                                                                                                                  aaataTAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaaga         <  1:271782/143‑3 (MQ=255)
                                                                                                                                                   aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg  <  2:563815/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTG  >  NC_000913/29347‑29641

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