Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA47.7% 36.4 / 23.1 44E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (7/16);  new base A (21/0);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508516‑3508811
                                                                                                                                                   |                                                                                                                                                      
gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATca                                                                                                                                                     >  1:229057/1‑150 (MQ=255)
      tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAAttaaaca                                                                                                                                                  >  1:27346/1‑143 (MQ=255)
       ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑GGAAAAAATTacacac                                                                                                                                                >  1:374912/1‑148 (MQ=255)
       ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGGAAAAATcaaa                                                                                                                                                   >  2:63729/1‑143 (MQ=255)
       ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATcaca                                                                                                                                                   >  2:348180/1‑145 (MQ=255)
              cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGC‑GGGAAAAATAACAAACAGTAGccccc                                                                                                                                      >  1:649220/1‑151 (MQ=255)
              cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTCAAACACGGGAGccccc                                                                                                                                      >  2:428517/1‑151 (MQ=255)
                gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACAAACgggaaacccccc                                                                                                                                    >  2:493224/1‑141 (MQ=255)
                  gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACGGGAGCCCCcac                                                                                                                                    >  1:179410/1‑149 (MQ=255)
                  gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACACCGGAGccccc                                                                                                                                      >  1:505265/1‑147 (MQ=255)
                       cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGTAGccccc                                                                                                                                      >  2:271606/1‑142 (MQ=255)
                                   gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                                  <  2:381998/150‑1 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGgcg                                                                                                                  >  1:555996/1‑147 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGAAAAAATCACACACGGTAGCCCCCCCAACAGAAAGCGGGGTGCGgac                                                                                                               >  2:298144/1‑150 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCCG‑GGGGAAAAACACAAACCGTAAACCCCCCAACCAAAAGCggggtggg                                                                                                                  >  2:552505/1‑145 (MQ=255)
                                        tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGca                                                                                                                <  2:264576/147‑1 (MQ=255)
                                         gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACGATTTTGTGGCGGTGGCGCGC‑GGGGAAAATAACAAACAGGAGACCCCACATCAAAAAACGTGGGGCGCacac                                                                                                             >  2:221789/1‑149 (MQ=255)
                                         gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACAcc                                                                                                            <  1:341854/150‑1 (MQ=255)
                                              aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                      <  1:681250/151‑1 (MQ=255)
                                               gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGTAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                     <  1:43506/151‑1 (MQ=255)
                                               gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                     <  1:290255/151‑1 (MQ=255)
                                                ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATTACACACCGTAGCCCCCACCCCACACAGCGCGGTGGGCCCAccc                                                                                                           >  2:591227/1‑143 (MQ=255)
                                                   gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTACACACGGGAGCCCCCACATCAGAAAGAGTGGGGCGCACACCCAGCTGc                                                                                                     >  1:604302/1‑147 (MQ=255)
                                                        ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACAACCGGTAGCCCCCCCCACAGAAAAAGTGGTGGGCACACCGAGCCGCGTTGCg                                                                                               >  1:278558/1‑148 (MQ=255)
                                                              gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCt                                                                                     <  2:336232/151‑1 (MQ=255)
                                                                 cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTg                                                                                   <  1:422568/150‑1 (MQ=255)
                                                                    aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGa                                                                               <  1:246345/151‑1 (MQ=255)
                                                                      cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGCGC‑GGGAAAAATCAAACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCCAGATGCGGAGCGATCC‑TCCCCTGCTGa                                                                               >  2:483668/1‑149 (MQ=255)
                                                                       tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACGGTAGACCCCACATCAGAAAGCGTGGGGCGCACACCGAGCCGCGCAGCGATCA‑TCCCTTTCTGAt                                                                              >  1:571990/1‑149 (MQ=255)
                                                                        cATAGGGTGTCTCCGGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTCGCCCCCTCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTtata                                                                             <  2:616771/149‑1 (MQ=255)
                                                                           aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCa                                                                        <  2:658967/151‑1 (MQ=255)
                                                                                 tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCCAGCGGCGTAGCGATCA‑TCCCCTGCTGATAAAACAGCACTg                                                                  >  1:146898/1‑151 (MQ=255)
                                                                                    ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTg                                                                  <  1:493411/148‑1 (MQ=255)
                                                                                    ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGca                                                               <  2:27341/151‑1 (MQ=255)
                                                                                     cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAATTCACACACGGTATCCCCCCCATCAGAAAGGGTGGGGCGCGCACCGAGCTGCGGAGCGATCA‑TCCCTTGCTGATAAACCCACACTgg                                                                 >  1:149335/1‑148 (MQ=255)
                                                                                         aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAAACACCGGAGCCCCCCCAACAGAAAACGGGGGGGGCACACCCAGCTGCGGAGCGATCC‑TCCCTTGCTGGTAAAACAGCACCGGcacc                                                             >  1:487260/1‑148 (MQ=255)
                                                                                              ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                     <  1:154008/151‑1 (MQ=255)
                                                                                                        aagaaaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGAAAAATAAACACACGGGAGACCCCACAATAGAAAGAGTGGGGGGCCCACACAGCTGCGTTGGGGTCACTCCCTT‑CTGGTTAAAAAGCACTCGCCACACTTCGCgggggggggg                                            >  2:599783/1‑143 (MQ=255)
                                                                                                                   gACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGCCCACCCCGCTGCGTAGCGATCA‑TTCCCTTCTGATAAAACCGCAACGGGACCACCTCTCGGGGATGGTGGCTGTGGccca                                 >  2:490345/1‑148 (MQ=255)
                                                                                                                      cATTTTGTGGCGGGGGCGCGC‑GGAAAAATTTACACACGGGAACCCCCCCATCAGAAAGCGGGGGGCGCACACCGCGCTGCGGAGCGATCA‑TCCCCTGCTGATAAACCAGCACTGGCACCACCTCGCGGGGATGGTGGCGGTGGCCAAAgg                              >  2:238971/1‑150 (MQ=255)
                                                                                                                                    ggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt               <  1:348287/151‑1 (MQ=255)
                                                                                                                                         gc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc          >  2:589518/1‑151 (MQ=255)
                                                                                                                                              aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc       >  2:427050/1‑150 (MQ=255)
                                                                                                                                              aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCa      <  1:66852/151‑1 (MQ=255)
                                                                                                                                                  aTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  >  1:409149/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508516‑3508811

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