Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 198,777 T→A 6.7% L1428Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947198,7770TA 6.7% 179.2 / 2.9 59L1428Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (2/2);  total (30/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 = 9.78e-01

GTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGG  >  NC_002947/198629‑198924
                                                                                                                                                    |                                                                                                                                                   
gTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTg                                                                                                                                                    <  1:578906/150‑1 (MQ=255)
   aGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTc                                                                                                                                                 <  1:252028/150‑1 (MQ=255)
     cACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGt                                                                                                                                               >  2:978321/1‑150 (MQ=255)
        cGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCAGGTCGTGag                                                                                                                                            <  1:553826/150‑2 (MQ=255)
          aGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGTCACGACCAGGTCGTGAgcca                                                                                                                                         <  2:991254/150‑2 (MQ=255)
           ggTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                         <  1:978321/150‑1 (MQ=255)
            gTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTg                                                                                                                                        >  1:198555/1‑150 (MQ=255)
             tCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                       <  2:125426/150‑1 (MQ=255)
                       gCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGc                                                                                                                             >  1:179118/1‑150 (MQ=255)
                          aGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                          >  2:146068/1‑150 (MQ=255)
                           ggTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                         >  1:486540/1‑150 (MQ=255)
                               cGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                      <  1:987364/133‑1 (MQ=255)
                               cGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                      >  2:987364/1‑133 (MQ=255)
                                gCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                    >  1:483838/1‑150 (MQ=255)
                                gCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                    <  2:198555/150‑1 (MQ=255)
                                   aGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                                 >  1:790532/1‑150 (MQ=255)
                                   aGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                                 >  1:826364/1‑150 (MQ=255)
                                       aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                             >  2:953609/1‑150 (MQ=255)
                                       aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                             <  1:800292/150‑1 (MQ=255)
                                        aCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAgg                                                                                                            <  1:953609/150‑1 (MQ=255)
                                             ggCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTc                                                                                                                      <  1:519872/135‑1 (MQ=255)
                                             ggCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTc                                                                                                                      >  2:519872/1‑135 (MQ=255)
                                                gATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGaa                                                                                                    >  1:1117120/1‑150 (MQ=255)
                                                  tatCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAAc                                                                                                  <  2:1117120/150‑1 (MQ=255)
                                                     acacAAAGTCTCCATCTCGTCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAg                                                                                               <  2:790532/147‑1 (MQ=17)
                                                      gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                              >  1:2512/1‑150 (MQ=255)
                                                         aaGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAGCGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                      >  1:1361040/1‑107 (MQ=255)
                                                         aaGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                      <  2:1361040/107‑1 (MQ=255)
                                                         aaGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCa                                                                                           <  1:146068/150‑1 (MQ=255)
                                                                  aTCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACt                                                                                  <  2:179118/150‑1 (MQ=255)
                                                                                 aCTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCac                                                                               >  2:163388/1‑138 (MQ=255)
                                                                                 aCTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCac                                                                               <  1:163388/138‑1 (MQ=255)
                                                                                   ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGa                                                                 >  2:694693/1‑150 (MQ=255)
                                                                                       gTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGAc                                                             <  1:694693/150‑1 (MQ=255)
                                                                                       gTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGAc                                                             <  2:483838/150‑1 (MQ=255)
                                                                                            cGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTa                                                        >  2:1028451/1‑150 (MQ=255)
                                                                                            cGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTa                                                        >  2:1090179/1‑150 (MQ=255)
                                                                                               gaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCa                                                     >  1:913580/1‑150 (MQ=255)
                                                                                                  cGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGa                                                  <  2:859319/150‑1 (MQ=255)
                                                                                                    tCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATg                                                >  2:581397/1‑150 (MQ=255)
                                                                                                       aaCCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGAGCCATGACTTCTATCAGGATGCgg                                             <  2:486540/150‑1 (MQ=255)
                                                                                                              cTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTc                                                                                                                           <  2:1196794/65‑1 (MQ=255)
                                                                                                              cTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTc                                                                                                                           >  1:1196794/1‑65 (MQ=255)
                                                                                                                tCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCa                                    <  1:1090179/150‑1 (MQ=255)
                                                                                                                  aCTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGc                                                                                                           >  2:1060968/1‑77 (MQ=255)
                                                                                                                  aCTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGc                                                                                                           <  1:1060968/77‑1 (MQ=255)
                                                                                                                     gtgCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTg                               <  1:581397/150‑1 (MQ=255)
                                                                                                                           aTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAAc                                                                                                  <  2:1413993/77‑1 (MQ=255)
                                                                                                                           aTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAAc                                                                                                  >  1:1413993/1‑77 (MQ=255)
                                                                                                                                 acCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcg                   >  2:1275040/1‑150 (MQ=255)
                                                                                                                                       cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTg             <  2:826364/150‑1 (MQ=255)
                                                                                                                                        tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAg                                                    >  1:1347280/1‑110 (MQ=255)
                                                                                                                                        tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAg                                                    <  2:1347280/110‑1 (MQ=255)
                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcgg                  <  2:30578/140‑1 (MQ=255)
                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcgg                  >  1:30578/1‑140 (MQ=255)
                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGc        >  1:6211/1‑150 (MQ=255)
                                                                                                                                            tCACGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGc        >  2:1263277/1‑150 (MQ=255)
                                                                                                                                            tCACGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCAGCAACTCGGCGGTGGACgg        >  1:840650/1‑149 (MQ=255)
                                                                                                                                                  ccTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTgg  >  2:1094990/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGG  >  NC_002947/198629‑198924

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