Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,857,970 G→T 11.0% T4N (ACT→AAT) ‡ PP_4272 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,857,9700GT11.0% 123.8 / 3.8 45T4N (ACT→AAT) ‡PP_4272membrane protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base T (3/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGT  >  NC_002947/4857826‑4858119
                                                                                                                                                |                                                                                                                                                      
cTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCaaa                                                                                                                                                   >  1:917378/1‑149 (MQ=255)
     cagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCATTCTGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaa                                                                                                                                             <  1:660051/150‑1 (MQ=255)
      agcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaag                                                                                                                                            <  1:19358/150‑1 (MQ=255)
              ccACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGc                                                                                                                                    >  2:946408/1‑150 (MQ=255)
              ccACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGc                                                                                                                                    >  1:377580/1‑150 (MQ=255)
               cacaGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCAAAAAAGTTGAAGTCAACTACAAGTCAAGGGcc                                                                                                                                   >  2:308441/1‑150 (MQ=255)
                     tacATCAGTAAAGTGACACCAAAGCAAATCGGCGATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAATCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAAt                                                                                                                             <  2:232218/148‑1 (MQ=255)
                       cATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                           <  1:1292023/150‑1 (MQ=255)
                       cATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                           <  1:789316/150‑1 (MQ=255)
                        aTCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCGTGCCGCGCAGCATAGAAGTTGAAGTCAGCTTCAAGTCAAGGGCCCGAAATGGa                                                                                                                          >  2:1157412/1‑150 (MQ=255)
                          cACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACa                                                                                                                        <  1:308441/150‑1 (MQ=255)
                                       acctggCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGc                                                                                                           >  2:895134/6‑150 (MQ=255)
                                        cAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                           >  2:1125554/1‑133 (MQ=255)
                                        cAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATgg                                                                                                                           <  1:1125554/133‑1 (MQ=255)
                                        aaaTGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCCCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGcc                                                                                                          <  2:568306/149‑1 (MQ=255)
                                               aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCaaa                                                                                                       >  1:1205114/1‑146 (MQ=255)
                                               aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCaaa                                                                                                       <  2:1205114/146‑1 (MQ=255)
                                                  cGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTg                                                                                                <  1:946408/150‑1 (MQ=255)
                                                        tGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTAc                                                                                          >  2:403666/1‑150 (MQ=255)
                                                             cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCg                                                                                                                  >  1:1211173/1‑121 (MQ=255)
                                                             cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCg                                                                                                                  <  2:1211173/121‑1 (MQ=255)
                                                              tAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaagtcaa                                                                                                                                        <  1:564219/98‑1 (MQ=255)
                                                              tAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtcaagtcaa                                                                                                                                        >  2:564219/1‑98 (MQ=255)
                                                                    cATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGacac                                                                              <  2:775838/150‑1 (MQ=255)
                                                                      tCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg                                                                            >  2:774272/1‑150 (MQ=255)
                                                                          aGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGtt                                                                        >  2:1244146/1‑150 (MQ=255)
                                                                                        gAACAAGCCCTCGAAGGTTGACGGTCATTACCGTGCCGCGCAGCATAGAAGTTGAAGTCAGCTTCACAGTAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGaaa                                                          <  1:1157412/150‑1 (MQ=255)
                                                                                        gAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGACTTCAACTtctat                                                                                                                                             >  2:1053355/1‑65 (MQ=38)
                                                                                        gAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGACTTCAACTtctat                                                                                                                                             <  1:1053355/68‑4 (MQ=38)
                                                                                           cAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGc                                                       <  1:1244146/150‑1 (MQ=255)
                                                                                           cAAGCCCGCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGc                                                       <  1:567834/150‑1 (MQ=255)
                                                                                              gCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTg                                                    >  2:361757/1‑150 (MQ=255)
                                                                                                 cTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTc                                                 <  2:377580/150‑1 (MQ=255)
                                                                                                  tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCa                                                >  2:236083/1‑150 (MQ=255)
                                                                                                    gAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTAAAGGGCCCGAAATGGACATTGCCGACGTCTCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGACGAAAGGCTTGCTCAAg                                              >  2:1366650/1‑150 (MQ=255)
                                                                                                                tCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCggg                                  >  1:817944/1‑150 (MQ=255)
                                                                                                                           ccGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaagg                       <  1:774272/150‑4 (MQ=255)
                                                                                                                             gcgcAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTc                     <  2:817944/150‑1 (MQ=255)
                                                                                                                                     ttgAAGTTGCCTTCAACTTCAAGTCAAGGGCCCGAGATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGttt             >  2:59822/3‑150 (MQ=255)
                                                                                                                                      tgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg                                                                            <  2:324015/85‑1 (MQ=255)
                                                                                                                                      tgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg                                                                            >  1:324015/2‑86 (MQ=255)
                                                                                                                                      aGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTg            >  1:45114/1‑150 (MQ=255)
                                                                                                                                          gTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGcc                                     <  2:450076/121‑1 (MQ=255)
                                                                                                                                          gTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGcc                                     >  1:450076/1‑121 (MQ=255)
                                                                                                                                             gAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGa     <  1:361757/150‑1 (MQ=255)
                                                                                                                                                gTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGAtgt  <  2:917378/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGT  >  NC_002947/4857826‑4858119

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