Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 68,340 G→T 9.4% intergenic (+28/+147) PP_0057 → / ← PP_0058 major facilitator family transporter/1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294768,3400GT 9.4% 160.0 / 4.4 53intergenic (+28/+147)PP_0057/PP_0058major facilitator family transporter/1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base T (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TATACCGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAG  >  NC_002947/68196‑68485
                                                                                                                                                |                                                                                                                                                 
tataCCGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTaaaa                                                                                                                                              <  1:136711/150‑1 (MQ=255)
       gggATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGAAATAAAAAAACCGGCTTTTACAGCCGGtt                                                                                                                                       <  2:71473/150‑1 (MQ=255)
        ggATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGTCATAAAAAAACCGGCTTTTACAGCCGGttt                                                                                                                                      <  1:822207/150‑1 (MQ=255)
           tatCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGtttttt                                                                                                                                   <  1:1002771/150‑1 (MQ=255)
              cttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTAt                                                                                                                                <  1:404229/150‑1 (MQ=255)
               ttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGtt                                                                                                                                       <  1:833012/142‑1 (MQ=255)
               ttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGtt                                                                                                                                       >  2:833012/1‑142 (MQ=255)
               ttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATg                                                                                                                               >  2:771849/1‑150 (MQ=255)
                tctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGc                                                                                                                              <  1:771849/150‑1 (MQ=255)
                    cGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCa                                                                                                                          >  2:1354922/1‑150 (MQ=255)
                          gtgGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCaaaaaaa                                                                                                                    <  1:1354922/150‑1 (MQ=255)
                                            gACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAg                                                                                                  <  2:1260463/150‑1 (MQ=255)
                                                 ggATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGaaaaa                                                                                             >  1:428269/1‑150 (MQ=255)
                                                       gCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCAt                                                                                       >  2:1178982/1‑150 (MQ=255)
                                                            gtggtTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTAt                                                                                                                                <  2:1226825/104‑1 (MQ=255)
                                                            gtggtTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTAt                                                                                                                                >  1:1226825/1‑104 (MQ=255)
                                                             tggtTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATGCATACAAAg                                                                                 >  1:292443/1‑150 (MQ=255)
                                                                ttGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGaaa                                                                              <  2:292443/150‑1 (MQ=255)
                                                                  gggATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGaa                                                                                                >  2:348986/1‑130 (MQ=255)
                                                                  gggATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGaa                                                                                                <  1:348986/130‑1 (MQ=255)
                                                                   ggATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAAt                                                                                           <  1:270785/134‑1 (MQ=255)
                                                                   ggATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAAt                                                                                           >  2:270785/1‑134 (MQ=255)
                                                                                  gagaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTATTTCATCGCACAAAAAAACTTATGCACGATTTTCAAACACCATTAATATAAACAAACAACTAGCAAATTCATac                                                            >  1:385973/1‑150 (MQ=255)
                                                                                   agaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAAt                                                                             <  1:1399006/132‑1 (MQ=255)
                                                                                   agaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAAt                                                                             >  2:1399006/1‑132 (MQ=255)
                                                                                   agaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTAATACAAAGAAATAACTAGCAAATTCATaca                                                           >  1:10382/1‑150 (MQ=255)
                                                                                        ccAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAg                                                      <  2:428269/150‑1 (MQ=255)
                                                                                         cAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGtt                                                                                                                                       <  1:1001458/68‑1 (MQ=255)
                                                                                         cAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGtt                                                                                                                                       >  2:1001458/1‑68 (MQ=255)
                                                                                          aaGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGat                                                    >  1:1365591/1‑150 (MQ=255)
                                                                                                    atatCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCaa                                          <  1:1178982/150‑1 (MQ=255)
                                                                                                    atatCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCaa                                          <  1:1271467/150‑1 (MQ=255)
                                                                                                    atatCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCaa                                          <  1:1206980/150‑1 (MQ=255)
                                                                                                     tatCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATAAACTTAGATATCAATTCAAg                                         >  2:635582/1‑150 (MQ=255)
                                                                                                         gACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCAtt                                     >  2:319844/1‑150 (MQ=255)
                                                                                                              aGTCTGATGTCCACATGCCATAAAACTCACGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATc                                <  1:361239/150‑1 (MQ=255)
                                                                                                                    aTGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTAAAGCATTTAATCCaaaaa                          >  1:63916/1‑150 (MQ=255)
                                                                                                                        ccACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGc                      >  2:154468/1‑150 (MQ=255)
                                                                                                                              gCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACaa                <  2:1365591/150‑1 (MQ=255)
                                                                                                                                  tAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTa                                   >  1:228168/1‑127 (MQ=255)
                                                                                                                                  tAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTa                                   <  2:228168/127‑1 (MQ=255)
                                                                                                                                  tAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACaa                >  2:1207416/1‑146 (MQ=255)
                                                                                                                                  tAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACaa                <  1:1207416/146‑1 (MQ=255)
                                                                                                                                     aaaaaCCGGCTTTTACAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATcc         >  2:957661/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGtt             >  2:1218260/1‑146 (MQ=255)
                                                                                                                                     aaaaaCCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGtt             <  1:1218260/146‑1 (MQ=255)
                                                                                                                                          ccGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTTATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGAACACAAGTTATCCACAGa    >  1:192203/1‑150 (MQ=255)
                                                                                                                                          ccGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAAAGTATGTGCGCTTCTCAGAAAAAGTTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATACACAGa    >  1:389240/1‑150 (MQ=255)
                                                                                                                                            ggCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATa                                                             <  1:780683/91‑1 (MQ=255)
                                                                                                                                            ggCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATa                                                             >  2:780683/1‑91 (MQ=255)
                                                                                                                                            ggCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAg  >  1:657328/1‑150 (MQ=255)
                                                                                                                                               tGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATc          >  1:80896/1‑139 (MQ=255)
                                                                                                                                               tGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATc          <  2:80896/139‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TATACCGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAG  >  NC_002947/68196‑68485

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