Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 334,605 A→T 13.0% L28* (TTG→TAG) ‡ PP_0276 ← Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947334,6050AT13.0% 141.8 / 13.7 54L28* (TTG→TAG) ‡PP_0276Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base T (4/3);  total (32/22)
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.21e-01

AGCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTC  >  NC_002947/334458‑334752
                                                                                                                                                   |                                                                                                                                                   
aGCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAg                                                                                                                                                   >  1:350426/1‑150 (MQ=255)
 gCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGt                                                                                                                                                  >  2:1254130/1‑150 (MQ=255)
 gCATGTCCATCATGGCGTCCAAGGTGAATTTATCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGt                                                                                                                                                  >  1:387321/1‑150 (MQ=255)
            atGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTGGAGGCCCCTTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCATTTTCatc                                                                                                                                       <  2:833911/150‑3 (MQ=255)
             tGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAgg                                                                                                                                      >  2:1920116/1‑150 (MQ=255)
                 gTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCAt                                                                                                                                  >  1:1313639/1‑150 (MQ=255)
                       ggTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCt                                                                                                                            >  2:1541014/1‑150 (MQ=255)
                        gTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCAttt                                                                                                                                            <  2:1467384/133‑1 (MQ=255)
                        gTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCAttt                                                                                                                                            >  1:1467384/1‑133 (MQ=255)
                             tttCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGt                                                                                                                      >  1:266732/1‑150 (MQ=255)
                                   gATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGc                                                                                                                <  1:561297/150‑1 (MQ=255)
                                     tCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGcgt                                                                                                              >  2:199244/1‑150 (MQ=255)
                                     tCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGcgt                                                                                                              <  1:1254130/150‑1 (MQ=255)
                                      cTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGcgtc                                                                                                             <  2:826436/150‑1 (MQ=255)
                                          ccATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCgtgg                                                                                                         <  2:350426/150‑1 (MQ=255)
                                           cATTTGCAAGCTCGGATAGGCGGGCTCGAATGACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCgtggg                                                                                                        >  1:990447/1‑150 (MQ=255)
                                                      tCGGATACGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                             <  1:1459650/150‑1 (MQ=255)
                                                       cGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAg                                                                                            <  1:650737/150‑1 (MQ=255)
                                                        ggATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGt                                                                                                  >  2:1735426/1‑143 (MQ=255)
                                                        ggATATGCGGGATTGGGTTACGTGCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGt                                                                                                  <  1:1735426/143‑1 (MQ=255)
                                                                         ttACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCAGCGTCGTTGGTGGCGTGGAAAGGGCTTTCGTGGAATTGaa                                                                          <  1:150929/150‑1 (MQ=255)
                                                                                  gCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCg                                                                 >  2:165479/1‑150 (MQ=255)
                                                                                          gCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCt                                                         >  2:969236/1‑150 (MQ=255)
                                                                                                        tttGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCt                                           >  1:1836105/1‑150 (MQ=255)
                                                                                                        tttGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCAATCGGCTTGCCTCGGCAGAAATATCt                                           >  1:1836478/1‑150 (MQ=255)
                                                                                                           gAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTaaa                                        >  1:1461021/1‑150 (MQ=255)
                                                                                                               cccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAAttt                                    >  1:746066/1‑150 (MQ=255)
                                                                                                                ccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGTTGGCGTTGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTg                                   >  2:837322/1‑150 (MQ=255)
                                                                                                                ccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCTTGGAATTGAAGCGCCATCGGCTTGCCCCGGCAGAAATATTTAAAAtttt                                   >  2:642386/1‑149 (MQ=255)
                                                                                                                 cATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                             <  1:626761/91‑1 (MQ=255)
                                                                                                                 cATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                             >  2:626761/1‑91 (MQ=255)
                                                                                                                  aTTTTTTTGATGCAGTCAGTGATGAAAATGGGCAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCg                                 >  2:551519/1‑150 (MQ=255)
                                                                                                                  aTTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAgc                                                                        <  2:237094/111‑1 (MQ=255)
                                                                                                                  aTTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAgc                                                                        >  1:237094/1‑111 (MQ=255)
                                                                                                                     tttttGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCGCGGCAGAAATATCTAAAATTTGCGATa                              >  1:1172449/1‑150 (MQ=255)
                                                                                                                         tGATGCAGTCAGTGATGACAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATc                          <  2:827655/150‑1 (MQ=255)
                                                                                                                                    gtGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAt               <  2:1836105/150‑1 (MQ=255)
                                                                                                                                    gtGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAt               <  2:1836478/150‑1 (MQ=255)
                                                                                                                                     tgatgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAg                                                                                            <  1:1846325/72‑1 (MQ=255)
                                                                                                                                     tgatgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAg                                                                                            >  2:1846325/1‑72 (MQ=255)
                                                                                                                                      gatgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGt                                                                                                  <  1:1117958/65‑1 (MQ=255)
                                                                                                                                      gatgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGt                                                                                                  >  2:1117958/1‑65 (MQ=255)
                                                                                                                                      cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCCGAAatat                                             >  2:1117642/3‑118 (MQ=255)
                                                                                                                                      cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAatat                                             <  1:1117642/116‑1 (MQ=255)
                                                                                                                                       ttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttgg            >  1:1840034/2‑150 (MQ=255)
                                                                                                                                       ttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttgg            >  2:829660/2‑150 (MQ=255)
                                                                                                                                        tgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttggt           >  1:1728834/1‑150 (MQ=255)
                                                                                                                                         gaAAATGGACAAGGCCATTTTCAAGGCCATTTTGCTGGTGGTATCTGCGTCGTGGGTGGCGTGGAGAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCGTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttggtt          <  2:655723/150‑1 (MQ=255)
                                                                                                                                            aaTGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGGAGGAAGGGCTTTTCTTTAAGTTAATCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGc       >  1:1178293/1‑150 (MQ=255)
                                                                                                                                                gACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGt   <  1:199244/150‑1 (MQ=255)
                                                                                                                                                gACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGt   <  1:1920116/150‑1 (MQ=255)
                                                                                                                                                 aCAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATATATCGTAGTGTGGATTGGTTGGCGAGTc  >  1:870/1‑150 (MQ=255)
                                                                                                                                                 aCAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAt               <  2:183940/137‑1 (MQ=255)
                                                                                                                                                 aCAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAt               >  1:183940/1‑137 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTC  >  NC_002947/334458‑334752

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: