Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,722,424 G→A 7.2% Q11* (CAG→TAG)  PP_5598 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,722,4240GA 7.2% 159.6 / 2.0 55Q11* (CAG→TAG) PP_5598hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/29);  new base A (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

TATGGGGCGACCAATCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATGGCACGCCTGCGGCGCCTGAG  >  NC_002947/3722283‑3722564
                                                                                                                                             |                                                                                                                                            
tATGGGGCGACCAATCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTa                                                                                                                                       <  1:124982/149‑1 (MQ=255)
gaTGGGGCGACCAATCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTa                                                                                                                                       <  8:88501/148‑1 (MQ=255)
    gggCGACCAATCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTAttct                                                                                                                                   >  8:201668/1‑149 (MQ=255)
            aaTCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACg                                                                                                                           >  2:147835/1‑149 (MQ=255)
             aTCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGt                                                                                                                          <  4:28813/149‑1 (MQ=255)
                agcaagcaTGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTaaa                                                                                                                       <  6:172964/149‑1 (MQ=255)
                 gcaagcaTGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTACTCAGTATTCTTCTCGACGTAAAc                                                                                                                      >  3:270694/1‑149 (MQ=255)
                   aagcaTGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGt                                                                                                                    >  4:324834/1‑149 (MQ=255)
                   aagcaTGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGt                                                                                                                    <  3:346918/149‑1 (MQ=255)
                       aTGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTtc                                                                                                                  <  8:7764/147‑1 (MQ=255)
                         gTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTtcatca                                                                                                              >  6:238037/1‑149 (MQ=255)
                            ccGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAggg                                                                                                           >  8:252112/1‑149 (MQ=255)
                             cgcgAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGc                                                                                                          >  4:363311/1‑149 (MQ=255)
                                 aGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTcc                                                                                                      <  4:77270/149‑1 (MQ=255)
                                     gggggCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCaaa                                                                                                   <  8:352467/148‑1 (MQ=255)
                                      ggggCATTGAGGCGGCTTGTAAGGCGGTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGACACCCGCCCGCCAGTCCGACCTGCGCACCCAACAGCGTCCCTGCTCAGAATTCTGCTCGACGAAAACGATGATCACGGCCTCCAAAg                                                                                                  >  4:290900/1‑148 (MQ=255)
                                      ggggCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCa                                                                                                                                          <  3:71862/108‑1 (MQ=255)
                                      ggggCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCa                                                                                                                                          >  4:71862/1‑108 (MQ=255)
                                       gggCATTGATGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCGACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCa                                                                                                <  1:46819/149‑1 (MQ=255)
                                              gAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATc                                                                                          <  3:384345/148‑1 (MQ=255)
                                                ggcggcTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCgg                                                                                        >  5:86221/1‑148 (MQ=255)
                                                 gcggcTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGTCGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTa                                                                                      <  4:374381/149‑1 (MQ=255)
                                                  cggcTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAg                                                                                     <  7:201668/149‑1 (MQ=255)
                                                  cggcTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAg                                                                                     <  3:324834/149‑1 (MQ=255)
                                                            gCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAg                                                                            >  8:201384/1‑148 (MQ=255)
                                                              cGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAAt                                                                         >  1:1560/1‑149 (MQ=255)
                                                                  cGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGaa                                                                      <  8:119583/148‑1 (MQ=255)
                                                                  cGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGaa                                                                      <  3:363311/148‑1 (MQ=255)
                                                                   gACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTACTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAg                                                                     <  4:270694/148‑1 (MQ=255)
                                                                    aCAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGc                                                                   <  8:331677/149‑1 (MQ=255)
                                                                      aaaCCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGt                                                                 >  5:205234/1‑149 (MQ=255)
                                                                      aaaCCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGt                                                                 >  4:298043/1‑149 (MQ=255)
                                                                        aCCGCTAGGGTTCGGCAGCCGTACTTAAACACACCCGCCAGTCATCCATGCGCACCCATCAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTcc                                                               <  7:222283/149‑1 (MQ=255)
                                                                          cGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCt                                                             >  1:371186/1‑149 (MQ=255)
                                                                                 cTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCa                                                      <  6:138746/149‑1 (MQ=255)
                                                                                     ggCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAgcg                                                   <  6:205234/148‑1 (MQ=255)
                                                                                      gCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCa                                                 <  5:238037/149‑1 (MQ=255)
                                                                                            gTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGaa                                           <  4:167861/149‑1 (MQ=255)
                                                                                                tGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTg                                                            >  7:319274/1‑128 (MQ=255)
                                                                                                tGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTg                                                            <  8:319274/128‑1 (MQ=255)
                                                                                                  aaaCCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAg                                                                     >  7:370882/1‑117 (MQ=255)
                                                                                                  aaaCCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAg                                                                     <  8:370882/117‑1 (MQ=255)
                                                                                                   aaCCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAgg                                    <  7:19539/149‑1 (MQ=255)
                                                                                                     cccacccGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGc                                  >  8:332468/1‑149 (MQ=255)
                                                                                                        acccGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAg                                                     <  2:218601/127‑1 (MQ=255)
                                                                                                        acccGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAg                                                     >  1:218601/1‑127 (MQ=255)
                                                                                                        acccGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGaa                                           <  4:176681/137‑1 (MQ=255)
                                                                                                        acccGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGaa                                           >  3:176681/1‑137 (MQ=255)
                                                                                                          ccGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTaa                             <  7:314323/149‑1 (MQ=255)
                                                                                                            gCCAGTCATCCATGCGCACCCAACAGCGTCCCTACTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACg                           >  4:271796/1‑149 (MQ=255)
                                                                                                                   aTCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATgg                     <  7:201384/148‑1 (MQ=255)
                                                                                                                           gcACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATGGCACGCCTgc            >  2:279911/1‑149 (MQ=255)
                                                                                                                               ccAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATGGCACGCCTgcggcg        >  7:337459/1‑149 (MQ=255)
                                                                                                                               ccAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATGGCACGCCTgcggcg        >  6:217207/1‑149 (MQ=255)
                                                                                                                                     gCGTCCCTACTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATGGCACGCCTGCGGCGCCTGAg  <  3:271796/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
TATGGGGCGACCAATCAGCAAGCATGTCCCGCGAGTGGGGGGCATTGAGGCGGCTTGTAAGCCGTTCGACAAACCGCTCGGCTTCGGCAGCCGTACTGAAACCCACCCGCCAGTCATCCATGCGCACCCAACAGCGTCCCTGCTCAGTATTCTTCTCGACGTAAACGTTCATCAGGGCCTCCAAAGCATGTATCGGTAGACAATGCAGAATGAAGGCGTCCCTGCCGGCAGCGCAGGTGAAGCTCAGGGCGTTAACGCCATGGCACGCCTGCGGCGCCTGAG  >  NC_002947/3722283‑3722564

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