Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,606,472 G→A 7.8% D328N (GAT→AAT)  PP_2284 → tail tubular protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,606,4720GA 7.8% 142.9 / 2.1 51D328N (GAT→AAT) PP_2284tail tubular protein B
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (2/2);  total (24/27)
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.63e-01

CTGACAGCACTGGCGTGTGGCGTGAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTACCGCACCACTGTGA  >  NC_002947/2606326‑2606619
                                                                                                                                                  |                                                                                                                                                   
cTGACAGCACTGGCGTGTGGCGTGAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGAt                                                                                                                                                   <  5:271996/149‑1 (MQ=255)
  gACAGCACTGGCGTGTGGCGTGAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGc                                                                                                                                                 >  1:23225/1‑149 (MQ=255)
   aCAGCACTGGCGTGTGGCGTGAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCg                                                                                                                                                <  4:273423/149‑1 (MQ=255)
        aCTGGCGTGTGGCGTGAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGacac                                                                                                                                           <  3:18952/149‑1 (MQ=255)
            gcgtgtggcgtgAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTGAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCaa                                                                                                                                        <  2:358273/148‑1 (MQ=255)
                   gcgtgAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAAGGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGAc                                                                                                                                 <  7:248401/148‑1 (MQ=255)
                     gtgAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCAATGCGGACACCAACGAAGAccc                                                                                                                               >  4:64884/1‑148 (MQ=255)
                     gtgAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCAATGCGGACACCAACGAAGAccc                                                                                                                               >  4:59995/1‑148 (MQ=255)
                     gtgAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCAATGCGGACACCAACGAAGAccc                                                                                                                               <  3:64884/148‑1 (MQ=255)
                     gtgAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCAATGCGGACACCAACGAAGAccc                                                                                                                               <  3:59995/148‑1 (MQ=255)
                            tgtTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCttt                                                                                                                        <  2:149723/148‑1 (MQ=255)
                            tgtTGGCCCTGGTGTGGCACTGGGCCTGGACGGCGCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCttt                                                                                                                        <  6:339875/148‑1 (MQ=255)
                              ttGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTtgtg                                                                                                                     <  3:252713/149‑1 (MQ=255)
                                    cTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCg                                                                                                                >  4:237637/1‑148 (MQ=255)
                                      ggTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGa                                                                                                              >  8:207138/1‑148 (MQ=255)
                                            gCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGaa                                                                                                                                    >  8:255341/1‑120 (MQ=255)
                                            gCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGaa                                                                                                                                    <  7:255341/120‑1 (MQ=255)
                                                  ggCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGAtg                                                                                                  <  1:140839/148‑1 (MQ=255)
                                                    ccTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGAc                                                                                                                                 <  5:349982/115‑1 (MQ=255)
                                                    ccTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGAc                                                                                                                                 >  6:349982/1‑115 (MQ=255)
                                                             cTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTAcc                                                                                      <  2:23225/149‑1 (MQ=255)
                                                                  ccATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAAc                                                                                 >  6:137747/1‑149 (MQ=255)
                                                                  ccATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAAc                                                                                 >  2:16301/1‑149 (MQ=255)
                                                                   cATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAAcc                                                                                >  4:46780/1‑149 (MQ=255)
                                                                       ccGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCt                                                                            <  6:222701/149‑1 (MQ=255)
                                                                         gCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTgg                                                                          <  7:207138/149‑1 (MQ=255)
                                                                                 tGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGCACACACTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGAGGCCCGGACAATCAACGATGTTGTCTTCTACAGGAACCGCCTGGGATACCtg                                                                  >  1:270715/1‑148 (MQ=255)
                                                                                   gTCCGTAACGCGGACGGCACCTTCGCCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGc                                                                >  6:71771/1‑149 (MQ=255)
                                                                                   gTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGAGAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGc                                                                >  6:1001/1‑149 (MQ=255)
                                                                                   gTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGc                                                                >  7:253882/1‑149 (MQ=255)
                                                                                      cGTAACGCGGACGGCACCTTCGCCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGa                                                             <  5:71771/149‑1 (MQ=255)
                                                                                       gTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGAt                                                            >  2:190229/1‑149 (MQ=255)
                                                                                       gTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGAt                                                            <  3:237637/149‑1 (MQ=255)
                                                                                         aaCGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGAt                                                            >  8:33386/1‑147 (MQ=255)
                                                                                                  ggCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGAt                                                 >  7:197150/1‑149 (MQ=255)
                                                                                                       cttcaccttcCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATTAGGCCGTGATTTtct                                            >  6:167145/1‑149 (MQ=255)
                                                                                                       cgatttcttcCAGGCTGCCACGTGGACGGACTGAGTGGGAGGCGATGCTGACACCAACGAAGACCCGTGCTTTGTGGGCCGGACTATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTtct                                            <  3:175990/143‑1 (MQ=255)
                                                                                                           accttcCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCATGATTTTCTCGGa                                        <  5:137747/149‑1 (MQ=255)
                                                                                                              ttcCAGGCAGCCACCTGCACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGACCTATGTGGGCCGGACAATCAACGATCTTGTCTACTACCGGAACCGCCAGGGATTCCTTGCCGATGAGGCCGTGATGTtct                                            >  4:175990/1‑142 (MQ=255)
                                                                                                                 cAGGCAGCCTCCTGGACGGACAGAGTGGCAGGCGATGCGGACACCTACGAAGACCGGTCCTTTGTGGGCCGGAAATTCATCGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCtgg                                  <  2:270715/149‑1 (MQ=255)
                                                                                                                        ccACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGa                                                                        >  4:287794/1‑104 (MQ=255)
                                                                                                                        ccACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGa                                                                        <  3:287794/104‑1 (MQ=255)
                                                                                                                        ccACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATAc                           <  8:253882/149‑1 (MQ=255)
                                                                                                                               gacggacAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACt                    >  7:168096/1‑149 (MQ=255)
                                                                                                                                  ggacAGAGTGGCAGGCGATGCTGACACCAACGAAGACCCGTCCTTTGTGGGCCGTACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACctct                 <  7:74893/149‑1 (MQ=255)
                                                                                                                                    acAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGGCGGACAATCAACGATTTTGTCTTCTACCGGAACCGCCTGTGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTAc               <  5:1001/149‑1 (MQ=255)
                                                                                                                                    acAGAGTGGCAGGCGATGCGGACACAAACGAAGACCAGTACTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTACTTGCCGATGAGGACGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTAc               >  5:182319/1‑149 (MQ=255)
                                                                                                                                      agagTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGCGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTACCg             <  2:336836/149‑1 (MQ=255)
                                                                                                                                        agTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTACCGca           >  5:99710/1‑149 (MQ=255)
                                                                                                                                         gTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTACCGcac          >  2:36840/1‑149 (MQ=255)
                                                                                                                                                 cGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGTAAATACTGGAACTTCTACCGCACCACTGTGa  <  1:36840/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGACAGCACTGGCGTGTGGCGTGAGTGTGTTGGCCCTGGTGTGGCACTGGGCCTGGACGGCTCAACCATGCCGCACACGCTGGTCCGTAACGCGGACGGCACCTTCACCTTCCAGGCAGCCACCTGGACGGACAGAGTGGCAGGCGATGCGGACACCAACGAAGACCCGTCCTTTGTGGGCCGGACAATCAACGATGTTGTCTTCTACCGGAACCGCCTGGGATTCCTTGCCGATGAGGCCGTGATTTTCTCGGAGTCGGGGAAATACTGGAACTTCTACCGCACCACTGTGA  >  NC_002947/2606326‑2606619

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