Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 1,995,953 T→C 7.4% I29M (ATA→ATG)  EO946_RS09955 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349431,995,9530TC 7.4% 162.9 / 3.6 54I29M (ATA→ATG) EO946_RS09955hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/35);  new base C (2/2);  total (17/37)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATCTCCTTTTATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATCATATTAAATCATAAAAAAGTGA  >  NZ_CP034943/1995815‑1996089
                                                                                                                                          |                                                                                                                                        
aTCTCCTTTTATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGtatta                                                                                                                                         >  1:2355466/1‑140 (MQ=255)
 tctcCTTTTATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGtattat                                                                                                                                        >  2:2286340/1‑140 (MQ=255)
        ttATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGt                                                                                                                                 <  2:1569061/140‑1 (MQ=255)
        ttATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGt                                                                                                                                 <  2:1282062/140‑1 (MQ=255)
        ttATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGt                                                                                                                                 <  2:616898/140‑1 (MQ=255)
        ttATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGt                                                                                                                                 <  2:651286/140‑1 (MQ=255)
        ttATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGt                                                                                                                                 <  2:649426/140‑1 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCGTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          <  1:1308132/140‑1 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          <  1:379775/140‑1 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          <  1:377840/140‑1 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          >  2:467189/1‑140 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          >  2:702270/1‑140 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          <  1:1445206/140‑1 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          >  2:872963/1‑140 (MQ=255)
               tAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCa                                                                                                                          <  1:1145413/140‑1 (MQ=255)
                  agCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGtt                                                                                                                       >  1:1169126/1‑140 (MQ=255)
                      tCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCtt                                                                                                                   <  1:872963/140‑1 (MQ=255)
                       cGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCTCAACAACATTATTAGTTCGTATCATACCTCGTTGCTTCAGTTCCTTa                                                                                                                  >  1:1014467/1‑140 (MQ=255)
                       cGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTa                                                                                                                  >  1:1006350/1‑140 (MQ=255)
                       cGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATCATACCTCGTTGCTTCAGTTCCTTa                                                                                                                  >  1:1595867/1‑140 (MQ=255)
                           atTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCTCAACAACATTATTAGTTCGTATCATACCTCGTTGCTTCAGTTCCTTAAGCa                                                                                                              <  2:1014467/140‑1 (MQ=255)
                           atTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCa                                                                                                              <  2:1006350/140‑1 (MQ=255)
                           atTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATCATACCTCGTTGCTTCAGTTCCTTAAGCa                                                                                                              <  2:1595867/140‑1 (MQ=255)
                              ttGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACg                                                                                                           <  2:2355466/140‑1 (MQ=255)
                                         ttGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGc                                                                                                <  1:632629/140‑1 (MQ=255)
                                             tGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATATAGCAAAg                                                                                            >  2:1790020/1‑140 (MQ=255)
                                             tGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAg                                                                                            >  2:577376/1‑140 (MQ=255)
                                                  cTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAAt                                                                                       >  1:2335235/1‑140 (MQ=255)
                                                                  cccTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAAtt                                                                       >  1:580430/1‑140 (MQ=255)
                                                                           tttATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATg                                                              >  1:424256/1‑140 (MQ=255)
                                                                              attatAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGaaa                                                           <  2:1088252/140‑1 (MQ=255)
                                                                              attatAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGaaa                                                           <  2:989486/140‑1 (MQ=255)
                                                                                      aaTTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAg                                                   <  1:1855356/140‑1 (MQ=255)
                                                                                           aTAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAac                                              <  2:1379321/140‑1 (MQ=255)
                                                                                           aTAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAac                                              <  2:760639/140‑1 (MQ=255)
                                                                                           aTAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAac                                              <  1:1496087/140‑1 (MQ=255)
                                                                                           aTAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAac                                              <  2:1985269/140‑1 (MQ=255)
                                                                                           aTAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAac                                              <  2:2289644/140‑1 (MQ=255)
                                                                                             aGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAacca                                            <  1:104968/140‑1 (MQ=255)
                                                                                             aGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAacca                                            <  1:581217/140‑1 (MQ=255)
                                                                                             aGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAacca                                            <  1:1548514/140‑1 (MQ=255)
                                                                                             aGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAacca                                            <  1:751492/140‑1 (MQ=255)
                                                                                                 aaGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATTAAACATTCAAGAAAACCACCCt                                        >  2:2152691/1‑140 (MQ=255)
                                                                                                     aTTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTtta                                    <  2:2402466/140‑1 (MQ=255)
                                                                                                     aTTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTtta                                    <  2:317607/140‑1 (MQ=255)
                                                                                                     aTTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTtta                                    <  1:21076/140‑1 (MQ=255)
                                                                                                     aTTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTtta                                    <  2:1756312/140‑1 (MQ=255)
                                                                                                     aTTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGGTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTtta                                    <  1:467189/140‑1 (MQ=255)
                                                                                                         tcCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTtattat                                <  1:702270/140‑1 (MQ=255)
                                                                                                           cAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTAttattt                              >  1:351068/1‑140 (MQ=255)
                                                                                                           cAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTAttattt                              >  1:1535159/1‑140 (MQ=255)
                                                                                                                 ccccAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATc                        <  1:577376/140‑1 (MQ=255)
                                                                                                                 ccccAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAACTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATc                        <  1:1790020/140‑1 (MQ=255)
                                                                                                                            ttattaGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATCATATTAAATCa             >  1:1680328/1‑140 (MQ=255)
                                                                                                                               ttaGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATCATATTAAATCATaa          <  2:424256/140‑1 (MQ=255)
                                                                                                                                     cGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATCATATTAAATCATAAAAAAGt    >  1:1199093/1‑140 (MQ=255)
                                                                                                                                       tattatACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATCATATTAAATCATAAAAAAGTGa  >  2:400843/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                        
ATCTCCTTTTATACTTAGAGCATCGATATTTTGAGTCCCTGTTGGTGCAGCTTGAAGTTTGGGGAGCCCTTTAGTTTTATTATAAAAATTAATAGCCAAGTATTCTCCAAGCTCCCCAACAACATTATTAGTTCGTATTATACCTCGTTGCTTCAGTTCCTTAAGCAACGAACTATAAAGCAAAGTTAATTCTAAATCATCTAATTTTTTTAATGAAACAGTCAAGAAAACCACCCTTTTATTATTTATAATCATATTAAATCATAAAAAAGTGA  >  NZ_CP034943/1995815‑1996089

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