Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,750,740 A→C 8.6% F87C (TTC→TGC)  PP_3315 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,750,7400AC 8.6% 140.4 / 2.5 46F87C (TTC→TGC) PP_3315membrane protein
Reads supporting (aligned to +/- strand):  ref base A (19/23);  new base C (2/2);  total (21/25)
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.96e-01

GCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACCGCGATCAC  >  NC_002947/3750592‑3750888
                                                                                                                                                    |                                                                                                                                                    
gcCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGa                                                                                                                                                      <  5:29352/149‑1 (MQ=255)
            cgcgACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAAt                                                                                                                                                    <  4:34539/139‑1 (MQ=255)
            cgcgACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAAt                                                                                                                                                    >  3:34539/1‑139 (MQ=255)
                  tttGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATc                                                                                                                                              >  1:346456/1‑139 (MQ=255)
                  tttGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATc                                                                                                                                              <  2:346456/139‑1 (MQ=255)
                    tGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATcagca                                                                                                                                    <  7:8359/147‑1 (MQ=255)
                     ggCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGAGCCACCCGATCAAGAGGACCGGCTGGCCCTGAATGAGGCGGCCATGGCGCACGCCAACGCGCCCGATGGACTGCGTGACTGCGATCGCGGTgagcctcg                                                                                                                                 >  4:293658/1‑145 (MQ=255)
                                 aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATcgt                                                                                                                                            <  7:337198/126‑2 (MQ=255)
                                 aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATcgt                                                                                                                                            >  8:337198/1‑125 (MQ=255)
                                 aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                                     <  3:55632/149‑1 (MQ=255)
                                 aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                                     >  7:60698/1‑149 (MQ=255)
                                     ggATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGaa                                                                                                                 <  8:321371/149‑1 (MQ=255)
                                               tCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTgg                                                                                                       <  6:338556/149‑1 (MQ=255)
                                                 aagATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGTCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGaa                                                                                                     <  8:18043/147‑1 (MQ=255)
                                                 aCGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGaa                                                                                                     <  6:67301/149‑1 (MQ=255)
                                                 aCGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGa                                                                                                      <  3:357723/148‑1 (MQ=255)
                                                 aCGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGa                                                                                                      <  7:242021/148‑1 (MQ=255)
                                                  cGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGCCCATACCGAACACATTCTgg                                                                                                       >  8:317592/1‑146 (MQ=255)
                                                   gATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATTAGCAGCGTCATGACCATACCGAACACATTCTgg                                                                                                       >  6:272548/1‑145 (MQ=255)
                                                                     gATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCt                                                                                                         >  7:284077/1‑125 (MQ=255)
                                                                         aaCCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGCTGGACCACCGGAtcttc                                                                             <  4:122063/149‑5 (MQ=255)
                                                                                  aaGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCgat                                                                    <  6:201337/149‑1 (MQ=255)
                                                                                    gagGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCgatga                                                                  >  3:278224/1‑149 (MQ=255)
                                                                                      ggACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATTAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATc                                                                <  5:272548/149‑1 (MQ=255)
                                                                                        aCCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcag                                                              <  7:198804/149‑1 (MQ=255)
                                                                                          cGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcagca                                                            <  8:60698/149‑1 (MQ=255)
                                                                                              ttGACCTGAATGATACTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCt                                                        >  6:127319/1‑149 (MQ=255)
                                                                                                gACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTg                                                       >  6:272537/1‑148 (MQ=255)
                                                                                                 aCCTGAATGATGCTGCCATGGCGCAAGACAACGCGCAGGATGGAATGCTTGAATTCGATCGCGATCAGCAGCGTCATGACGATACCGAACACTTTCTGGAAGATGCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCATGCTGaa                                                     <  1:386373/149‑1 (MQ=255)
                                                                                                 aCCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTg                                                       >  5:335616/1‑147 (MQ=255)
                                                                                                  ccTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGaa                                                     <  7:45996/148‑1 (MQ=255)
                                                                                                   cTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGa                                                   >  1:83602/1‑149 (MQ=255)
                                                                                                      aaTGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCg                                                                      >  8:102429/1‑127 (MQ=255)
                                                                                                      aaTGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCg                                                                      <  7:102429/127‑1 (MQ=255)
                                                                                                             ctgcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCTGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                                         <  1:271688/149‑1 (MQ=255)
                                                                                                               gcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                                         >  6:264232/1‑147 (MQ=255)
                                                                                                               gcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCAGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                                         >  4:375534/1‑147 (MQ=255)
                                                                                                                  aTGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATc                                              <  4:273696/139‑1 (MQ=255)
                                                                                                                  aTGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATc                                              >  3:273696/1‑139 (MQ=255)
                                                                                                                       gcAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTAc                                >  5:141398/1‑148 (MQ=255)
                                                                                                                        cAAGACAACGCGCCCGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCa                              <  3:293658/149‑1 (MQ=255)
                                                                                                                          aGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGc                            <  5:272537/149‑1 (MQ=255)
                                                                                                                                   gcACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcagc                                                             <  5:267099/107‑1 (MQ=255)
                                                                                                                                   gcACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcagc                                                             >  6:267099/1‑107 (MQ=255)
                                                                                                                                         tGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAacc             >  6:272482/1‑149 (MQ=255)
                                                                                                                                             aTGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAacca            >  5:197098/1‑146 (MQ=255)
                                                                                                                                             aTGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAacca            <  6:197098/146‑1 (MQ=255)
                                                                                                                                                    aaTTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACCGCGATCa   <  5:127319/148‑1 (MQ=255)
                                                                                                                                                    aaTTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACCGCGATCAc  >  6:34102/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACCGCGATCAC  >  NC_002947/3750592‑3750888

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