Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,383,254 A→T 7.7% L16Q (CTG→CAG)  PP_1204 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,383,2540AT 7.7% 137.3 / 3.4 52L16Q (CTG→CAG) PP_1204lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base T (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

ACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTA  >  NC_002947/1383107‑1383398
                                                                                                                                                    |                                                                                                                                                
acacGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAg                                                                                                                                                 >  3:108514/1‑149 (MQ=255)
    gTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcg                                                                                                                                             >  4:244542/1‑149 (MQ=255)
             gAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGccccc                                                                                                                                     >  8:23079/1‑148 (MQ=255)
             gAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGcccc                                                                                                                                      >  1:79002/1‑147 (MQ=255)
             gAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGcccc                                                                                                                                      <  2:79002/147‑1 (MQ=255)
              aaGCTCCCTCACCATCGGCTCGCCGACATCGCTCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCaa                                                                                                                                   <  6:266525/149‑1 (MQ=255)
               aGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                        >  4:271488/1‑143 (MQ=255)
               aGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                        <  3:271488/143‑1 (MQ=255)
                 cTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTGGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAggg                                                                                                                                <  4:68262/149‑1 (MQ=255)
                 cTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAggg                                                                                                                                <  5:299316/149‑1 (MQ=255)
                        aCCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGTACATCTGGTTG‑GCGCAGCCGCTTACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGt                                                                                                                         <  7:23079/149‑1 (MQ=255)
                         ccTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGt                                                                                                                         >  4:147465/1‑148 (MQ=255)
                         ccTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTa                                                                                                                        >  4:117096/1‑149 (MQ=255)
                         ccTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCAGTAATGTTAGGTTGATAGAACCCCCTGACGAGCAGGGGCCCCCCCAATGGGGCGCCCCCAAGGGTAAAgg                                                                                                                          >  1:111951/1‑148 (MQ=255)
                              ggCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCg                                                                                                                   <  1:226558/149‑1 (MQ=255)
                               gCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTAATTCGc                                                                                                                  >  1:151801/1‑149 (MQ=255)
                                   gCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTg                                                                                                               >  1:17688/1‑148 (MQ=255)
                                   gCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGGCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCtt                                                                                                                <  5:155784/147‑1 (MQ=255)
                                     cGACATCGATGTGCTGTGTTTGCTCGTGCAGTTGGGGCTCGACGGGCTCCTCATGGTCGCTTCGCGGGGACATCTGGTTG‑GCGAAGCCGCTGACTAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTAc                                                                                                            <  4:37629/149‑1 (MQ=255)
                                     cGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACACTGCGGCGCCCCCAAGGGGAAAGGTACTTCGCTTGTa                                                                                                             >  5:45571/1‑148 (MQ=255)
                                                  cAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCAc                                                                                                    <  3:261804/144‑1 (MQ=255)
                                                      gTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCAtt                                                                                           <  6:304099/149‑1 (MQ=255)
                                                       tATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTg                                                                                          <  2:94593/149‑1 (MQ=255)
                                                              gAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTAt                                                                                     >  8:136763/1‑147 (MQ=255)
                                                                      tGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGaa                                                                           >  4:214279/1‑149 (MQ=255)
                                                                           tCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATc                                                                                    <  5:83719/135‑1 (MQ=255)
                                                                           tCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATc                                                                                    >  6:83719/1‑135 (MQ=255)
                                                                                    tcctcATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAACGTACTTAGCTTGTACATGATCACACGGGCATAGGGTATCAGCGGCGAACGCGGGCCTggcg                                                              >  6:271078/1‑145 (MQ=255)
                                                                                         aTGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAg                                                          <  3:214279/147‑1 (MQ=255)
                                                                                           gggcgctgcgctgGTACATCTGGTGG‑GCGCAGCCGCGGACGAGATGGGCGCCGCACAGGGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTc                                                       <  6:45571/145‑1 (MQ=255)
                                                                                            gTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAg                                                     >  3:231985/1‑149 (MQ=255)
                                                                                               gctgcgctgGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCa                                                   <  3:244542/148‑1 (MQ=255)
                                                                                               gctgcgctgGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCAc                                                  >  1:130767/1‑149 (MQ=255)
                                                                                                   cgctgGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACtt                                                <  3:147465/147‑1 (MQ=255)
                                                                                                   cgctgGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACtt                                                <  3:117096/147‑1 (MQ=255)
                                                                                                     cggggACATCGGGTTG‑GCGCTGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCACGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAtt                                            <  4:37695/147‑1 (MQ=255)
                                                                                                        ggACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAg                                         >  6:5717/1‑149 (MQ=255)
                                                                                                         gACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGc                                        <  2:130767/149‑1 (MQ=255)
                                                                                                          aCATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCg                                       <  8:31938/149‑1 (MQ=255)
                                                                                                           cATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCgg                                      <  4:108514/149‑1 (MQ=255)
                                                                                                              cTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCaa                                   >  5:306051/1‑149 (MQ=255)
                                                                                                                    g‑gCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGc                                >  7:234937/1‑146 (MQ=255)
                                                                                                                    g‑gCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGc                                <  8:234937/146‑1 (MQ=255)
                                                                                                                      gcgcAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                            >  6:128597/1‑149 (MQ=255)
                                                                                                                          aGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGt                        <  5:5717/149‑1 (MQ=255)
                                                                                                                                 ctcgtgCAGGGCGCCGCACCGTGCGGGGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCt                  <  1:160095/143‑1 (MQ=255)
                                                                                                                                   ccttgggggCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcgg               >  5:182672/7‑149 (MQ=255)
                                                                                                                                            cgccgcACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCt      <  6:182672/149‑1 (MQ=255)
                                                                                                                                            cgccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCt      <  6:306051/149‑1 (MQ=255)
                                                                                                                                             gccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCt                     <  6:150446/133‑1 (MQ=255)
                                                                                                                                             gccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCt                     >  5:150446/1‑133 (MQ=255)
                                                                                                                                                gcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTa  <  2:111951/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTG‑GCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTA  >  NC_002947/1383107‑1383398

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