Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793901,7170A.71.2% 72.4 / 70.8 58coding (437/465 nt)SAUSA300_RS04425iron‑sulfur cluster assembly scaffold protein NifU
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base . (37/0);  minor base C (15/0);  total (58/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GAAGCATTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAAAGCAGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATATTAGAC  >  NC_007793/901620‑901780
                                                                                                 |                                                               
taaGCATTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAg                                                               >  1:833951/2‑100 (MQ=255)
caaGCATTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAg                                                               >  1:2981130/2‑100 (MQ=255)
 aaGCATTGCAAGGTGTATCTCAATTCCTAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAga                                                              >  1:4319455/1‑99 (MQ=255)
 aaGCATTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGATACAga                                                              >  1:3611602/1‑99 (MQ=255)
   gCATTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatc                                                            >  1:1123719/1‑97 (MQ=255)
     gctGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGATGGCaagatcgga                                                          >  1:3111972/3‑93 (MQ=255)
     cttGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAGGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcgg                                                          >  1:709044/2‑95 (MQ=255)
     aTTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAGGGTACAgatcgg                                                          >  1:1585831/1‑95 (MQ=255)
      ttGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGATGGTACAgatcgga                                                         >  1:984275/1‑94 (MQ=255)
      ttGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAGATCGGAAg                                                         >  1:4565587/1‑100 (MQ=255)
      ttGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAGATCGGAAg                                                         >  1:4331167/1‑100 (MQ=255)
      ttGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcgga                                                         >  1:1063501/1‑94 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAGGGTaagatcggaag                                                        >  1:1960912/1‑91 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTaagatcggaag                                                        >  1:4670555/1‑91 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAGATCGGAAga                                                        >  1:1068379/1‑99 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcggaa                                                        >  1:4729983/1‑93 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcggaa                                                        >  1:212787/1‑93 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcggaa                                                        >  1:1511168/1‑93 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcggaa                                                        >  1:2715010/1‑93 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcggaa                                                        >  1:16628/1‑93 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTACAgatcggaa                                                        >  1:5631136/1‑93 (MQ=255)
       tGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGATAGATCGGAAga                                                        >  1:400642/1‑99 (MQ=255)
        gCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAGATCGGAAgag                                                       >  1:4720101/1‑98 (MQ=255)
                          tccAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAAAGCAGAAGGTACGACTGAAGAagaatag                                     >  1:4669177/2‑100 (MQ=255)
                                          tgtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGat         >  1:2094691/1‑100 (MQ=255)
                                          tgtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGat         >  1:125652/1‑100 (MQ=255)
                                          tgtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGat         >  1:2347255/1‑100 (MQ=255)
                                          tgtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGGACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGat         >  1:3943617/1‑100 (MQ=255)
                                           ttgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:3193087/2‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGTAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:4593809/1‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:3829714/1‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:5336184/1‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:2533951/1‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:2224275/1‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:4115808/1‑100 (MQ=255)
                                           gtgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTCGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:5198238/1‑100 (MQ=255)
                                           ctgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:5063897/2‑100 (MQ=255)
                                           ctgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:3422969/2‑100 (MQ=255)
                                           atgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGata        >  1:3856274/2‑100 (MQ=255)
                                            tgCCACATTAGCTTGGGAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGatat       >  1:4732189/1‑100 (MQ=255)
                                            tgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGatat       >  1:2586311/1‑100 (MQ=255)
                                            tgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGatat       >  1:3885281/1‑100 (MQ=255)
                                            tgCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGatat       >  1:3874719/1‑100 (MQ=255)
                                            agCCACATTAGCTTGGAAAACATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGatat       >  1:3640132/2‑100 (MQ=255)
                                             tccACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:1777770/2‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:1001270/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:841473/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:124446/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:676206/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:545856/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:4232835/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:3364098/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:1834256/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:2798754/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:2726282/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:2523474/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:2227677/1‑100 (MQ=255)
                                             gCCACATTAGCTTGGAAAGCATTGGAAAAAGGCACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:4725057/1‑100 (MQ=255)
                                             atcacaTTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATAtt      >  1:2119662/3‑100 (MQ=255)
                                              ccACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATATTa     >  1:4689947/1‑100 (MQ=255)
                                                 caTTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAA‑‑‑‑‑‑‑‑‑‑‑‑AGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATATTAGAc  >  1:387754/1‑100 (MQ=255)
                                                                                                 |                                                               
GAAGCATTGCAAGGTGTATCTCAATTCCCAGCTCGTATTAAATGTGCCACATTAGCTTGGAAAGCATTGGAAAAAGGTACTGTTGCTAAAGAAGGTAAAGCAGAAGGTACGACTGAAGAAGAATAGAATGCTGTTAATCATAGATAATTTTGATATTAGAC  >  NC_007793/901620‑901780

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