Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,176,020 A→T 9.5% S272T (TCA→ACA)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,0200AT 9.5% 137.7 / 6.1 42S272T (TCA→ACA) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base T (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACAC  >  NC_000913/3175872‑3176163
                                                                                                                                                    |                                                                                                                                               
ggAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGa                                                                                                                                                 >  4:211462/1‑149 (MQ=255)
ggAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGa                                                                                                                                                 >  1:4018/1‑149 (MQ=255)
 gAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGaa                                                                                                                                                <  1:165524/149‑1 (MQ=255)
 gAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGaa                                                                                                                                                >  2:467781/1‑149 (MQ=255)
         cggcggCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATc                                                                                                                                        <  4:11058/149‑1 (MQ=255)
          ggcggcAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCa                                                                                                                                       <  2:126612/149‑1 (MQ=255)
           gcggcAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAg                                                                                                                                      <  3:211462/149‑1 (MQ=255)
            cggcAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGTAGCGGTATCAgg                                                                                                                                     >  4:31407/1‑149 (MQ=255)
            cggcAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGTAGCGGTATCAgg                                                                                                                                     >  4:147942/1‑149 (MQ=255)
                aaCTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTgg                                                                                                                                 >  3:177060/1‑149 (MQ=255)
                   tGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGaaa                                                                                                                              >  3:84052/1‑149 (MQ=255)
                    gCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAAc                                                                                                                             <  4:177060/149‑1 (MQ=255)
                      ggAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACgt                                                                                                                           <  3:385820/149‑1 (MQ=255)
                      ggAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACgt                                                                                                                           <  3:259007/149‑1 (MQ=255)
                      ggAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACgt                                                                                                                           <  1:467781/149‑1 (MQ=255)
                                                   ttCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCggtggt                                                                                              <  4:84052/149‑1 (MQ=255)
                                                       gccagccaGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTc                                                                                                                      >  1:207766/1‑121 (MQ=255)
                                                       gccagccaGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTc                                                                                                                      <  2:207766/121‑1 (MQ=255)
                                                            ccagTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGcc                                                                                     <  3:87865/149‑1 (MQ=255)
                                                                             tCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGtt                                                                                                                                   >  2:433670/1‑86 (MQ=255)
                                                                             tCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGtt                                                                                                                                   <  1:433670/86‑1 (MQ=255)
                                                                                 ttttGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGt                                                                >  1:99858/1‑149 (MQ=255)
                                                                                          gcgcGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAACGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGcc                                                       >  2:343189/1‑149 (MQ=255)
                                                                                                      gTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATgtt                                           >  2:411973/1‑147 (MQ=255)
                                                                                                           aaCCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCa                                      >  3:328858/1‑149 (MQ=255)
                                                                                                                tgtAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTa                                 >  2:167719/1‑149 (MQ=255)
                                                                                                                    aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGAt                                              >  4:351440/1‑132 (MQ=255)
                                                                                                                    aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGAt                                              <  3:351440/132‑1 (MQ=255)
                                                                                                                    aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGt                             >  3:224565/1‑149 (MQ=255)
                                                                                                                    aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGt                             >  3:401921/1‑149 (MQ=255)
                                                                                                                       aaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTgg                          <  2:4018/149‑1 (MQ=255)
                                                                                                                                  gACATCAGTAGCCTTCTGTAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGc               <  3:147942/149‑1 (MQ=255)
                                                                                                                                  gACATCAGTAGCCTTCTGTAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGc               <  3:31407/149‑1 (MQ=255)
                                                                                                                                  gACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGc               <  4:463297/149‑1 (MQ=255)
                                                                                                                                   aCATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCt              >  4:355735/1‑149 (MQ=255)
                                                                                                                                   aCATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCt              >  4:38397/1‑149 (MQ=255)
                                                                                                                                     aTCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa            >  4:290396/1‑149 (MQ=255)
                                                                                                                                       cAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTaaa          <  3:121123/149‑1 (MQ=255)
                                                                                                                                         gTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTg                           <  1:289149/130‑1 (MQ=255)
                                                                                                                                         gTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTg                           >  2:289149/1‑130 (MQ=255)
                                                                                                                                               ttCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGcacac  >  1:232236/1‑149 (MQ=255)
                                                                                                                                               ttCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGcacac  >  2:270537/1‑149 (MQ=255)
                                                                                                                                               ttCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGcacac  >  4:313810/1‑149 (MQ=255)
                                                                                                                                               ttCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGcacac  >  1:187975/1‑149 (MQ=255)
                                                                                                                                                   gAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCa                                                          >  2:121479/1‑89 (MQ=255)
                                                                                                                                                   gAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCa                                                          <  1:121479/89‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACAC  >  NC_000913/3175872‑3176163

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