Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,436,026 C→T V180V (GTG→GTA)  cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,436,0260CT96.0% 84.9 / ‑4.1 25V180V (GTG→GTA) cpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (11/13);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

ATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCG  >  NC_000913/4435881‑4436168
                                                                                                                                                 |                                                                                                                                              
aTAGCCCGGAATGTGCCAGAAAGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTacga                                                                                                                                            <  2:265382/150‑1 (MQ=255)
   gCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGc                                                                                                                                         >  1:255063/1‑150 (MQ=255)
         aaTGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAatat                                                                                                                                   <  2:235294/150‑1 (MQ=255)
            gtgCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGc                                                                                                                                >  1:286528/1‑150 (MQ=255)
            gtgCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGc                                                                                                                                >  2:121221/1‑150 (MQ=255)
                   aaCGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCt                                                                                                                          <  2:244647/149‑1 (MQ=255)
                       acaacaacaTCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGc                                                                                                                     <  2:127908/150‑1 (MQ=255)
                                                   tttCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCa                                                                                         >  1:267916/1‑150 (MQ=255)
                                                   tttCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCa                                                                                         >  1:102194/1‑150 (MQ=255)
                                                                    cgcACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTttaatta                                                                       <  2:60595/151‑1 (MQ=255)
                                                                          gTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACg                                                                 >  2:79163/1‑151 (MQ=255)
                                                                          gTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACg                                                                 >  2:297501/1‑151 (MQ=255)
                                                                          gTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACg                                                                 >  2:229904/1‑151 (MQ=255)
                                                                                  aaTATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAAc                                                         <  1:137038/151‑1 (MQ=255)
                                                                                   atatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACa                                                        <  1:32137/151‑1 (MQ=255)
                                                                                   atatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACa                                                        >  1:132114/1‑151 (MQ=255)
                                                                                    tatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACAt                                                       <  2:140222/151‑1 (MQ=255)
                                                                                                                    gCTTTATCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTgg                       >  2:233827/1‑151 (MQ=255)
                                                                                                                         aTCCCAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCAttt                  >  1:40842/1‑151 (MQ=255)
                                                                                                                           cccAGCCCATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTAc                <  1:249647/151‑1 (MQ=255)
                                                                                                                                cccATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAg           <  2:155747/151‑1 (MQ=255)
                                                                                                                                  cATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAggg         <  1:44838/151‑1 (MQ=255)
                                                                                                                                  cATGATTTGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAggg         <  2:255061/151‑1 (MQ=255)
                                                                                                                                       tttGTGGTGGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAAttt    <  2:308886/151‑1 (MQ=255)
                                                                                                                                          gtggtgGTACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCg  >  2:197252/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCG  >  NC_000913/4435881‑4436168

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