Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,568,744 C→A 11.8% V283F (GTC→TTC)  PP_2252 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,568,7440CA11.8% 84.1 / 3.8 34V283F (GTC→TTC) PP_2252hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/19);  new base A (2/2);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

CGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTG  >  NC_002947/2568604‑2568889
                                                                                                                                            |                                                                                                                                                 
cggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGACCACCCGCAGATGGTCGAACAACTGCgg                                                                                                                                          <  2:321653/150‑1 (MQ=255)
   cgcACTGGCCCGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGggtg                                                                                                                                       <  1:389777/150‑1 (MQ=255)
          gCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACgg                                                                                                                                <  2:750029/150‑1 (MQ=255)
               cAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCCTCTGCGGGTGGTCACGGTCATg                                                                                                                           <  2:348280/150‑1 (MQ=255)
                     gTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcat                                                                                                                     <  2:112967/150‑1 (MQ=255)
                      tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatc                                                                                                                    >  1:346661/1‑150 (MQ=255)
                        cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGTTGGTCACGGTCATGGGTCATCAc                                                                                                                  >  2:412265/1‑150 (MQ=255)
                        cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                                  >  1:168102/1‑150 (MQ=255)
                                      aCAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTg                                                                                                    >  1:342434/1‑150 (MQ=255)
                                            aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGTTGGTCACGTTCATGGGTCATCACCCCGTAACCGGCTGCGGCaa                                                                                              <  1:787843/150‑1 (MQ=255)
                                             aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAAt                                                                                             <  1:704560/150‑1 (MQ=255)
                                                 aCCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaa                                                                                         >  2:647974/1‑150 (MQ=255)
                                                                                ccATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCg                                                          <  2:346661/150‑1 (MQ=255)
                                                                                 cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGa                                                         <  1:785699/150‑1 (MQ=255)
                                                                                      aCTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCCGTTGTTCGACCCTCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCa                                                    <  1:235351/150‑1 (MQ=255)
                                                                                            gCCATTCCTGTTTCAACGCCACGTTCAGGATAACCCGCAGTTGTTCCACCACCTGCGGGTGGTCAATGTCATGGGTCATCACTCCGTAGCCGGCTCCGGCAATCGAAAACATCGCCTCTGCCATCTTGTCCTCTGCCGATGCCACGGCAc                                              <  2:89622/150‑1 (MQ=255)
                                                                                            gCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCAc                                              >  2:344088/1‑150 (MQ=255)
                                                                                               aTTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTgg                                           <  2:843015/150‑1 (MQ=255)
                                                                                                    tGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATc                                      <  2:780630/150‑1 (MQ=255)
                                                                                                            cGCCACGTTCAGGATCGCCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgg                              <  1:647974/150‑1 (MQ=255)
                                                                                                               cACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                           <  2:168102/150‑1 (MQ=255)
                                                                                                               cACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                           >  2:374266/1‑150 (MQ=255)
                                                                                                                            aCCCGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg                 <  1:522490/147‑1 (MQ=255)
                                                                                                                            aCCCGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg                 >  2:522490/1‑147 (MQ=255)
                                                                                                                             cccGCAGTTGCTTGACCCTCTGCGTTTGATCAAGGTCACGGGTCATCCCCCCGTAACCGCCTGCGGCAATCGAAATCATCCCCCCTGCCCCCTCGCCCCCTCCCCCTGCCACGGCACTCGCGATCCCCCGTGGTCCCAGGGTGTCGCCaa             <  2:32013/150‑1 (MQ=255)
                                                                                                                              ccGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGTCAAg            >  1:591066/1‑150 (MQ=255)
                                                                                                                                    ttgttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                        <  2:586611/100‑1 (MQ=255)
                                                                                                                                    ttgttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                        >  1:586611/1‑100 (MQ=255)
                                                                                                                                      gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTCGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGc    >  2:599371/1‑150 (MQ=255)
                                                                                                                                       ttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCGCCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCATCTACACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCt   >  1:692427/1‑150 (MQ=255)
                                                                                                                                       ttCGACCATCTGCGGGTGGTCACGGTCATGCGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCGTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCt   <  2:518919/150‑1 (MQ=255)
                                                                                                                                        tCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTg  <  2:465329/150‑1 (MQ=255)
                                                                                                                                           aCCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCg                                                                                                  >  2:561708/1‑51 (MQ=255)
                                                                                                                                           aCCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCg                                                                                                  <  1:561708/51‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTG  >  NC_002947/2568604‑2568889

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