Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,057 A→G 8.3% P32P (CCA→CCG)  PP_3788 → non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,316,0570AG 8.3% 135.9 / 2.8 48P32P (CCA→CCG) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base G (2/2);  total (26/22)
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.10e-01

CGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCC  >  NC_002947/4315910‑4316206
                                                                                                                                                   |                                                                                                                                                     
cGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATc                                                                                                                                                     >  1:1295313/1‑150 (MQ=255)
       tcatcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCa                                                                                                                                              <  2:1295313/150‑1 (MQ=255)
        catcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGccg                                                                                                                                             >  2:477066/1‑148 (MQ=255)
          tcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGCCGct                                                                                                                                           >  1:1334893/1‑149 (MQ=255)
          tcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCg                                                                                                                                           >  2:217222/1‑150 (MQ=255)
           caGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCTCTACCCATCGAATGCAGCgg                                                                                                                                          >  2:1301636/1‑150 (MQ=255)
             gggAATCAGTCTCGTACGCCACACAGCACGGAGCGCAGCATGAACCAGGCAGTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCa                                                                                                                                        >  2:414688/1‑150 (MQ=255)
              ggAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCAt                                                                                                                                       >  2:426042/1‑150 (MQ=255)
                     gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGCCGCTGCATTcgatgg                                                                                                                                <  2:1334893/150‑4 (MQ=255)
                     gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatgg                                                                                                                                <  1:296476/150‑4 (MQ=255)
                      tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGCCGCTGCATTcgatggg                                                                                                                               <  1:477066/150‑5 (MQ=255)
                      tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                               >  2:63955/1‑150 (MQ=255)
                                 acacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                    <  2:913725/150‑1 (MQ=255)
                                   acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGa                                                                                                                  <  2:1419980/150‑1 (MQ=255)
                                   acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGACATTCGACGAGGACCAGCGCTTGa                                                                                                                  >  1:1047689/1‑150 (MQ=255)
                                           gAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGc                                                                                                          >  1:627154/1‑150 (MQ=255)
                                              cgcAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGaa                                                                                                       <  1:865014/150‑1 (MQ=255)
                                                   caTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAg                                                                                                  <  1:426042/150‑1 (MQ=255)
                                                      gAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGAt                                                                                               <  1:998709/150‑1 (MQ=255)
                                                       aaCCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATc                                                                                              <  2:627154/150‑1 (MQ=255)
                                                         ccAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGa                                                                                            >  2:1037221/1‑150 (MQ=255)
                                                               aTTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCt                                                                                      >  1:853165/1‑150 (MQ=255)
                                                                  aCTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTcac                                                                                   >  2:41444/1‑150 (MQ=255)
                                                                            aCGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCt                                                                         >  2:737254/1‑150 (MQ=255)
                                                                                 aGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCt                                                                    <  1:1037221/150‑1 (MQ=255)
                                                                                    gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCg                                                                                                                                           <  2:708451/76‑1 (MQ=255)
                                                                                    gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCg                                                                                                                                           >  1:708451/1‑76 (MQ=255)
                                                                                        gcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGaa                                                             <  1:414688/150‑1 (MQ=255)
                                                                                          ctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGTAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACt                                                           >  1:104510/1‑150 (MQ=255)
                                                                                           tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGaa                                                                                           <  2:1093978/117‑1 (MQ=255)
                                                                                           tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGaa                                                                                           >  1:1093978/1‑117 (MQ=255)
                                                                                           tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCt                                                                    <  1:665701/140‑1 (MQ=255)
                                                                                           tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCt                                                                    >  2:665701/1‑140 (MQ=255)
                                                                                              cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCa                                                                 <  2:1313524/140‑1 (MQ=255)
                                                                                              cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCa                                                                 >  1:1313524/1‑140 (MQ=255)
                                                                                                cGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGCGCACTGATCGCCCTGCTGGCCACCAACGTCCGAACTGCTCACACTGCACGCTCAGCTGCATGAACTCGGCGt                                                     >  1:1236384/1‑150 (MQ=255)
                                                                                                cGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGt                                                     <  1:63955/150‑1 (MQ=255)
                                                                                                   gCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTgcg                                                  >  1:491242/1‑150 (MQ=255)
                                                                                                         ctcGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCgg                                            <  1:1301636/150‑1 (MQ=255)
                                                                                                                gTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAAGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTAGGCGTGCGCAGCGGGGAtgtg                                     >  2:85201/1‑150 (MQ=255)
                                                                                                                        cATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                             <  1:41444/150‑1 (MQ=255)
                                                                                                                         aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCa                                                                 <  1:515019/113‑1 (MQ=255)
                                                                                                                         aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCa                                                                 >  2:515019/1‑113 (MQ=255)
                                                                                                                              gCGTATCATCCACCGCTACCCATCGAATGCAGCGACATTCGACGAGGACCAGCGCTTGACCTACAGCGAAATAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGc                       <  2:1047689/150‑1 (MQ=255)
                                                                                                                                   tcatcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAg                  >  2:426970/1‑150 (MQ=255)
                                                                                                                                      tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGc               <  1:737254/150‑1 (MQ=255)
                                                                                                                                        cACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCaa             >  2:203740/1‑150 (MQ=255)
                                                                                                                                          ccGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACt           >  1:820790/1‑150 (MQ=255)
                                                                                                                                                   aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATcc  >  1:520923/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCC  >  NC_002947/4315910‑4316206

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