Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,637,1330GC49.0% ‑2.8 / 70.1 52A277G (GCC→GGC) ssuB_2Aliphatic sulfonates import ATP‑binding proteinSsuB
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base C (14/11);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCAGTTGCCAGACGCCGAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGACGCGCGCCCATCACGA  >  ACB122/3637002‑3637268
                                                                                                                                   |                                                                                                                                       
tgcCAGTTGCCAGACGCCGAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTgg                                                                                                                                    >  8:319932/1‑137 (MQ=255)
 gcCAGTTGCCAGACGCCGAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTggg                                                                                                                                   <  5:148807/137‑1 (MQ=255)
      ttGCCAGACGCCGAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAgg                                                                                                                              <  1:106880/137‑1 (MQ=255)
      ttGCCAGACGCCGAGCAGCAAGGAGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAgg                                                                                                                              <  7:319932/137‑1 (MQ=255)
         ccAGACGCCGAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCAATGGGCCAGGAGt                                                                                                                           >  7:162488/1‑137 (MQ=255)
           aGACGCCTAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGATGGGTCGCCATTGGCCTGGGCCAGGAGTCg                                                                                                                         <  6:102753/137‑1 (MQ=255)
           aGACGCCGAGCAGCAAGTCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGTAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCg                                                                                                                         <  2:332457/137‑1 (MQ=255)
                    gcagcaAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTg                                                                                                                >  6:182537/1‑137 (MQ=255)
                       ccaAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTgcgc                                                                                                             <  5:58457/136‑1 (MQ=255)
                                   ccAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGAt                                                                                                 <  7:332211/137‑1 (MQ=255)
                                          tggtgAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGcc                                                                                          >  8:195431/1‑137 (MQ=255)
                                             tgAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGATAGGTGCGCGCATTGCTCGATGAACGCCGGt                                                                                       >  6:141947/1‑137 (MQ=255)
                                              taatcGATGCTTATTACAGATGGATGGCGTAGCCTGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGGCAGGAGTCGGGACAGGTGCGCGCATTGATCGATGAACGCCGGTg                                                                                      <  2:238800/133‑1 (MQ=255)
                                                               aGACGGAAGGCGTCGCCAGGGAGTAGAGGGTCGGCGCATGATGAGGCGCGCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTACGGTACGCCtcg                                                                     <  2:299425/137‑1 (MQ=255)
                                                               aGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCtcg                                                                     >  6:326421/1‑137 (MQ=255)
                                                                        gCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTcc                                                            >  7:77876/1‑137 (MQ=255)
                                                                         cGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCTTTGGCCTggg                                                                                                                                   <  8:121654/65‑1 (MQ=255)
                                                                         cGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTggg                                                                                                                                   >  7:121654/1‑65 (MQ=255)
                                                                         cGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCa                                                           >  7:17185/1‑137 (MQ=255)
                                                                         cGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCa                                                           >  7:75426/1‑137 (MQ=255)
                                                                         cGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCa                                                           >  5:147939/1‑137 (MQ=255)
                                                                           tAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGt                                                         <  5:182537/137‑1 (MQ=255)
                                                                                aGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGcc                                                    >  6:271507/1‑137 (MQ=255)
                                                                                aGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTTGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGcc                                                    <  3:10785/137‑1 (MQ=255)
                                                                                 gggAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGAGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCg                                                   <  2:185561/137‑1 (MQ=255)
                                                                                   gAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGATAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTc                                                 <  5:141947/137‑1 (MQ=255)
                                                                                   gAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGGCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTc                                                 <  7:141763/137‑1 (MQ=255)
                                                                                     gTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGCCCGGTGCGCGCATTGCTCGATGCACGCCGGTGAAGTGCGGTACGCCTAGTCGCGTTCCAGTGGGgcggcgg                                               >  2:352555/1‑136 (MQ=255)
                                                                                                 cTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGGTTATTAAATCgg                                   >  1:139816/1‑137 (MQ=255)
                                                                                                 cTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCgg                                   >  2:201472/1‑137 (MQ=255)
                                                                                                 cTCATGATGATGCTCTCAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGCAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCTGCGCTGATTACATCgg                                   >  4:3509/1‑137 (MQ=255)
                                                                                                   cATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATAGGCg                                 >  6:6265/1‑137 (MQ=255)
                                                                                                          atgCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAgg                                                                                                                              <  1:65950/37‑1 (MQ=255)
                                                                                                          atgCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAgg                                                                                                                              >  2:65950/1‑37 (MQ=255)
                                                                                                              tctcAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGAcc                      >  5:314544/1‑137 (MQ=255)
                                                                                                              tctcAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCtcg                                                                     >  3:66627/1‑90 (MQ=255)
                                                                                                              tctcAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCtcg                                                                     <  4:66627/90‑1 (MQ=255)
                                                                                                              tctcAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGa                                >  1:253241/1‑127 (MQ=255)
                                                                                                              tctcAGGAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGa                                <  2:253241/127‑1 (MQ=255)
                                                                                                              tctcAGGAAGGGTCGCCATTGACCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTAAATCGGCCATCACTCGAcc                      >  2:64930/1‑137 (MQ=255)
                                                                                                                   ggAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGAc                 <  4:157580/137‑1 (MQ=255)
                                                                                                                   ggAAGGGTCGCCATTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGAc                 >  3:157580/1‑137 (MQ=255)
                                                                                                                     aaGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGCcgc               >  3:42396/1‑137 (MQ=255)
                                                                                                                     aaGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGAcgc               <  7:195431/137‑1 (MQ=255)
                                                                                                                     aaGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGAcgc               >  6:109324/1‑137 (MQ=255)
                                                                                                                             ccATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGAcgc               >  4:287965/1‑129 (MQ=255)
                                                                                                                             ccATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGAcgc               <  3:287965/129‑1 (MQ=255)
                                                                                                                               aTTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATc                                     <  7:113135/105‑1 (MQ=255)
                                                                                                                               aTTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATc                                     >  8:113135/1‑105 (MQ=255)
                                                                                                                               aTTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATc                              <  8:57425/112‑1 (MQ=255)
                                                                                                                               aTTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATc                              >  7:57425/1‑112 (MQ=255)
                                                                                                                               aTTGCCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGACGCGCGCCCATCa     >  5:285248/1‑137 (MQ=255)
                                                                                                                                  ggCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGACGCGCGCCCATCACGa  <  1:201472/137‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TGCCAGTTGCCAGACGCCGAGCAGCAAGGCGCCAACCAGTAGTGGTGAAACGATGCGTAGTACAGAAGGATGGCGTAGCCAGGGAGTAGAGGGTCGGCTCATGATGATGCTCTCAGGAAGGGTCGCCATTGGCCTGGGCCAGGAGTCGGGACAGGTGCGCGCATTGCTCGATGAACGCCGGTGAAGTGCGGTACGCCTCGTCGCGTTCCAGTGGGCCGTCGATGATTACATCGGCGATCACTCGACCGGGACGCGCGCCCATCACGA  >  ACB122/3637002‑3637268

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