Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,740,4670TG21.9% 61.4 / 15.4 43intergenic (+32/+39)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base G (6/3);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG  >  NC_002947/5740320‑5740613
                                                                                                                                                    |                                                                                                                                                   
gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGt                                                                                                                                                     <  3:150668/148‑1 (MQ=255)
gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTg                                                                                                                                                    <  6:95079/149‑1 (MQ=255)
           gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGa                                                                                                                                        <  8:41106/149‑1 (MQ=255)
            aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGGGCTTG‑CACGCGaa                                                                                                                                       <  1:231818/149‑1 (MQ=255)
             tgggATGGGGAATGGGCAAAGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGCTGAAAGCCGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTTGGCTGAGCGTACCGGGGTTGCCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGaa                                                                                                                                       <  1:38495/146‑1 (MQ=255)
               tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAACCCTTCGCGTGCAAG‑CCCGCGAAgg                                                                                                                                     <  2:247343/148‑1 (MQ=255)
                  tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGtt                                                                                                                                  <  2:77629/148‑1 (MQ=255)
                  tGGGGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTGTGCCAAGGTGATGCGCCAGCTTTTGGGCTGAGCGTACCGGGCTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTg                                                                                                                                 <  2:262555/149‑1 (MQ=255)
                   ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAg                                                                                                                                      <  6:168065/143‑1 (MQ=255)
                   ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAg                                                                                                                                      >  5:168065/1‑143 (MQ=255)
                       aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTgcaac                                                                                                                             <  3:162475/148‑1 (MQ=255)
                         tGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTgcaacgca                                                                                                                          <  5:176442/149‑1 (MQ=255)
                          gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGGGctt                                                                                                                                                 >  5:280692/1‑124 (MQ=255)
                          gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGGGctt                                                                                                                                                 <  6:280692/126‑3 (MQ=255)
                                            cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑GCTTTAAGCCCATAGCCTGCAAC‑CCCGCCCAGGGTTGGAACGCCACCCCAAAAAAcaaaatg                                                                                                        >  2:207127/1‑142 (MQ=255)
                                                        gCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTTGAGCCCATCGCGTGCAAG‑CCCGCGCAGGGGTGGAACGCCACCCCAAAAAACTCAGTGGATTTTTACTAtt                                                                                           >  5:92484/1‑149 (MQ=255)
                                                            atgAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAGCCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                        <  1:99757/147‑1 (MQ=255)
                                                            aGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTTCAGCCCATAGCGCCCAAG‑CCCGCGCAGGGTTTCCACGCAACCCCCAAAAAAAAACTGCATTTTTTATAttttga                                                                                       >  2:291886/1‑146 (MQ=255)
                                                            aGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTAACGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAGTCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                        <  1:95076/149‑1 (MQ=255)
                                                              gAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGGTTGAGCACATAGCGGGGCAA‑CCCGCGCAGGGTTTGAAACCAACCCCACAAAAAACAGTGCAGTTTTTATAtttttttt                                                                                     >  8:263537/1‑144 (MQ=255)
                                                                aGATGAAAGCGGGCGGGCGTTTGCCCAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGAt                                                                                    <  8:300767/148‑1 (MQ=255)
                                                                 gATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTTCAGCCCATAGCGTGGAAG‑CCCGCGCAGGGTTGGAACGCCACCCCAAAAAACTCACTGGCTTTTTTCTATTTgggggtcg                                                                                  >  6:331351/1‑142 (MQ=255)
                                                                                   tttGCCAAGGTGTCGCGCCAGCTATTGGGCTGAGCGCACCGCGGTTG‑GGTGCCACCCCTACGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAg                                                                <  7:170267/149‑1 (MQ=255)
                                                                                     tGCCAAGGGGTTGCCCCTGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTGGCGTCCACG‑CCGGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                               <  8:218986/148‑1 (MQ=255)
                                                                                      gCCAAGGGGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgccc                                                                           <  7:194872/135‑1 (MQ=255)
                                                                                       ccAAGGTGTTGCGCAAGCTTTTGGGCTGAGCGTACCGGGGTTT‑CGTTTGAGCCCATAGCGGGCCAG‑CCCCCCAACGGTTGGAAACCAACCCCAAAAAAAACACTGGATTTTTTATATTTTGGGAAAGGCCCCCCCCCAGATAAAACGGc                                                            >  2:70807/1‑149 (MQ=255)
                                                                                           ggTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGCGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                        <  1:291886/149‑1 (MQ=255)
                                                                                                 gcgcCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCCAG‑CCCGCGAAGGGTTTCCACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATAAAGCGGCTGGGGTCATa                                                  >  5:182278/1‑149 (MQ=255)
                                                                                                           ttGGGCTGTGGGTACGGGGGTTG‑GTTTGCAGGCGTTCGCGTGAAAG‑CCGTCTAAGGGTTTAAACGCAACACCAAAAAACTCAGTGGATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgc                                           <  3:336452/146‑1 (MQ=255)
                                                                                                                   aGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGccc                                 <  6:92484/148‑1 (MQ=255)
                                                                                                                   aGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                <  8:111745/149‑1 (MQ=255)
                                                                                                                         ccGGGGTTG‑CGTTGCAGCCCTTCGCGGGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATa                                                  <  7:17372/125‑1 (MQ=255)
                                                                                                                            gggTTG‑CGTTGCAGCCCTTAGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGGTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                       >  7:11442/1‑149 (MQ=255)
                                                                                                                              gttg‑cgttgcAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAATAACTCAGTGCATTTTTACTATTTGCTGTTCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                     >  7:174936/1‑149 (MQ=255)
                                                                                                                              gttg‑cgttgcAACCCTTCGCGGGCTTG‑CACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAAGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                     >  8:165408/1‑149 (MQ=37)
                                                                                                                               ttg‑cgttgcAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTATCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                     >  4:16175/1‑148 (MQ=255)
                                                                                                                               tgg‑cgttgcAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                    <  8:84478/147‑1 (MQ=255)
                                                                                                                                         aacccTTCGCGGGCTTG‑CACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc           >  2:120797/3‑149 (MQ=255)
                                                                                                                                         aGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc           >  4:216126/1‑149 (MQ=255)
                                                                                                                                           cccTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc          >  1:305732/1‑148 (MQ=255)
                                                                                                                                             cTTCGCGGGCTTG‑CACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg       >  3:201428/1‑149 (MQ=39)
                                                                                                                                              ttCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                    <  2:82001/135‑1 (MQ=255)
                                                                                                                                              ttCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                    >  1:82001/1‑135 (MQ=255)
                                                                                                                                                  cgTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg  <  6:182278/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTG‑CGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG  >  NC_002947/5740320‑5740613

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: ‑