Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,973,7770AC36.5% 61.3 / 13.9 53L161R (CTG→CGG) cobNAerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base A (13/20);  new base C (19/0);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCG  >  ACB143/3973648‑3973903
                                                                                                                                 |                                                                                                                               
agcGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGGACCAGTCCTCGAGTGGCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGAGCTGCTAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCAtt                                                                                                                           <  4:182632/136‑1 (MQ=255)
  cGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCCTtgg                                                                                                                         >  3:349211/1‑136 (MQ=255)
  cGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCCTtgg                                                                                                                         >  4:12085/1‑136 (MQ=255)
  cGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATtgg                                                                                                                         >  2:63342/1‑136 (MQ=255)
    gTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGGCCCAGGCATAGTCGCGGCCCCCCCCTtggcc                                                                                                                       >  2:327325/1‑135 (MQ=255)
      gCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATtggcg‑gg                                                                                                                    >  7:257173/1‑136 (MQ=255)
        cACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGcc                                                                                                                  <  1:53218/136‑1 (MQ=255)
        cACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGcc                                                                                                                  <  2:13010/136‑1 (MQ=255)
          cTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGCG‑GGCCAg                                                                                                                >  7:256003/1‑136 (MQ=255)
                 ctCGGGGTGCCAGTGTGGGTACCAGTCCTCTAGTGTCGCGCTGCCCTTGGCTGTGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGtt                                                                                                          <  5:247569/135‑1 (MQ=255)
                 ctCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGtt                                                                                                          <  7:248581/135‑1 (MQ=255)
                    ggggTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGaa                                                                                                       >  2:34690/1‑135 (MQ=255)
                    ggggTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCAGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCCCCCATGGGG‑GGGCCGGCAGTTGaa                                                                                                       >  2:102464/1‑135 (MQ=255)
                              tgtgGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCCTTGGCT‑GGCCAGGCAGGTGAACAGGTTGATgg                                                                                            >  4:283960/1‑136 (MQ=255)
                                   taaacAGTCCTCGAGAGTCGCGCTGACCTTATCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCCCCATTGGCT‑GGCCAGGCACTTGAACAGGTTTAAGGCGTTg                                                                                       >  8:55364/5‑136 (MQ=255)
                                   gTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCTGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTg                                                                                       >  2:227993/1‑136 (MQ=255)
                                   gTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCCCCCTTGGGT‑GGGCAGGGACTTTAACAGGTTGAAGGCGttt                                                                                       >  4:349244/1‑135 (MQ=255)
                                   gTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTg                                                                                       >  1:16781/1‑136 (MQ=255)
                                   gTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCAATTGGT‑GGGCAGGGAGGTGAAAAGGTTTATTGCGttt                                                                                       >  3:253644/1‑135 (MQ=255)
                                     aCCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCCCCCATTGCG‑GGGCAGGCAGTTGAACAGGTTGATGGGGTTgtc                                                                                     >  5:108376/1‑134 (MQ=255)
                                      ccAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGcc                                                                                    <  4:25990/136‑1 (MQ=255)
                                       cAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCCCCATTTGCT‑GGCCAGGCAGTTGAACACGTTGATGGCGTTGGCCt                                                                                   >  7:248637/1‑136 (MQ=255)
                                            cTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCAGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTgcc                                                                               <  1:102464/135‑1 (MQ=255)
                                                 gtgtCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTg                                                                         <  6:35691/136‑1 (MQ=255)
                                                      gcgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCa                                                                    <  1:34690/136‑1 (MQ=255)
                                                      gcgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCa                                                                    <  6:323517/136‑1 (MQ=255)
                                                       cgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGGCCAACCAATTGCTGGGCCCGGCCG‑TGGACAGGGTGTTTGTGTTTGGCTTTGCGCCCCGgcgccg                                                                   >  4:132564/1‑134 (MQ=255)
                                                       cgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATTGCGTTGGGCTTGCCGCCCTGGCGCAg                                                                   >  1:282420/1‑136 (MQ=255)
                                                       cgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAg                                                                   >  7:100421/1‑136 (MQ=255)
                                                       cgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGCT‑GGCCAGGCAGTTGAACACGTTTATGGCGTTGGCCTTTCCGCCCTGGCGCAg                                                                   >  5:135885/1‑136 (MQ=255)
                                                              cTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGc                                                            <  3:132564/136‑1 (MQ=255)
                                                                    ggggTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGc                                                      <  8:168379/136‑1 (MQ=255)
                                                                         gATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGCG‑GGGCAGGCAGTTGAAAAGGTTGATTGCGTTGGCCTTGCCGCCCGGGCGCAGGTAATTCCAAACCCGCTc                                                 >  4:163451/1‑136 (MQ=255)
                                                                         gATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTc                                                 >  8:185851/1‑136 (MQ=255)
                                                                         gATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCACGCATTGGCA‑GGGCAGGGAGTTTGAAAGGTTGATTGCGTTTGCCATGCCGCCCTGGCGCAGGGAATTCCACACCCGccg                                                 >  7:234894/1‑134 (MQ=255)
                                                                              aCCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCCCCCTTTGGT‑GGGCCGGCAGGTGTACAGGTTGATGGCGTTGGGCGTTCCCCCCCGGGGCCGGGAATGCCCAACCCGCGCGGccg                                            >  5:199790/1‑135 (MQ=255)
                                                                              aCCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTTGCT‑GGGCAGGCCGTTGTAAAGGTTGATTGCGTTTGCCATTCCGCCCAGGCGCAGGGAATTCCAAAGCCGCTAGGccg                                            >  3:351023/1‑135 (MQ=255)
                                                                                   cgTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCCTTGGCG‑GGCCAGGCAGTTGTACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTTATGCCAAAGCCGCTCGGCCTGCTCa                                       >  5:76633/1‑136 (MQ=255)
                                                                                    gTGCGCGGCAAAGGCTGCGGCGCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCAc                                      <  6:18850/136‑1 (MQ=255)
                                                                                     tGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGCT‑GGGCCGGCAGTTGAACAGGTTGATGGGGTTGGCCTTGCCGCCCTGGCGCCGGGAATGCCAAAGCCGCTCGGCCGGCTCAcc                                     >  3:317020/1‑136 (MQ=255)
                                                                                     tGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCCTTGGGG‑GGGCAGGGAGGTGATCAGGGTGATGGGGTTGGCCTTGCCGCCCGGGCGCAGGTAATGCCCAAGCCGCGCGGCCGGCTCGcc                                     >  3:325290/1‑136 (MQ=255)
                                                                                        gcgGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACcgcg                                  <  2:16781/136‑1 (MQ=255)
                                                                                         cgGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCgg                                 <  2:33196/136‑1 (MQ=255)
                                                                                         cgGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCgg                                 <  4:70434/136‑1 (MQ=255)
                                                                                             aaaGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCCGCCATTGGGT‑GGGCAGGGCGTTGAACAGGTTTATGGCGTTGGGCTTGCCGCCCTGGCGCAGGGAATGCCAAAGCCGCGCGGCCTGCTCCCCGCGGCCgg                             >  4:308275/1‑135 (MQ=255)
                                                                                                       gCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATTGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGcg                   >  5:255079/1‑135 (MQ=255)
                                                                                                       gCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCt                   <  3:131824/136‑1 (MQ=255)
                                                                                                        cTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTg                  <  2:27217/136‑1 (MQ=255)
                                                                                                         tCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCggg                 >  2:31723/1‑136 (MQ=255)
                                                                                                                    aTAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGGAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGCCGTTGCCCAGGCGGGTCAGCtccggg      >  1:251982/1‑136 (MQ=255)
                                                                                                                       ggcgcgGCCCAGCCAATGGCT‑GGCCAGGCAGTTGAACAGGGTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCtccgggtcc   <  4:234981/134‑1 (MQ=255)
                                                                                                                       gTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCtccgggtcc   <  7:185851/136‑1 (MQ=255)
                                                                                                                        tCGCGGCCCCGCCATTGGCT‑GGGCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCCGGCTGGGCAGCtccgggtccg  >  3:348378/1‑136 (MQ=255)
                                                                                                                                 |                                                                                                                               
AGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCT‑GGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCG  >  ACB143/3973648‑3973903

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