Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 516,305 T→C 100% T169T (ACT→ACC bioB → biotin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE516,3050TC81.7% 109.8 / 13.6 71T169T (ACT→ACCbioBbiotin synthase
Reads supporting (aligned to +/- strand):  ref base T (12/1);  new base C (57/0);  total (70/1)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAACCTGGACACCTCGCCGGAGTTTTACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGG  >  minE/516254‑516367
                                                   |                                                              
caaccTGGACACCTCGCCGGAGTTTTACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGa                                               <  1:2792318/69‑1 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGaaaaa                                  >  1:986053/1‑68 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGaaaaa                                  >  1:3303099/1‑68 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGaaaaa                                  >  1:3358487/1‑68 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:1291676/1‑69 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:49690/1‑69 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:525544/1‑69 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:2123235/1‑69 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:2522401/1‑69 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:2618858/1‑69 (MQ=255)
              cGCCGGAGTTTTACGGCAATATCATCACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAg                                 >  1:2841410/1‑69 (MQ=255)
                        ttACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGc                       >  1:1574158/1‑69 (MQ=255)
                        ttACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGc                       >  1:2981990/1‑69 (MQ=255)
                        ttACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGc                       >  1:3091323/1‑69 (MQ=255)
                        ttACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGc                       >  1:3118471/1‑69 (MQ=255)
                        ttACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGc                       >  1:2116112/1‑69 (MQ=255)
                             gCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCggg                   >  1:1026760/1‑68 (MQ=255)
                             gCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCggg                   >  1:2725367/1‑68 (MQ=255)
                             gCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCggg                   >  1:1859776/1‑68 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCgaga                          >  1:809136/1‑56 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:952745/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:118431/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:2502942/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:2873515/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:1311774/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:1290234/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:867352/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:158393/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:728250/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:186274/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:1829198/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaaa             >  1:179740/1‑69 (MQ=255)
                                  atcatcACCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTaa              >  1:220578/1‑68 (MQ=255)
                                  atcatcACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCaaa             >  1:771058/1‑69 (MQ=255)
                                  atcatcACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCaaa             >  1:2316898/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGctg                     >  1:335652/1‑57 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAg            >  1:924879/1‑68 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:518803/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2911560/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2875191/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:3226949/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2559395/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2391428/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:1824573/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:177456/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2763728/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2284590/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGt           >  1:2410389/1‑69 (MQ=255)
                                    catcaCCACCCGCACCTATCAGGAACGTCATGATACGCTGGAAAAAGTGCGAGAAGCTGGGa                  >  1:638483/1‑62 (MQ=255)
                                    catcaCCACCCGCACATATCAGGAACGTCTTGATACGCTGGAAAAAGTGAGAGAAGCTGGGAtt                >  1:1672651/1‑63 (MQ=255)
                                       caccacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTg        >  1:2541013/1‑69 (MQ=255)
                                       caccacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTg        >  1:1258678/1‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:835675/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:83141/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:487252/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:1580494/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:1428991/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:729876/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:133067/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:1153245/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:1871857/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2108761/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2141110/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2171489/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:3267887/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2394844/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2480938/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2911204/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2611706/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2690297/2‑69 (MQ=255)
                                             ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg  >  1:2740157/2‑69 (MQ=255)
                                                   |                                                              
CAACCTGGACACCTCGCCGGAGTTTTACGGCAATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGG  >  minE/516254‑516367

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