Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,025 G→T 10.2% A36A (GCG→GCT)  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,0250GT10.2% 131.0 / 6.0 49A36A (GCG→GCT) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base T (2/3);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

ATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCGTTGTC  >  NC_002947/4315883‑4316172
                                                                                                                                              |                                                                                                                                                   
aTCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTc                                                                                                                                               >  2:109438/1‑149 (MQ=255)
   aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcga                                                                                                                                             <  1:867/148‑1 (MQ=255)
     tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatgg                                                                                                                                          <  6:93136/149‑4 (MQ=255)
     tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCATGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGCCGCTGCATTcgatg                                                                                                                                           >  8:102889/1‑146 (MQ=255)
     tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCATGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGCCGCTGCATTcgatg                                                                                                                                           <  7:102889/148‑3 (MQ=255)
          tACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                                      >  6:262424/1‑148 (MQ=255)
            cGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAg                                                                                                                                   >  8:23539/1‑149 (MQ=255)
                  acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTg                                                                                                                             >  3:85148/1‑149 (MQ=255)
                    aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                           <  8:20426/149‑1 (MQ=255)
                                       aCCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATc                                                                                                        >  6:295945/1‑149 (MQ=255)
                                          aGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGaa                                                                                                     <  5:295945/149‑1 (MQ=255)
                                             cATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTg                                                                                                  <  4:85148/149‑1 (MQ=255)
                                                     aCGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTAACCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACt                                                                                          >  1:116802/1‑149 (MQ=255)
                                                       gATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTg                                                                                         >  5:56081/1‑148 (MQ=255)
                                                             tggAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTc                                                                                  <  1:109438/148‑1 (MQ=255)
                                                             gggAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                              <  6:32715/105‑1 (MQ=255)
                                                             gggAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                              >  5:32715/1‑105 (MQ=255)
                                                                    cctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCCCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCa                                                                           >  1:231286/1‑149 (MQ=255)
                                                                             cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCgg                                                                  >  3:62770/1‑149 (MQ=255)
                                                                               cGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCgg                                                                                >  5:91328/1‑133 (MQ=255)
                                                                               cGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCgg                                                                                <  6:91328/133‑1 (MQ=255)
                                                                                   cACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTgc                                                             <  4:62770/148‑1 (MQ=255)
                                                                                      ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTc                                                                                  <  8:161673/124‑1 (MQ=255)
                                                                                      ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTc                                                                                  >  7:161673/1‑124 (MQ=255)
                                                                                          cGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTgc                                                             >  1:254276/1‑141 (MQ=255)
                                                                                          cGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTgc                                                             <  2:254276/141‑1 (MQ=255)
                                                                                           gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCgggg                                                    <  5:262424/149‑1 (MQ=255)
                                                                                           gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCggg                                                     <  3:276783/148‑1 (MQ=255)
                                                                                           gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCgg                                                      <  3:263527/147‑1 (MQ=255)
                                                                                                  gAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTgg                                              <  2:116802/148‑1 (MQ=255)
                                                                                                  gAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGt                                             >  1:206321/1‑149 (MQ=255)
                                                                                                   aaTACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGtt                                            >  4:149314/1‑149 (MQ=255)
                                                                                                       cATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCg                                         <  7:23539/148‑1 (MQ=255)
                                                                                                                tATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGa                               >  2:296241/1‑149 (MQ=255)
                                                                                                                tATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGa                               >  3:47414/1‑149 (MQ=255)
                                                                                                                tATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGa                               <  3:149314/149‑1 (MQ=255)
                                                                                                                 atcatcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAAt                                   >  2:180644/1‑144 (MQ=255)
                                                                                                                 atcatcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAAt                                   <  1:180644/144‑1 (MQ=255)
                                                                                                                          cGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACt                     >  5:287652/1‑149 (MQ=255)
                                                                                                                           gCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTg                    <  1:296241/149‑1 (MQ=255)
                                                                                                                             tACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgcg                  >  3:68906/1‑149 (MQ=255)
                                                                                                                               cccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTACACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTgg                                              <  1:226031/119‑1 (MQ=255)
                                                                                                                               cccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTACACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTgg                                              >  2:226031/1‑119 (MQ=255)
                                                                                                                                tcgtCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTa               >  8:24190/4‑149 (MQ=255)
                                                                                                                                ccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTa               >  7:108503/1‑149 (MQ=255)
                                                                                                                                 cATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCg                                                                   <  7:76828/96‑1 (MQ=255)
                                                                                                                                 cATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCg                                                                   >  8:76828/1‑96 (MQ=255)
                                                                                                                                  aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc             <  6:287652/149‑1 (MQ=255)
                                                                                                                                             cGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCGTTGTc  >  5:193417/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCGTTGTC  >  NC_002947/4315883‑4316172

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