Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,642,356 T→C 8.5% E283E (GAA→GAG)  PP_3209 ← isoxanthopterin deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,642,3560TC 8.5% 128.5 / 2.9 47E283E (GAA→GAG) PP_3209isoxanthopterin deaminase
Reads supporting (aligned to +/- strand):  ref base T (24/19);  new base C (2/2);  total (26/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 = 9.84e-01

CGATCGAGACCGACATGCCTGCGGCCTCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCAAATGGCTGTG  >  NC_002947/3642235‑3642499
                                                                                                                         |                                                                                                                                               
cGATCGAGACCGACATGCCTGCGGCCTCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTgc                                                                                                                      <  5:97967/149‑1 (MQ=255)
     gagaCCGACATGCCTGCGGCCTCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGaa                                                                                                                 >  7:230423/1‑149 (MQ=255)
                       gCCTCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAg                                                                                               >  5:285980/1‑149 (MQ=255)
                       gCCTCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAg                                                                                               >  1:311049/1‑149 (MQ=255)
                        ccTCCATGGCAACGATATCTGAGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGc                                                                                              >  3:188834/1‑149 (MQ=255)
                          tCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCa                                                                                            >  7:65615/1‑149 (MQ=255)
                          tCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCa                                                                                            >  7:100702/1‑149 (MQ=255)
                             aTGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTc                                                                                         >  3:209523/1‑149 (MQ=255)
                              tGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGCTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCATGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTcc                                                                                        >  3:143697/1‑149 (MQ=255)
                                   aCGATATCTGAGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCt                                                                                    <  4:188834/148‑1 (MQ=255)
                                          cTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCcgcg                                                                               >  7:25799/1‑146 (MQ=255)
                                            gCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATCTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGaa                                                                          >  8:147587/1‑149 (MQ=255)
                                                  ccACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCNGAGCCATTCCTGCTCCGCGCAGAAcaccac                                                                    <  1:30232/149‑1 (MQ=255)
                                                  ccACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAAcaccac                                                                    <  1:24821/149‑1 (MQ=255)
                                                   cACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACg                                                                   >  2:19819/1‑149 (MQ=255)
                                                         ctaGCCGGCCATTGCTTTGTTGGCAGTGGGGGATGCCGGTGCCGGTGGCTCCCAGCAGCTGGATTTCTTCGGCTCTCAGCTTGACCAGGTGCGGGAACCACAGGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACTGGCgag                                                             <  4:143697/147‑1 (MQ=255)
                                                          aaGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCgaga                                                            <  2:6074/149‑1 (MQ=255)
                                                            gCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAgg                                                          <  8:100702/149‑1 (MQ=255)
                                                               ggCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGc                                                                                                    >  1:248133/1‑104 (MQ=255)
                                                               ggCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGc                                                                                                    <  2:248133/104‑1 (MQ=255)
                                                                 ccATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTg                                                     <  2:243927/149‑1 (MQ=255)
                                                                  cATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCTCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTgg                                                    <  5:69396/149‑1 (MQ=255)
                                                                  cATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAgg                                                          <  3:179782/143‑1 (MQ=255)
                                                                  cATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAgg                                                          >  4:179782/1‑143 (MQ=255)
                                                                  cATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTgg                                                    <  3:315023/149‑1 (MQ=255)
                                                                  cATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTgg                                                    <  1:19819/149‑1 (MQ=255)
                                                                    ttGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGc                                                   >  4:159069/1‑148 (MQ=255)
                                                                        tttGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCa                                              >  2:319984/1‑149 (MQ=255)
                                                                            tgGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATgg                                           >  4:4935/1‑148 (MQ=255)
                                                                                 aGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGTCGCTCAGCTTGACCAGGTGCGCGAAc                                                                                                                <  6:330490/74‑1 (MQ=255)
                                                                                 aGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAAc                                                                                                                >  5:330490/1‑74 (MQ=255)
                                                                                    gggCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATCTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTc                                  <  7:147587/149‑1 (MQ=255)
                                                                                      gCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTg                                <  1:319984/149‑1 (MQ=255)
                                                                                       cGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTgg                               <  4:209523/149‑1 (MQ=255)
                                                                                       cGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTgg                               >  6:208450/1‑149 (MQ=255)
                                                                                             cggtgccggtgGCGCCCAGCAGCTGGATCTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTc                                  <  4:78656/140‑1 (MQ=255)
                                                                                             cggtgccggtgGCGCCCAGCAGCTGGATCTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTc                                  >  3:78656/1‑140 (MQ=255)
                                                                                               gtgccggtgGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCAc                       >  5:135261/1‑149 (MQ=255)
                                                                                                  ccggtgGCTCCCAGAAGCTGGATTTCTTCGGCGCTCAGCTTGAACAGGTGCGCGAAACACACGTCCGGCCCCACCCATTCCTGCTCCGCGCAGAACACCACGGGAGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGt                    >  2:786/1‑149 (MQ=255)
                                                                                                      tgGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCg                <  3:159069/149‑1 (MQ=255)
                                                                                                          gcCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTc             >  2:194841/1‑148 (MQ=255)
                                                                                                           cccAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCaa           <  8:25799/149‑1 (MQ=255)
                                                                                                            ccAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCaaa          >  5:289435/1‑149 (MQ=255)
                                                                                                             cagcagCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCaaa          <  4:86946/148‑1 (MQ=255)
                                                                                                             cagcagCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCaaa          >  3:86946/1‑148 (MQ=255)
                                                                                                                 agcTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCAAATGGCt     >  3:106241/1‑149 (MQ=255)
                                                                                                                      gATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCAAATGGCtgtg  >  2:330279/1‑147 (MQ=255)
                                                                                                                         |                                                                                                                                               
CGATCGAGACCGACATGCCTGCGGCCTCCATGGCAACGATATCTGCGATACCACTGCCAAGCCGGCCATTGCTTTGTGGGCAGTGGGCGATGCCGGTGCCGGTGGCGCCCAGCAGCTGGATTTCTTCGGCGCTCAGCTTGACCAGGTGCGCGAACCACACGTCCGGCCCGAGCCATTCCTGCTCCGCGCAGAACACCACGGGCGAGAGGCCGTGGCGATCATGGACTTCGTTCTGGTACTCCACGGTTTCGCTCAAATGGCTGTG  >  NC_002947/3642235‑3642499

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: ‑