Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,268,928 T→A 8.2% E139D (GAA→GAT) ‡ PP_3740 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,268,9280TA 8.2% 140.4 / 3.1 49E139D (GAA→GAT) ‡PP_3740MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base A (2/2);  total (23/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.41e-01

TCACCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCA  >  NC_002947/4268783‑4269070
                                                                                                                                                 |                                                                                                                                              
tCACCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGc                                                                                                                                            >  1:551896/1‑150 (MQ=255)
  aCCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTg                                                                                                                                          <  2:639483/150‑1 (MQ=255)
   ccAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTgg                                                                                                                                         >  2:186676/1‑150 (MQ=255)
      gTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGAACGCTGGcgc                                                                                                                                      <  2:613688/150‑3 (MQ=255)
         cGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCTCCAATGCCCCAGCAGCCAGACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGAACGCTGGcgcggc                                                                                                                                   <  2:285276/150‑6 (MQ=255)
         cGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATAtt                                                                                                                                   <  2:192694/150‑1 (MQ=255)
              acgacgCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCATCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTg                                                                                                                              <  1:42029/150‑1 (MQ=255)
               cgacgCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTgg                                                                                                                             <  2:761924/150‑1 (MQ=255)
                   gccgccAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAAt                                                                                                                         >  2:698394/1‑150 (MQ=255)
                   gccgccAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAAt                                                                                                                         <  1:556590/150‑1 (MQ=255)
                           aCAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATtgc                                                                                                                 >  1:570969/1‑150 (MQ=255)
                                  tACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGc                                                                                                          >  1:622267/1‑150 (MQ=255)
                                        tCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGAt                                                                                                    <  2:633792/150‑1 (MQ=255)
                                            cAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCGCCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCa                                                                                                <  2:570969/150‑1 (MQ=255)
                                            cAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCa                                                                                                >  2:486448/1‑150 (MQ=255)
                                            cAGCCAAAATGCCCCAGAACGTGTACCGCCAGCAAACCGCCAAGCGTCGCACCCAGCGCATACCCAGTGGACTGCAGGCTCACCTCCAAGCCGCGCCAGCGTTCGCTTGCATATTCACTAGCAATCACATTGCTGCTGGCCAGGATGCCa                                                                                                <  1:44346/150‑1 (MQ=255)
                                               ccccAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCg                                                                                             <  2:596391/150‑1 (MQ=255)
                                               ccccAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCg                                                                                             <  2:596355/150‑1 (MQ=255)
                                                cccAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACAtt                                                                                                                   <  2:158880/127‑1 (MQ=255)
                                                cccAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACAtt                                                                                                                   >  1:158880/1‑127 (MQ=255)
                                                  cAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATg                                                                                          <  2:622267/150‑1 (MQ=255)
                                                         ccAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGTATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGcc                                                                                   >  1:159152/1‑150 (MQ=255)
                                                                       gCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCTGTCCATGACGCGTTCGATCTCCGACCCTGATCAGCTGGCCAGGATGCCCCCGATGCCCAGCCCGTTCAGGCcgctc                                                                     >  1:101084/1‑148 (MQ=255)
                                                                        ccAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCgca                                                                    >  1:230890/1‑150 (MQ=255)
                                                                             aaaCCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCagcagc                                                               >  2:780217/1‑150 (MQ=255)
                                                                               aCCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGcc                                                                         <  1:784617/138‑1 (MQ=255)
                                                                               aCCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGcc                                                                         >  2:784617/1‑138 (MQ=255)
                                                                                   ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCGAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGcc                                                                                   >  1:348067/1‑124 (MQ=255)
                                                                                   ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCGAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGcc                                                                                   <  2:348067/124‑1 (MQ=255)
                                                                                   ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGa                                                                                                     <  2:9179/106‑1 (MQ=255)
                                                                                   ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGa                                                                                                     >  1:9179/1‑106 (MQ=255)
                                                                                   ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGt                                                         >  2:743530/1‑150 (MQ=255)
                                                                                   ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCATCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGt                                                         >  2:76417/1‑150 (MQ=255)
                                                                                    cAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGtt                                                        <  1:743530/150‑1 (MQ=255)
                                                                                          tCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCt                                                                                                                                           <  2:510815/61‑1 (MQ=255)
                                                                                          tCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCt                                                                                                                                           >  1:510815/1‑61 (MQ=255)
                                                                                                 cAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGcc                                                                                        <  2:844390/105‑1 (MQ=255)
                                                                                                 cAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGcc                                                                                        >  1:844390/1‑105 (MQ=255)
                                                                                                  aGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTgg                                          <  1:672444/150‑1 (MQ=255)
                                                                                                          gCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGc                                  <  2:101084/150‑1 (MQ=255)
                                                                                                                ggACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCagagag                            <  1:486448/150‑1 (MQ=255)
                                                                                                                  aCTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGagca                          >  1:153364/1‑150 (MQ=255)
                                                                                                                      cAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCa                                 >  1:522277/1‑139 (MQ=255)
                                                                                                                      cAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCa                                 <  2:522277/139‑1 (MQ=255)
                                                                                                                        ggCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGc                    <  1:698394/150‑1 (MQ=255)
                                                                                                                           tCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGa                 >  2:204317/1‑150 (MQ=255)
                                                                                                                                     aatatgcCAGCGAACGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAg       >  2:727213/6‑150 (MQ=255)
                                                                                                                                       tatgcCAGCGAACGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAgcg     >  2:149599/4‑150 (MQ=255)
                                                                                                                                          gcCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCa  <  2:153364/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCACCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCA  >  NC_002947/4268783‑4269070

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