Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,183,617 T→C 9.1% A465A (GCT→GCC)  PP_3680 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,183,6170TC 9.1% 122.4 / 2.7 44A465A (GCT→GCC) PP_3680hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/17);  new base C (2/2);  total (25/19)
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.60e-01

GTTCGAGAAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTC  >  NC_002947/4183474‑4183764
                                                                                                                                               |                                                                                                                                                   
gttcgAGAAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCAc                                                                                                                                                 >  1:260230/1‑148 (MQ=255)
    gagaAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAggg                                                                                                                                            >  2:163321/1‑149 (MQ=255)
        aGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgc                                                                                                                                        <  2:260230/149‑1 (MQ=255)
          ccGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgcgc                                                                                                                                      >  3:111075/1‑149 (MQ=255)
          ccGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgcgc                                                                                                                                      >  3:106494/1‑149 (MQ=255)
            gTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCCCCCCTAGTGGAgcgatc                                                                                                                                    >  1:313678/1‑146 (MQ=255)
            gTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCCCCCCTAGTGGAgcgatc                                                                                                                                    <  2:313678/149‑4 (MQ=255)
                  gCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCaa                                                                                                                              >  4:250736/1‑149 (MQ=255)
                          cGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCt                                                                                                                      <  3:250736/149‑1 (MQ=255)
                              ctctACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTggc                                                                                                                   <  8:62271/148‑1 (MQ=255)
                                       gCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGaa                                                                                                         >  3:88572/1‑149 (MQ=255)
                                                 tggtACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGt                                                                                               <  1:163321/149‑1 (MQ=255)
                                                    tACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTc                                                                                                                                <  1:145440/113‑1 (MQ=255)
                                                    tACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTc                                                                                                                                >  2:145440/1‑113 (MQ=255)
                                                               ggggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGgctg                                                                                                              >  4:82905/1‑120 (MQ=255)
                                                               ggggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGgctg                                                                                                              <  3:82905/120‑1 (MQ=255)
                                                                 ggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGAccc                                                                               >  1:89639/1‑149 (MQ=255)
                                                                       tGAACAAAACCTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATAGCTACAATAGGCGGGCGCCATGTCAATCACATCTGTCGGCTGCTGAACGCGCTGGGTAAACCTCATATCACCCTGCTCg                                                                         >  4:264386/1‑149 (MQ=255)
                                                                        gaagaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTGat                                                                                      >  8:37657/1‑135 (MQ=255)
                                                                         aagaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGAc                                                                       >  5:85467/1‑149 (MQ=255)
                                                                          agaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGAcc                                                                      <  5:275836/149‑1 (MQ=255)
                                                                                         tCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCAc                                                                                                                                                 >  4:76432/1‑59 (MQ=255)
                                                                                         tCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCAc                                                                                                                                                 <  3:76432/59‑1 (MQ=255)
                                                                                              tCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                  <  4:106494/149‑1 (MQ=255)
                                                                                              tCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                  <  4:111075/149‑1 (MQ=255)
                                                                                                 ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                  <  3:119649/146‑1 (MQ=255)
                                                                                                 ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                  >  4:119649/1‑146 (MQ=255)
                                                                                                 ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGt                                               >  1:139210/1‑149 (MQ=255)
                                                                                                      gCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTg                                          >  2:268010/1‑149 (MQ=255)
                                                                                                            ggAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCg                                    >  7:149298/1‑149 (MQ=255)
                                                                                                                 tCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACCCCTCATATGACCCTGTTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGtgtg                                >  8:189026/1‑148 (MQ=255)
                                                                                                                 tCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGCCCCTGCTCGACCTAGACGTCGGGCGCCATCAAGGTGGCTGGGGGCGtgtg                                >  4:10908/1‑148 (MQ=255)
                                                                                                                  cccGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGGCAATCACTTCTGGCGGCTGCGGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGGGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGaa                              <  2:188066/149‑1 (MQ=255)
                                                                                                                   ccGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAAc                             <  4:88572/149‑1 (MQ=255)
                                                                                                                         aCGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGtt                       <  3:264386/149‑1 (MQ=255)
                                                                                                                          cGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                      <  1:268010/149‑1 (MQ=255)
                                                                                                                                 gCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACt                 >  2:315628/1‑147 (MQ=255)
                                                                                                                                    gTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCt            <  2:139210/149‑1 (MQ=255)
                                                                                                                                       catggcGCCCCCCTAGTGGAGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                      >  5:36847/6‑136 (MQ=255)
                                                                                                                                       catggcGCCCCCCTAGTGGAGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                      <  6:36847/131‑1 (MQ=255)
                                                                                                                                         gATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATg       >  7:99220/1‑149 (MQ=255)
                                                                                                                                             gCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTAtt   >  3:155777/1‑149 (MQ=255)
                                                                                                                                             gCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTAtt   <  6:85467/149‑1 (MQ=255)
                                                                                                                                              cTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACATCATATGACCCTGCTCGACCTAGACGTGGGGCGACATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTc  >  4:282598/1‑149 (MQ=255)
                                                                                                                                               tCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTc  >  6:153661/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GTTCGAGAAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTC  >  NC_002947/4183474‑4183764

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