Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,665,425 C→T 8.3% A97V (GCA→GTA)  PP_3229 → aliphatic sulfonate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,665,4250CT 8.3% 129.2 / 3.1 48A97V (GCA→GTA) PP_3229aliphatic sulfonate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

TCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCAAAGCCGCTGGCAGCCACTATCTGT  >  NC_002947/3665279‑3665574
                                                                                                                                                  |                                                                                                                                                     
tctcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAgg                                                                                                                                                    >  2:774993/1‑150 (MQ=255)
tctcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAgg                                                                                                                                                    >  2:345314/1‑150 (MQ=255)
tctcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAgg                                                                                                                                                    >  2:345110/1‑150 (MQ=255)
  tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgcg                                                                                                                                                  >  2:1018780/1‑150 (MQ=255)
    gAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGcc                                                                                                                                                <  2:763497/150‑1 (MQ=255)
      aCAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCaa                                                                                                                                              <  2:75399/150‑1 (MQ=255)
           gcCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGa                                                                                                                                         <  1:1018780/150‑1 (MQ=255)
            ccTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGAc                                                                                                                                        <  2:1025388/150‑1 (MQ=255)
                   gCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAGCCTTGGAACTGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTc                                                                                                                                 <  1:134708/150‑1 (MQ=255)
                            aTTCGTGTGACCTGGCACTAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGaa                                                                                                                        >  2:565899/1‑150 (MQ=255)
                              tCGCTTGACGTGGCACGAGTTCCCCCGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCATGCCTGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCCCCCGGGAAGc                                                                                                                      <  2:278466/150‑1 (MQ=255)
                                gcgTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCAc                                                                                                                    >  1:741562/1‑150 (MQ=255)
                                  gTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgc                                                                                                                                                   <  1:185296/117‑1 (MQ=255)
                                  gTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgc                                                                                                                                                   >  2:185296/1‑117 (MQ=255)
                                       cTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCAGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTCAGTCCTGGAAGCACCTTCGcc                                                                                                             >  1:335497/1‑150 (MQ=255)
                                       cTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGAAGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGcc                                                                                                             <  1:917272/150‑1 (MQ=255)
                                           cACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCTGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGAcc                                                                                                         <  2:263641/150‑1 (MQ=255)
                                               aGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGccc                                                                                                     >  1:559348/1‑150 (MQ=255)
                                                 ttCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGTAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAg                                                                                                   >  2:1024178/1‑150 (MQ=255)
                                                 ttCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGTAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAg                                                                                                   >  2:1024203/1‑150 (MQ=255)
                                                            ccTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGCTCCTGGt                                                                                        >  2:45132/1‑150 (MQ=255)
                                                            ccTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGt                                                                                        <  1:345314/150‑1 (MQ=255)
                                                            ccTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGt                                                                                        <  1:345110/150‑1 (MQ=255)
                                                                  gctgctGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCcggc                                                                                  >  1:276922/1‑150 (MQ=255)
                                                                   ctgctgGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCcggcg                                                                                 <  2:559348/150‑1 (MQ=255)
                                                                     gctgGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGa                                                                               <  2:177162/150‑1 (MQ=255)
                                                                       tgGAGGCCCTGAGCCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACt                                                                             >  2:960874/1‑150 (MQ=255)
                                                                       tgGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACt                                                                             >  2:369777/1‑150 (MQ=255)
                                                                            gCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCg                                                                        <  2:276922/150‑1 (MQ=255)
                                                                                     cTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGTAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGAcc                                                               <  1:1024178/150‑1 (MQ=255)
                                                                                     cTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGTAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGAcc                                                               <  1:1024203/150‑1 (MQ=255)
                                                                                              gTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGAc                                                      <  1:905355/150‑1 (MQ=255)
                                                                                              gTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGAc                                                      <  1:484076/150‑1 (MQ=255)
                                                                                               tCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGAcc                                                     >  2:1001781/1‑150 (MQ=255)
                                                                                                cGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCt                                                    >  2:920541/1‑150 (MQ=255)
                                                                                                     tctcGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAg                                               >  1:377286/1‑150 (MQ=255)
                                                                                                         ggcgtATGTTGCTGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAAACCCTGGCCGACCTCAAAGGCaa                                           <  2:550573/150‑1 (MQ=255)
                                                                                                           cggACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGc                                         >  1:625684/1‑150 (MQ=255)
                                                                                                                gTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGtgtg                                    >  1:99369/1‑150 (MQ=255)
                                                                                                                   gCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCg                                 <  2:741562/150‑1 (MQ=255)
                                                                                                                    ccGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCgg                                <  2:99369/150‑1 (MQ=255)
                                                                                                                       acacCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGa                             <  1:1001781/150‑1 (MQ=255)
                                                                                                                           cGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGACTGGGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCaa                         >  1:619881/1‑150 (MQ=255)
                                                                                                                                    gtTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCAAAGCCGCTgg                >  1:281681/1‑150 (MQ=255)
                                                                                                                                     ttCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACTGCTCAGGCGATCCTGGTGCCTGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCAAAGCCGCTGGc               >  2:143209/1‑150 (MQ=255)
                                                                                                                                        aCTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCAGACCTCAAAGGCAGGCGTGTGGCGGTGACCAAAGCCGCTGGCAGc            <  1:140538/150‑1 (MQ=255)
                                                                                                                                         cTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCAAAGCCGCTGGCAGcc           >  1:1066342/1‑150 (MQ=255)
                                                                                                                                          tCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGAACGCGCGTGCAGCTGTCCAGAGACACTgccagcca          >  2:1087164/1‑150 (MQ=255)
                                                                                                                                                  cAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGGGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCAAAGCCGCTGGCAGCCACTATCTGt  <  1:565899/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGCCGGCGGACTCACCGCTCAAGACCCTGGCCGACCTCAAAGGCAAGCGTGTGGCGGTGACCAAAGCCGCTGGCAGCCACTATCTGT  >  NC_002947/3665279‑3665574

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