Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 844,722 G→A 7.5% L41F (CTT→TTT)  PP_0727 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947844,7220GA 7.5% 150.8 / 2.2 53L41F (CTT→TTT) PP_0727hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

GCAAGATCAAGCAAGACCCTCTCAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCG  >  NC_002947/844577‑844866
                                                                                                                                                 |                                                                                                                                                
gCAAGATCAAGCAAGACCCTCTCAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTa                                                                                                                                               >  5:12137/1‑149 (MQ=255)
      tCAAGCAAGACCCTCTCAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCgg                                                                                                                                         >  2:224914/1‑149 (MQ=255)
            aagACCCTCTCAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAg                                                                                                                                   >  7:76192/1‑149 (MQ=255)
                   tctcAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACAt                                                                                                                            >  8:226554/1‑149 (MQ=255)
                     tcAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTg                                                                                                                          >  5:115250/1‑149 (MQ=255)
                       aGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTg                                                                                                                          <  8:76192/147‑1 (MQ=255)
                       aGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTg                                                                                                                          <  1:224914/147‑1 (MQ=255)
                       aGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGt                                                                                                                        >  1:239434/1‑149 (MQ=255)
                                   tgtgtCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCaa                                                                                                            <  2:239434/149‑1 (MQ=255)
                                   tgtgtCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCaa                                                                                                            >  8:179719/1‑149 (MQ=255)
                                     tgtCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAAcc                                                                                                          <  6:88399/149‑1 (MQ=255)
                                       tCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCGTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAAcc                                                                                                          >  3:234699/1‑147 (MQ=255)
                                        cGGCTGTCCCACTCAGGGCCTGGCAACTACCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATGCAACCCAg                                                                                                       >  8:14958/1‑149 (MQ=255)
                                          gCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAgcg                                                                                                     <  4:302244/149‑1 (MQ=255)
                                                ccACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTc                                                                                               <  5:7286/149‑1 (MQ=255)
                                                      aGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAgcgc                                                                                         >  4:19985/1‑149 (MQ=255)
                                                      aGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAgcgc                                                                                         <  4:234699/149‑1 (MQ=255)
                                                            tGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCg                                                                                   >  4:26483/1‑149 (MQ=255)
                                                              gCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGc                                                                                  >  7:113911/1‑148 (MQ=255)
                                                                     cctccACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTg                                                                          >  2:138611/1‑149 (MQ=255)
                                                                      ctccACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTgg                                                                         >  3:146056/1‑149 (MQ=255)
                                                                      ctccACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTgg                                                                         <  7:14958/149‑1 (MQ=255)
                                                                      ctccACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAAGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTgg                                                                         >  2:158098/1‑149 (MQ=255)
                                                                      ctccACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAAGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATACCTgg                                                                         >  2:162916/1‑149 (MQ=255)
                                                                             ggTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCTCATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCgg                                                                  <  4:146056/149‑1 (MQ=255)
                                                                               tAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGTGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCTCGTCCCGCCATTCCTGGCGTTCGGCg                                                                <  1:271293/149‑1 (MQ=255)
                                                                               tAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCg                                                                >  6:175575/1‑149 (MQ=255)
                                                                                aaGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGt                                                               >  4:67905/1‑149 (MQ=255)
                                                                                 accAGGAGGCCCTGAGGGCGCAATGATTGGACTACCGCCATGCTGACTCTTGGATCGTTTGGAAGGGACAGCGGGGGCAGCCCAAATTGGTCTACGTATTCAACCGAGCTCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGtt                                                              <  7:260528/147‑1 (MQ=255)
                                                                                      gaggCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACt                                                         >  3:24663/1‑149 (MQ=255)
                                                                                        ggCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGc                                                       <  5:175575/149‑1 (MQ=255)
                                                                                           ccTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCaa                                                                                                            <  8:164110/93‑1 (MQ=255)
                                                                                            cTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAgc                                                                                           >  4:85639/1‑109 (MQ=255)
                                                                                            cTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAgc                                                                                           <  3:85639/109‑1 (MQ=255)
                                                                                            cTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGa                                                   <  1:138611/149‑1 (MQ=255)
                                                                                            cTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGa                                                   <  3:26483/149‑1 (MQ=255)
                                                                                            cTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGa                                                   <  8:113911/149‑1 (MQ=255)
                                                                                                gggCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCAt                                               <  7:179719/149‑1 (MQ=255)
                                                                                                       aTGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGtt                                                                     >  2:85971/1‑120 (MQ=255)
                                                                                                                   aCGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCAtt                                                                              >  6:81768/1‑99 (MQ=255)
                                                                                                                   aCGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCaaaa                              <  3:19985/147‑1 (MQ=255)
                                                                                                                     gCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCAt                          <  4:24663/149‑1 (MQ=255)
                                                                                                                     gCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCAt                          <  7:226554/149‑1 (MQ=255)
                                                                                                                       cATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGc                        <  6:115250/149‑1 (MQ=255)
                                                                                                                       cATGCTGACTCTTGGATCGTTTGGAAAGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGc                        <  1:158098/149‑1 (MQ=255)
                                                                                                                       cATGCTGACTCTTGGATCGTTTGGAAAGTACAGCGGGGTCAGCCCACATTGGGCTAATAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGc                        <  1:162916/149‑1 (MQ=255)
                                                                                                                                     gATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGtata           <  3:67905/148‑1 (MQ=255)
                                                                                                                                     gATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGtata           <  6:12137/148‑1 (MQ=255)
                                                                                                                                      aTCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATAcg         >  8:107960/1‑149 (MQ=255)
                                                                                                                                             ggAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCg  >  6:291876/1‑149 (MQ=255)
                                                                                                                                             ggAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCg  >  6:251398/1‑149 (MQ=255)
                                                                                                                                               aaGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATg                         <  7:185906/124‑1 (MQ=255)
                                                                                                                                               aaGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATg                         >  8:185906/1‑124 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAAGATCAAGCAAGACCCTCTCAGTGACGTCGCGTGTGTCGGCTGTCCCACTCAGGGCCTGGCAACTCCCTCCACCGGTAAGCAGGAGGCCCTGAGGGCGCAATGATTGGACTAACGCCATGCTGACTCTTGGATCGTTTGGAAGGTACAGCGGGGTCAGCCCACATTGGTCTAAGAATTCAACCCAGCGCTGGTCGAGCGCGTCCCGCCATTCCTGGCGTTCGGCGTTGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCG  >  NC_002947/844577‑844866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑