Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,740,4670TG25.5% 53.0 / 29.2 47intergenic (+32/+39)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (14/21);  new base G (7/5);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGG  >  NC_002947/5740323‑5740609
                                                                                                                                                |                                                                                                                                                
aGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                              <  1:359809/149‑1 (MQ=255)
     cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcg                                                                                                                                         <  7:78451/149‑1 (MQ=255)
         aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGaa                                                                                                                                     <  2:182557/149‑1 (MQ=39)
          tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                    <  2:260898/149‑1 (MQ=39)
            tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACTGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                                  >  4:334727/1‑149 (MQ=38)
            tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACTGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                                  <  3:334727/149‑1 (MQ=38)
            tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                                  <  1:166537/149‑1 (MQ=255)
              aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                     <  4:126414/144‑1 (MQ=255)
              aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                     >  3:126414/1‑144 (MQ=255)
              aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGGTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGTGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                <  5:161290/149‑1 (MQ=255)
              aTGGTGAATGGGCAATGGGTGGGGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                <  8:198197/149‑3 (MQ=39)
               tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTCTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCTTAGCGTGCAAGCCCGCGCAGGGttt                                                                                                                               >  4:111491/1‑148 (MQ=255)
                  tGATTGGGCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAATATTAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                            <  8:234905/149‑1 (MQ=255)
                   gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTCGGTTCCGCCCTTAGAGGGGCAGCCCGCGCAGGGGTgga                                                                                                                             >  7:345517/1‑145 (MQ=255)
                    aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                           <  8:214976/148‑1 (MQ=255)
                    aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                           <  8:134863/148‑1 (MQ=255)
                     aTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCTTCGCGTGCAAGCCCCCCAAGGGTTgcaacgc                                                                                                                         >  8:112059/1‑149 (MQ=255)
                                  ggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGGAGCCCGTAGCGTGGAAGCCCGCGCAGGGTTGGAACGCCAccccggtgcgc                                                                                                            >  1:24257/1‑142 (MQ=255)
                                     cggcCCGGGGCGGCATGCCCAGGAATATGAAAGCGGGCGGGGGTTTGCCATGGTGTTGCGCCTGCTTTTGGGCTGAGCTTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCa                                                                                                         <  8:321366/147‑1 (MQ=255)
                                                   aTGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGGGTTGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTa                                                                                          <  7:112059/149‑1 (MQ=255)
                                                           gAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATAGCCTGGAAGCCACCGAAAGGTTGCAAAGCAACCCCAAAAAACTCAGT‑GGCTTTTTAATAttttgag                                                                                   >  6:361254/1‑144 (MQ=255)
                                                                                tttGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGATGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑g                                                            <  3:111491/149‑1 (MQ=255)
                                                                                     cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATCGCGTGGAAGCCCGCGCAGGGTTGGAACGCAACCCCAAAAAAAAGAGT‑GGAATTTTTATATTTTCTGATAGGACCCCCGCCACATAAA‑GAGGCt                                                       >  6:153406/1‑149 (MQ=255)
                                                                                      aaGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTTCAGCCCCTAGCGGGGCAGCCCCCCGACGGTTGGAAACCCACCCCAAAAAAAAAAGTGGAATTTTTATAATTTG‑TTATAGGCCGCCCCCCACATAAA‑CCggggg                                                      >  4:118866/1‑146 (MQ=255)
                                                                                        ggTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGgt                                                    <  3:118866/149‑1 (MQ=255)
                                                                                        ggTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGgt                                                    <  7:334808/149‑1 (MQ=255)
                                                                                                gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGGAGCCCATAGCGTGGAAACCCCCGCAAGGTTTGAACGCCACCCCCAAAAAAAAAAT‑GGATTTTTTGTATTTTGTGATAGGACGCCCCCCACAACAA‑ACGGGTGgtgtgctaac                                            >  2:87749/1‑143 (MQ=255)
                                                                                                   aGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGc                                                        <  5:347210/134‑1 (MQ=255)
                                                                                                                        gggtTTGCGTTGCTGACCTTAGCGTGCAAGCCCTCGAATGGTTGCAACGCAACCCCAAATAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTTTCATACCCCTGCTGCCGCGCCCTGGAcagca                    <  5:361254/149‑1 (MQ=255)
                                                                                                                         gggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                   <  1:54908/149‑1 (MQ=255)
                                                                                                                          ggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTAAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                  >  4:180077/1‑149 (MQ=255)
                                                                                                                           gttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                  <  5:42582/148‑1 (MQ=255)
                                                                                                                           gttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                  >  6:42582/1‑148 (MQ=255)
                                                                                                                            ttgcgttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                >  4:241596/1‑149 (MQ=37)
                                                                                                                            atgcgttgcTGCCCTTCGGGTGCAAGCCCGCTATGGGTTCCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                <  1:209349/148‑1 (MQ=255)
                                                                                                                                gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctcc            >  6:209828/1‑148 (MQ=38)
                                                                                                                                 ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAACGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa            <  1:87749/149‑1 (MQ=255)
                                                                                                                                 ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                     >  7:323157/1‑139 (MQ=255)
                                                                                                                                 ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                     <  8:323157/139‑1 (MQ=255)
                                                                                                                                   gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc         >  1:282777/1‑149 (MQ=255)
                                                                                                                                    cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAGCTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt        >  4:86861/1‑149 (MQ=255)
                                                                                                                                       cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCctgct                                      >  4:267709/1‑116 (MQ=255)
                                                                                                                                       cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCctgct                                      <  3:267709/116‑1 (MQ=255)
                                                                                                                                       cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc     <  3:114495/149‑1 (MQ=255)
                                                                                                                                        ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc      >  4:73881/1‑147 (MQ=255)
                                                                                                                                         cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg   >  3:197191/1‑149 (MQ=39)
                                                                                                                                          ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCCCGCTCCCAgg  >  6:304280/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGT‑GCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAA‑GCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGG  >  NC_002947/5740323‑5740609

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