Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 543,160 C→T 7.5% T150I (ACA→ATA)  mazG → possible tetrapyrrole methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730543,1600CT 7.5% 156.5 / 2.6 54T150I (ACA→ATA) mazGpossible tetrapyrrole methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/19);  new base T (2/2);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

AGCAAAGGACAATACTGATATAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAGCGATGGTG  >  CP000730/543026‑543297
                                                                                                                                      |                                                                                                                                         
aGCAAAGGACAATACTGATATAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTaa                                                                                                                                     >  2:427928/1‑141 (MQ=255)
      ggACAATACTGATATAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAg                                                                                                                               >  2:280827/1‑141 (MQ=255)
            tACTGATATAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAAAACTTAATGTTAGAACGCa                                                                                                                         >  2:160418/1‑141 (MQ=255)
                gATATAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACa                                                                                                                     >  2:334331/1‑141 (MQ=255)
                 atatAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACAt                                                                                                                    <  1:433788/141‑1 (MQ=255)
                   ataGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTg                                                                                                                  >  2:17725/1‑141 (MQ=255)
                   ataGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTg                                                                                                                  <  1:423758/141‑1 (MQ=255)
                      gATGTGAAAGTTTTAGGTGGTAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAATACTTAATGTTAGAACGCATACATTGAtt                                                                                                               >  1:183653/1‑141 (MQ=255)
                          tgAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGc                                                                                                           >  2:216708/1‑141 (MQ=255)
                          tgAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGc                                                                                                           >  2:216698/1‑141 (MQ=255)
                                     ggTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGt                                                                                                >  1:91341/1‑141 (MQ=255)
                                          gAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAAt                                                                                           <  1:140914/141‑1 (MQ=255)
                                               gCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTtg                                                                                      >  2:443571/1‑141 (MQ=255)
                                               gCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTtg                                                                                      >  1:311026/1‑141 (MQ=255)
                                                 tttATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTtgct                                                                                    >  1:135015/1‑141 (MQ=255)
                                                  ttATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTtgctg                                                                                   >  1:256641/1‑141 (MQ=255)
                                                   tATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTtgctgc                                                                                  <  2:20371/141‑1 (MQ=255)
                                                   tATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAACGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTtgctgc                                                                                  >  2:391263/1‑141 (MQ=255)
                                                        atgatgTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAAtt                                                                             <  2:373906/141‑1 (MQ=255)
                                                          gatgTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTg                                                                           >  2:362981/1‑141 (MQ=255)
                                                             gtgtTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGaaa                                                                        >  2:205592/1‑141 (MQ=255)
                                                               gtTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAATACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAAt                                                                      >  2:261500/1‑141 (MQ=255)
                                                               gtTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAAt                                                                      <  1:291531/141‑1 (MQ=255)
                                                               gtTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAAt                                                                      >  2:174982/1‑141 (MQ=255)
                                                                   gAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACt                                                                  <  1:306111/141‑1 (MQ=255)
                                                                     aGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACttt                                                                >  1:357281/1‑141 (MQ=255)
                                                                      gCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTa                                                               >  1:144138/1‑141 (MQ=255)
                                                                         gTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAAGGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATg                                                            >  1:456338/1‑141 (MQ=255)
                                                                         gTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATg                                                            >  1:30996/1‑141 (MQ=255)
                                                                            aaTGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGaa                                                         <  2:151661/141‑1 (MQ=255)
                                                                               gTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGa                                                      <  2:106281/141‑1 (MQ=255)
                                                                                  gATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGttat                                                   >  2:243079/1‑141 (MQ=255)
                                                                                  gATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGatat                                                   >  2:10024/1‑141 (MQ=255)
                                                                                  gATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGatat                                                   >  1:296324/1‑141 (MQ=255)
                                                                                        aaTGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCtgat                                             >  1:39600/1‑141 (MQ=255)
                                                                                        aaTGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCtgat                                             <  2:217489/141‑1 (MQ=255)
                                                                                               gCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTAcc                                      >  1:156808/1‑141 (MQ=255)
                                                                                                cTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTAccc                                     >  2:5589/1‑141 (MQ=255)
                                                                                                 ttCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTGCGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCt                                    <  1:28789/141‑1 (MQ=255)
                                                                                                    acacTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGtt                                 <  2:135015/141‑1 (MQ=255)
                                                                                                     gacTGTTAGAGGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAAGGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCGGATGATTACCCTGTTc                                <  2:1350/140‑1 (MQ=255)
                                                                                                       cTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCaa                              <  2:453316/141‑1 (MQ=255)
                                                                                                            aGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGt                         >  1:354183/1‑141 (MQ=255)
                                                                                                            aGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGt                         <  1:427928/141‑1 (MQ=255)
                                                                                                              aTGCGACATCATTACAAGAAGTAATACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCa                       <  2:446433/141‑1 (MQ=255)
                                                                                                                 cGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTg                    >  2:72812/1‑141 (MQ=255)
                                                                                                                  gACATCATTACAAGAAGTAATACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTgg                   <  2:278854/141‑1 (MQ=255)
                                                                                                                     atcatTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGc                >  2:68483/1‑141 (MQ=255)
                                                                                                                         ttACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGa            >  1:103127/1‑141 (MQ=255)
                                                                                                                           aCAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAg          <  1:72812/141‑1 (MQ=255)
                                                                                                                           aCAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAg          <  1:375923/141‑1 (MQ=255)
                                                                                                                           aCAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAg          <  1:19235/141‑1 (MQ=255)
                                                                                                                                 agTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAGCGATgg    >  1:266016/1‑141 (MQ=255)
                                                                                                                                   tAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAGCGATGGTg  <  2:404006/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                         
AGCAAAGGACAATACTGATATAGATGTGAAAGTTTTAGGTGGGAAAAGCTTTATTGATGATGTGTTTGAAGCAGTTAATGTAGATCCAAATGATGGCTTCACACTGTTAGATGCGACATCATTACAAGAAGTAACACTTAATGTTAGAACGCATACATTGATTACGCAAGTTTATAGTGCAATGGTTGCTGCTAATTTGAAAATCACTTTAATGGAACGATATCCTGATGATTACCCTGTTCAAATTGTCACTGGTGCACGAAGCGATGGTG  >  CP000730/543026‑543297

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