Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,953,657 Δ1 bp 8.9% coding (841/1383 nt) PP_2586 ← putative amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,953,6570A. 8.9% 135.9 / 7.5 45coding (841/1383 nt)PP_2586putative amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (17/24);  new base . (2/2);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGTCGCGGCCCATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTC  >  NC_002947/2953513‑2953805
                                                                                                                                                |                                                                                                                                                     
cGTCGCGGCCCATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATcc                                                                                                                                                   <  8:302863/149‑1 (MQ=255)
cGTCGCGGCCCATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATcc                                                                                                                                                   <  5:116192/149‑1 (MQ=255)
   cgcgGCCCCTGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGa                                                                                                                                                <  4:57799/149‑1 (MQ=255)
        cccATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGc                                                                                                                                          <  7:220567/149‑1 (MQ=255)
         ccATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCg                                                                                                                                         <  8:8591/149‑1 (MQ=255)
          cATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCg                                                                                                                                         >  7:275297/1‑148 (MQ=255)
              gCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGATCCGGATTCGGCgtcgg                                                                                                                                     <  6:126982/148‑5 (MQ=255)
                             gACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGaa                                                                                                                     <  3:66395/149‑1 (MQ=255)
                             gACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGaa                                                                                                                     <  8:189657/149‑1 (MQ=255)
                                acGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAacac                                                                                                                  <  8:275297/149‑1 (MQ=255)
                                  gCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAAACGGGGAccacca                                                                                                                >  6:166826/1‑149 (MQ=255)
                                   cTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCAt                                                                                                               <  5:134254/149‑1 (MQ=255)
                                      gCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATc                                                                                                              <  4:114461/147‑1 (MQ=255)
                                           ggACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTgc                                                                                                         <  5:244115/147‑1 (MQ=255)
                                               acCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCAGCCCCAATCGGGGAACACCATCTGTGCGCAAt                                                                                                    >  8:163148/1‑148 (MQ=255)
                                                 cATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAAt                                                                                                    >  5:225861/1‑146 (MQ=255)
                                                  aTGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCg                                                                                                                                                  >  5:195598/1‑100 (MQ=255)
                                                  aTGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCg                                                                                                                                                  <  6:195598/100‑1 (MQ=255)
                                                           gCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGAcaccac                                                                                          >  6:197039/1‑146 (MQ=255)
                                                             gAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGa                                                                                      >  6:155564/1‑148 (MQ=255)
                                                                  agcCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCACCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGccc                                                                                                                                >  4:51190/1‑101 (MQ=255)
                                                                  agcCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCACCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGccc                                                                                                                                <  3:51190/101‑1 (MQ=255)
                                                                     ccGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAgg                                                                             >  4:308772/1‑149 (MQ=255)
                                                                                 tagCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAAc                                                                 <  6:225861/149‑1 (MQ=255)
                                                                                         aaGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGc                                                         >  5:267072/1‑149 (MQ=255)
                                                                                         aaGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGc                                                         <  6:267072/149‑1 (MQ=255)
                                                                                            gCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGc                                                         >  6:226313/1‑146 (MQ=255)
                                                                                                ttACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAAc                                                                 >  8:184393/1‑134 (MQ=255)
                                                                                                ttACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAAc                                                                 <  7:184393/134‑1 (MQ=255)
                                                                                                      gcagcTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACg                                            >  5:226533/1‑149 (MQ=255)
                                                                                                              ccgccCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCg                                    <  5:197039/149‑1 (MQ=255)
                                                                                                                    aCCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGt                              >  1:270328/1‑149 (MQ=255)
                                                                                                                    aCCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGt                              >  4:259938/1‑149 (MQ=255)
                                                                                                                       cGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCgg                           <  6:226533/149‑1 (MQ=255)
                                                                                                                          cgcATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGtt                        <  7:214039/149‑1 (MQ=255)
                                                                                                                              tCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTg                    <  5:226313/149‑1 (MQ=255)
                                                                                                                              tCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTg                    <  3:308772/149‑1 (MQ=255)
                                                                                                                                 cGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATGGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGtt                 >  5:255438/1‑149 (MQ=255)
                                                                                                                                       ccGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCa           <  5:155564/149‑1 (MQ=255)
                                                                                                                                         ttcGCCGATCCGGATTCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTAc                      >  7:249111/3‑136 (MQ=255)
                                                                                                                                         ttcGCCGATCCGGATTCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTAc                      <  8:249111/134‑1 (MQ=255)
                                                                                                                                         ttcGCCGATCCGGATTCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGAACAACGACACCGCCAGCATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCATCAGc         >  1:363947/3‑149 (MQ=255)
                                                                                                                                          aCGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCgg                                   <  8:132127/122‑1 (MQ=255)
                                                                                                                                          aCGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCgg                                   >  7:132127/1‑122 (MQ=255)
                                                                                                                                             ccGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATGGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATg     <  6:255438/149‑1 (MQ=255)
                                                                                                                                                aaTCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTc  <  7:163148/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CGTCGCGGCCCATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGATCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTC  >  NC_002947/2953513‑2953805

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: