Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,206,999 G→A 7.5% I259I (ATC→ATT)  PP_3691 ← DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,206,9990GA 7.5% 149.7 / 2.1 53I259I (ATC→ATT) PP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base A (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

TTCTGCGTACTGAGCGTTGTTTCAAATAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCG  >  NC_002947/4206853‑4207135
                                                                                                                                                  |                                                                                                                                         
ttCTGCGTACTGAGCGTTGTTTCAAATAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGa                                                                                                                                          <  3:325983/148‑1 (MQ=255)
     cGTACTGAGCGTTGTTTCAAATAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCg                                                                                                                                    <  5:322515/149‑1 (MQ=255)
      gTACTGAGCGTTGTTTCAAATAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCgg                                                                                                                                   <  2:73248/149‑1 (MQ=255)
                       aaaTAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTACAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCt                                                                                                                  >  3:116058/1‑149 (MQ=255)
                        aaTAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTc                                                                                                                 >  7:174876/1‑149 (MQ=255)
                               ccAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAggg                                                                                                            >  5:284185/1‑147 (MQ=255)
                                  gCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAgcg                                                                                                       >  5:130969/1‑149 (MQ=255)
                                  gCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAgcg                                                                                                       <  8:162238/149‑1 (MQ=255)
                                       cGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑aa                                                                                                 <  8:174876/149‑1 (MQ=255)
                                        gCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑aa                                                                                                 >  7:45801/1‑148 (MQ=255)
                                            cTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AAcagca                                                                                            <  5:192592/149‑1 (MQ=255)
                                             atCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCATCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCCGGTTGAGCCCGATTTGCTCAAGGCCAGCGCAG‑AAcagcag                                                                                           >  8:132116/2‑149 (MQ=255)
                                                aCCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGt                                                                                        >  2:277064/1‑149 (MQ=255)
                                                  ccGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGa                                                                                                                              <  6:149650/110‑1 (MQ=255)
                                                  ccGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACAGCGAGTGAATACGATGACCGGGTTGa                                                                                                                              >  5:149650/1‑110 (MQ=255)
                                                                ccgTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTa                                                                        <  1:277064/149‑1 (MQ=255)
                                                                        tcatcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGc                                                                                                                   <  1:213181/99‑1 (MQ=255)
                                                                        tcatcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGc                                                                                                                   >  2:213181/1‑99 (MQ=255)
                                                                           tcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAggg                                                                                                            >  6:267170/1‑103 (MQ=255)
                                                                           tcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAggg                                                                                                            <  5:267170/103‑1 (MQ=255)
                                                                               ttAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTAt                                                           >  3:262809/1‑147 (MQ=255)
                                                                                      aGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCt                                                  <  6:284185/149‑1 (MQ=255)
                                                                                       gTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTg                                                 <  6:130969/149‑1 (MQ=255)
                                                                                          gTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAg                                                                           <  7:90348/120‑1 (MQ=255)
                                                                                          gTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAg                                                                           >  8:90348/1‑120 (MQ=255)
                                                                                          gTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGt                                                                         >  5:161089/1‑122 (MQ=255)
                                                                                          gTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGt                                                                         <  6:161089/122‑1 (MQ=255)
                                                                                          gTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCt                                              >  6:360252/1‑149 (MQ=255)
                                                                                            gAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACAATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCa                                                         >  5:36214/1‑136 (MQ=255)
                                                                                            gAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACAATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCa                                                         <  6:36214/136‑1 (MQ=255)
                                                                                              aCTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTggg                                          <  3:214674/149‑1 (MQ=255)
                                                                                              aCTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAAGACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTggg                                          <  8:154819/149‑1 (MQ=255)
                                                                                                    ctGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGt                                                             >  5:83483/1‑124 (MQ=255)
                                                                                                    ctGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGt                                                             <  6:83483/124‑1 (MQ=255)
                                                                                                     tGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACg                                     <  8:45801/147‑1 (MQ=255)
                                                                                                     tGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACAATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTa                                       <  5:286464/145‑1 (MQ=255)
                                                                                                     tGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACAATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTa                                       >  6:286464/1‑145 (MQ=255)
                                                                                                            aTCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAAACAGCAGGGTGAGATCCATGCAGGTAGCCAAAc                                                                 <  1:229817/113‑1 (MQ=255)
                                                                                                            aTCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAAc                                                                 >  2:229817/1‑112 (MQ=255)
                                                                                                            aTCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtttt                             <  7:132116/148‑1 (MQ=255)
                                                                                                                    tCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTgaga                                                                                    <  2:73254/85‑1 (MQ=255)
                                                                                                                    tCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTgaga                                                                                    >  1:73254/1‑85 (MQ=255)
                                                                                                                      cTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAg                                                                           >  7:35196/1‑92 (MQ=255)
                                                                                                                      cTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAg                                                                           <  8:35196/92‑1 (MQ=255)
                                                                                                                      cTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGt                  >  4:315237/1‑149 (MQ=255)
                                                                                                                      cTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGt                  >  1:105189/1‑149 (MQ=255)
                                                                                                                         cGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGtt                 <  2:105189/147‑1 (MQ=255)
                                                                                                                         cGAAAGAGTGACCGCTAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCGAGGGCAGCGCAG‑AACAGCAGGGTCAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGTCTTGGGCTACGCACCTTTTGACCTGACTGtt                 <  4:116058/147‑1 (MQ=255)
                                                                                                                                    cgtagAGTGAATACGATGAGAGGGTTTAGCCCGATTTTCTCCAGGGCAGCGCAG‑AACAGCAGTGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGTCTACTCACCTTTTGACCTGTCTGTTCAAAGctggctgg    <  4:262809/145‑1 (MQ=255)
                                                                                                                                     cgcgAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCtggctggc   >  4:21440/1‑149 (MQ=255)
                                                                                                                                      gcgAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTggcg  <  5:360252/149‑1 (MQ=255)
                                                                                                                                           tGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtttt                             >  6:218708/1‑117 (MQ=255)
                                                                                                                                           tGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtttt                             <  5:218708/117‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
TTCTGCGTACTGAGCGTTGTTTCAAATAGAACCAGCTCCCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAG‑AACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCG  >  NC_002947/4206853‑4207135

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